Implementation#

We want to solve the PDE-constrained shape optimization problem

\[ \min_{\Omega \subset \mathsf{D}} J(u) := \int_\Omega |\nabla u_\Omega - \nabla u_{ref}|^2 \, dx \]

subject to \((\Omega,u_\Omega)\) satisfying

\[ \int_\Omega \nabla u_\Omega \cdot \nabla v + u_\Omega v \, dx = \int_\Omega f v \, dx \quad \text{for all } v \in H_0^1(\Omega), \]

where \(\Omega \subset \mathbb{R}^2\) and \(u_d, f \in C^1(\mathbb{R}^2)\).

from ngsolve import *
from netgen.occ import *
from ngsolve.webgui import Draw
import matplotlib.pyplot as plt
from netgen.geom2d import SplineGeometry
import numpy as np
from ngsolve.nonlinearsolvers import Newton
maxH = 0.05
geo = OCCGeometry(Circle((0, 0), 0.9).Face(), dim=2).GenerateMesh(maxh=maxH)
mesh = Mesh(geo)

vol = 3e-2

a = 0.6
b = 0.4

ud = 1 - ((x-0.01)**2/a**2 + (y-0.01)**2/b**2)

f = (2/a**2 + 2/b**2) + ud

grad_f = CoefficientFunction( (f.Diff(x), f.Diff(y) ) )
grad_ud = CoefficientFunction( (ud.Diff(x), ud.Diff(y) ) )

#Draw(ud, mesh, "target_function")
fes = H1(mesh, order=2, dirichlet=".*")
u, v = fes.TnT()
gfu = GridFunction(fes)
scene_gfu = Draw (gfu, mesh, "state")

a = BilinearForm(fes, symmetric=True)
a += grad(u)*grad(v)*dx
a += u*v*dx

fstate = LinearForm(fes)
fstate += f*v*dx
def SolveStateEquation():
    rhs = gfu.vec.CreateVector()
    rhs.data = fstate.vec - a.mat * gfu.vec
    update = gfu.vec.CreateVector()
    update.data = a.mat.Inverse(fes.FreeDofs()) * rhs
    gfu.vec.data += update
a.Assemble()
fstate.Assemble()
SolveStateEquation()
scene_gfu.Redraw()

Adjoint equation#

The adjoint p satisfies

\[ \int_\Omega \nabla w \cdot \nabla p + w p \, dx = -2 \int_\Omega (\nabla u_\Omega - \nabla u_{ref}) w \quad \text{for all } w \in H_0^1(\Omega), \]
def Cost(u):
    cost = InnerProduct(grad(u) - grad_ud, grad(u) - grad_ud) * dx
    cost += InnerProduct(CF(vol), CF(1)) * dx

    return cost



p, w = fes.TnT()
gfp = GridFunction(fes)
scene_gfp = Draw (gfp, mesh, "adjoint")

fadjoint = LinearForm(fes)
#fadjoint += -1*Cost(gfu).Diff(gfu,w)
fadjoint += -2* InnerProduct(grad(gfu) - grad_ud, grad(w)) * dx
def SolveAdjointEquation():
    rhs = gfp.vec.CreateVector()
    rhs.data = fadjoint.vec - a.mat.T * gfp.vec
    update = gfp.vec.CreateVector()
    update.data = a.mat.Inverse(fes.FreeDofs()).T * rhs
    gfp.vec.data += update
fadjoint.Assemble()
SolveAdjointEquation()
scene_gfp.Redraw()

The shape derivative of the shape function

\[ \mathcal{J}(\Omega) := \int_\Omega |\nabla u - \nabla u_d|^2 \, dx \]

at \(\Omega\) in direction \(X\) is given by

\[\begin{split} \begin{aligned} DJ(\Omega)(X) &= \int_\Omega \operatorname{div}(X)|\nabla u - \nabla u_{\mathrm{ref}}|^2 - 2 \int_\Omega (\nabla u - \nabla u_{\mathrm{ref}}) \nabla^2 u_{\mathrm{ref}} \cdot X\\ &\quad + \int_\Omega \operatorname{div}(X) u p \\ &\quad + \int_\Omega \bigl(\operatorname{div}(X)I - \partial X - \partial X^\top\bigr) \nabla u \cdot \nabla p \\ &\quad - \int_\Omega (\nabla f \cdot X)\,p\\ & \quad - \int_\Omega \operatorname{div}(X)fp . \end{aligned} \end{split}\]

We now assemble this shape derivative in NGSolve as follows:

def hesse(u):
    return CoefficientFunction(
        (u.Diff(x).Diff(x), u.Diff(x).Diff(y),
         u.Diff(y).Diff(x), u.Diff(y).Diff(y)), 
        dims=(2,2)
    )

#VEC = H1(mesh, order=3, dirichlet=".*", dim=2)
VEC = H1(mesh, order=1, dim = 2)
X, PHI = VEC.TnT()

dJOmega = LinearForm(VEC)
dJOmega += InnerProduct((div(PHI)*(grad(gfu)- grad_ud)), grad(gfu)- grad_ud) * dx
dJOmega += -2* InnerProduct(hesse(ud)*PHI, grad(gfu)- grad_ud) * dx
dJOmega += div(PHI)* gfu *gfp *dx
dJOmega += InnerProduct(div(PHI)*grad(gfu), grad(gfp)) * dx
dJOmega += -InnerProduct(PHI.Deriv()*grad(gfu), grad(gfp)) * dx
dJOmega += -InnerProduct(grad(gfu), PHI.Deriv()*grad(gfp)) * dx
dJOmega += -div(PHI)* f * gfp * dx
dJOmega += -grad_f * PHI*gfp * dx
# added volumen constraint
dJOmega += -vol * div(PHI)* dx 
  • Next, we want to find a vector field \(X\) which yields a decrease of the objective function, i.e. we want to find a vector field \(X\) such that

    \[ d\mathcal{J}(\Omega; X) < 0. \]
  • This can be achieved by solving an auxiliary boundary value problem of the type

    \[ \text{Find } X \in H : \qquad b(X, \Phi) = d\mathcal{J}(\Omega; \Phi) \quad \text{for all } \Phi \in H, \]

    for a suitable Hilbert space \(H\) (e.g. \(H = H^1(\Omega)^2\)). Here, \(b(\cdot,\cdot) : H \times H \to \mathbb{R}\) is a positive definite bilinear form which we are free to choose. Then \(-X\) is a descent direction since

    \[ d\mathcal{J}(\Omega; -X) = - d\mathcal{J}(\Omega; X) = - b(X, X) < 0. \]
  • First we start with the bilinear form

    \[ b(X, Y) = \int_\Omega \nabla X : \nabla Y + X \cdot Y \, dx . \]
  • or

    \[ \tilde{b}(X,Y) = b(X, Y) + \gamma \int BX \cdot BY dx \]

    with $\( B := \begin{pmatrix} -\partial_x & \partial_y \\ \partial_y & \partial_x \end{pmatrix}.\)$

  • we can also choose X as $\( \min _{X \in H^1(\Omega)} \frac{1}{p} \int _\Omega |\nabla X|^p - DJ(\Omega)(X) \)$

    since this is a continous functional on \(H^1(\Omega)'\), it is equivalent to finding \(X \in H^1(\Omega)\) sucht that,

    \[ \int _\Omega |\nabla X | ^{p-2} \nabla X \cdot \nabla \Phi = DJ(\Omega)(\Phi) \qquad \forall \Phi \in H^1(\Omega) \]
def pose_deformation_equation(gamma=None, p = None): 
    b = BilinearForm(VEC)
    # standard case
    if not p:
        b += InnerProduct(grad(X),grad(PHI))*dx + InnerProduct(X,PHI)*dx

    if gamma:
        # standard case with gamma regularization
        print('Using gamma =', gamma)
        b = BilinearForm(VEC)
        b += (1.0/gamma)*InnerProduct(grad(X),grad(PHI))*dx + (1.0/gamma)*InnerProduct(X,PHI)*dx
    

        b += gamma * (PHI.Deriv()[0,0] - PHI.Deriv()[1,1])*(X.Deriv()[0,0] - X.Deriv()[1,1]) *dx
        b += gamma * (PHI.Deriv()[1,0] + PHI.Deriv()[0,1])*(X.Deriv()[1,0] + X.Deriv()[0,1]) *dx


        #print("expression tree\n")
        #print(b)
        return b
    
    elif p:
        
        # solving minimization problem
        print('Using p =', p)
        eps = 1E-5
    
        grad_norm = sqrt(InnerProduct(grad(X), grad(X)) + eps**2)
        

        c = BilinearForm(VEC)
        c += (grad_norm**(p-2)) * InnerProduct(grad(X), grad(PHI)) * dx
        c += InnerProduct(X,PHI)*dx
        
        # -DJ(Omega)(PHI)
        c += -InnerProduct((div(PHI)*(grad(gfu)- grad_ud)), grad(gfu)- grad_ud) * dx
        c += 2* InnerProduct(grad(gfu)- grad_ud, hesse(ud)*PHI) * dx
        c += -div(PHI)* gfu *gfp *dx
        c += -InnerProduct(div(PHI)*grad(gfu), grad(gfp)) * dx
        c += InnerProduct(PHI.Deriv()*grad(gfu), grad(gfp)) * dx
        c += InnerProduct(grad(gfu), PHI.Deriv()*grad(gfp)) * dx
        c += div(PHI)* f * gfp * dx
        c += grad_f*PHI*gfp * dx
        c += vol * div(PHI)*dx 
        
        # c += -div(PHI)* f * gfp * dx
        # c += -grad_f*PHI*gfp * dx
        
        return c

        
    #print("expression tree\n")
    #print(b)
    return b
   

def SolveDeformationEquation(linear=True):
    if linear:
        #print("solving linear deformation equation")
        rhs = gfX.vec.CreateVector()
        rhs.data = dJOmega.vec - b.mat * gfX.vec
        update = gfX.vec.CreateVector()
        update.data = b.mat.Inverse(VEC.FreeDofs()) * rhs
        gfX.vec.data += update
        #scene_gfset.Redraw()
   
    elif not linear:
        #print("solving nonlinear deformation equation with Newton")
        Newton(c, gfX, maxerr = 1E-8, freedofs=VEC.FreeDofs(), maxit = 1000, printing=True)

        #print("Deformation norm:", Norm(gfX.vec)
        #scene_gfset.Redraw()
        

checking the mesh quality#

We are estimating the meshquality by taking the elementwise (approximated) diameter

def check_mesh_quality():
    h = specialcf.mesh_size
    gf_h = GridFunction(L2(mesh, order=0))
    gf_h.Set(h)

    qual = gf_h.vec.FV().NumPy()
    qual /= maxH
    
    print("=================Mesh Quality Check=================")
    print("max quality:", np.max(qual), " min quality:", np.min(qual), " mean quality:", np.mean(qual))
    print("===================================================")
 
    #return Integrate(h, mesh, element_wise=True)
def optimize_shape(linear = True, LineSearch = True):
    gfset.Set((0,0))
    mesh.SetDeformation(gfset)
    scene.Redraw()
    a.Assemble()
    fstate.Assemble()
    SolveStateEquation()


    iter_max = 600
    Jold = Integrate(Cost(gfu), mesh)
    converged = False

    # input("Press enter to start optimization")
    print('Cost at initial design', Integrate (Cost(gfu), mesh))
    for k in range(iter_max):
        mesh.SetDeformation(gfset)
        scene.Redraw()

        print('cost at iteration', k, ': ', Jold)

        a.Assemble()
        fstate.Assemble()
        SolveStateEquation()

        fadjoint.Assemble()
        SolveAdjointEquation()
        if linear:
            b.Assemble()
        #plt.spy(b.mat.ToDense())
        dJOmega.Assemble()
        SolveDeformationEquation(linear=linear)

        #scale = 0.01 / Norm(gfX.vec)
        
        scale = 1E-4

        gfsetOld.vec[:] = gfset.vec[:]
        gfset.vec.data -= scale * gfX.vec

        Jnew = Integrate(Cost(gfu), mesh)
        if not LineSearch:
            if Jnew >= Jold:
                print("No descent found for current step size")
                converged = True

        #print("gradient norm", np.linalg.norm(dJOmega.vec.FV().NumPy()))
        if LineSearch:
            while Jnew > Jold and scale > 1e-12:
                #input('a')
                scale = min(1E-3, scale/2)
                print("scale = ", scale)
                #print("scale without normalization factor = ", scale* Norm(gfX.vec))
                if scale <= 1e-12:
                    converged = True
                    break

                gfset.vec.data = gfsetOld.vec - scale * gfX.vec         # scale times -X is descent direction
                mesh.SetDeformation(gfset)                              # applying deformation and solve new state

                a.Assemble()                                           # note that during line search we only have to solve state to check if update yields lower costs              
                fstate.Assemble()
                SolveStateEquation()
                Jnew = Integrate(Cost(gfu), mesh)
                print("New cost during line search: ", Jnew)
        
        #print("scaling parameter:", scale)
        #print("descet direction value", -1*np.dot(dJOmega.vec.FV().NumPy(), gfX.vec.FV().NumPy()))

        if converged==True:
            print("No more descent can be found")
            break
        Jold = Jnew

        Redraw(blocking=True)

    #print("Norm of gset", Integrate(gfset**2, mesh))

linear#

a1 = check_mesh_quality()
gfX = GridFunction(VEC)
b = pose_deformation_equation()

gfset = GridFunction(VEC)
gfset.Set((0,0))
gfsetOld = GridFunction(VEC)

scene = Draw(gfset,mesh,"gfset")
SetVisualization (deformation=True)



optimize_shape(linear = True)
mesh.SetDeformation(gfset)

a2 = check_mesh_quality()

# gfX2 = GridFunction(VEC)
# gfX2.vec.data =
# print("=====")
# print(np.dot(dJOmega.vec.FV().NumPy(), gfX.vec.FV().NumPy()))
# print(Integrate(Cost(gfu), mesh))
=================Mesh Quality Check=================
max quality: 1.452321790073649  min quality: 0.5768477774464785  mean quality: 0.8831324760843785
===================================================
Cost at initial design 13.96238810275433
cost at iteration 0 :  13.96238810275433
cost at iteration 1 :  13.735550739621843
cost at iteration 2 :  13.526642840331323
cost at iteration 3 :  13.322392417508466
cost at iteration 4 :  13.12265849375249
cost at iteration 5 :  12.927305689331154
cost at iteration 6 :  12.736203943744316
cost at iteration 7 :  12.549228253853594
cost at iteration 8 :  12.36625842744546
cost at iteration 9 :  12.187178851181772
cost at iteration 10 :  12.011878271971462
cost at iteration 11 :  11.840249590867716
cost at iteration 12 :  11.67218966866377
cost at iteration 13 :  11.507599142420847
cost at iteration 14 :  11.34638225221715
cost at iteration 15 :  11.188446677459389
cost at iteration 16 :  11.03370338214665
cost at iteration 17 :  10.88206646851788
cost at iteration 18 :  10.733453038557334
cost at iteration 19 :  10.587783062867315
cost at iteration 20 :  10.44497925645389
cost at iteration 21 :  10.304966961000465
cost at iteration 22 :  10.167674033236612
cost at iteration 23 :  10.03303073903283
cost at iteration 24 :  9.900969652880633
cost at iteration 25 :  9.771425562437086
cost at iteration 26 :  9.64433537783695
cost at iteration 27 :  9.519638045494702
cost at iteration 28 :  9.397274466134974
cost at iteration 29 :  9.277187416810861
cost at iteration 30 :  9.159321476681768
cost at iteration 31 :  9.043622956339174
cost at iteration 32 :  8.930039830481821
cost at iteration 33 :  8.818521673753528
cost at iteration 34 :  8.709019599570668
cost at iteration 35 :  8.601486201774575
cost at iteration 36 :  8.495875498956943
cost at iteration 37 :  8.392142881313857
cost at iteration 38 :  8.290245059893905
cost at iteration 39 :  8.190140018113214
cost at iteration 40 :  8.091786965419237
cost at iteration 41 :  7.995146292990532
cost at iteration 42 :  7.900179531368054
cost at iteration 43 :  7.806849309918279
cost at iteration 44 :  7.715119318036037
cost at iteration 45 :  7.624954267997505
cost at iteration 46 :  7.53631985938323
cost at iteration 47 :  7.449182744990578
cost at iteration 48 :  7.363510498164074
cost at iteration 49 :  7.279271581473583
cost at iteration 50 :  7.1964353166745845
cost at iteration 51 :  7.114971855889361
cost at iteration 52 :  7.034852153950404
cost at iteration 53 :  6.956047941851222
cost at iteration 54 :  6.878531701251429
cost at iteration 55 :  6.802276639988421
cost at iteration 56 :  6.727256668547503
cost at iteration 57 :  6.653446377446716
cost at iteration 58 :  6.580821015495286
cost at iteration 59 :  6.509356468885111
cost at iteration 60 :  6.4390292410776775
cost at iteration 61 :  6.369816433451991
cost at iteration 62 :  6.30169572667845
cost at iteration 63 :  6.2346453627871155
cost at iteration 64 :  6.168644127900511
cost at iteration 65 :  6.1036713356005805
cost at iteration 66 :  6.039706810904678
cost at iteration 67 :  5.97673087482172
cost at iteration 68 :  5.91472432946657
cost at iteration 69 :  5.853668443706879
cost at iteration 70 :  5.793544939320945
cost at iteration 71 :  5.734335977645499
cost at iteration 72 :  5.67602414669258
cost at iteration 73 :  5.618592448716386
cost at iteration 74 :  5.562024288211479
cost at iteration 75 :  5.506303460325568
cost at iteration 76 :  5.45141413966889
cost at iteration 77 :  5.3973408695056415
cost at iteration 78 :  5.344068551311047
cost at iteration 79 :  5.291582434680406
cost at iteration 80 :  5.239868107575666
cost at iteration 81 :  5.188911486897126
cost at iteration 82 :  5.138698809366406
cost at iteration 83 :  5.089216622710051
cost at iteration 84 :  5.040451777131078
cost at iteration 85 :  4.992391417058067
cost at iteration 86 :  4.9450229731612945
cost at iteration 87 :  4.898334154625466
cost at iteration 88 :  4.852312941669703
cost at iteration 89 :  4.806947578305664
cost at iteration 90 :  4.762226565324642
cost at iteration 91 :  4.718138653505423
cost at iteration 92 :  4.6746728370347395
cost at iteration 93 :  4.631818347132968
cost at iteration 94 :  4.589564645876579
cost at iteration 95 :  4.547901420211695
cost at iteration 96 :  4.506818576150781
cost at iteration 97 :  4.466306233146395
cost at iteration 98 :  4.426354718635686
cost at iteration 99 :  4.386954562749661
cost at iteration 100 :  4.348096493181397
cost at iteration 101 :  4.309771430207505
cost at iteration 102 :  4.271970481857837
cost at iteration 103 :  4.234684939228223
cost at iteration 104 :  4.197906271931019
cost at iteration 105 :  4.161626123679323
cost at iteration 106 :  4.125836307999768
cost at iteration 107 :  4.090528804069933
cost at iteration 108 :  4.055695752676047
cost at iteration 109 :  4.02132945228679
cost at iteration 110 :  3.98742235523977
cost at iteration 111 :  3.953967064036446
cost at iteration 112 :  3.9209563277422292
cost at iteration 113 :  3.8883830384882834
cost at iteration 114 :  3.856240228071364
cost at iteration 115 :  3.8245210646493395
cost at iteration 116 :  3.793218849528267
cost at iteration 117 :  3.7623270140387124
cost at iteration 118 :  3.7318391164983775
cost at iteration 119 :  3.7017488392583173
cost at iteration 120 :  3.672049985829624
cost at iteration 121 :  3.64273647808877
cost at iteration 122 :  3.6138023535586794
cost at iteration 123 :  3.5852417627631556
cost at iteration 124 :  3.557048966652588
cost at iteration 125 :  3.52921833409841
cost at iteration 126 :  3.501744339454925
cost at iteration 127 :  3.474621560184901
cost at iteration 128 :  3.44784467454924
cost at iteration 129 :  3.421408459356327
cost at iteration 130 :  3.3953077877712023
cost at iteration 131 :  3.369537627181565
cost at iteration 132 :  3.344093037119673
cost at iteration 133 :  3.3189691672374315
cost at iteration 134 :  3.294161255334572
cost at iteration 135 :  3.2696646254369646
cost at iteration 136 :  3.245474685924152
cost at iteration 137 :  3.221586927704939
cost at iteration 138 :  3.1979969224390046
cost at iteration 139 :  3.174700320803726
cost at iteration 140 :  3.151692850804505
cost at iteration 141 :  3.1289703161277425
cost at iteration 142 :  3.106528594534672
cost at iteration 143 :  3.0843636362954583
cost at iteration 144 :  3.062471462662174
cost at iteration 145 :  3.040848164379228
cost at iteration 146 :  3.019489900230858
cost at iteration 147 :  2.998392895623941
cost at iteration 148 :  2.9775534412056017
cost at iteration 149 :  2.9569678915144437
cost at iteration 150 :  2.9366326636644047
cost at iteration 151 :  2.916544236060737
cost at iteration 152 :  2.8966991471463954
cost at iteration 153 :  2.877093994178978
cost at iteration 154 :  2.8577254320364656
cost at iteration 155 :  2.838590172051993
cost at iteration 156 :  2.819684980875429
cost at iteration 157 :  2.8010066793627577
cost at iteration 158 :  2.7825521414908816
cost at iteration 159 :  2.7643182932984156
cost at iteration 160 :  2.746302111851069
cost at iteration 161 :  2.7285006242307546
cost at iteration 162 :  2.710910906549088
cost at iteration 163 :  2.693530082982694
cost at iteration 164 :  2.6763553248313885
cost at iteration 165 :  2.6593838495978726
cost at iteration 166 :  2.642612920088303
cost at iteration 167 :  2.626039843533671
cost at iteration 168 :  2.6096619707309516
cost at iteration 169 :  2.593476695203691
cost at iteration 170 :  2.577481452381752
cost at iteration 171 :  2.561673718799414
cost at iteration 172 :  2.546051011311209
cost at iteration 173 :  2.5306108863257855
cost at iteration 174 :  2.515350939056333
cost at iteration 175 :  2.5002688027880033
cost at iteration 176 :  2.4853621481613435
cost at iteration 177 :  2.470628682471653
cost at iteration 178 :  2.4560661489834916
cost at iteration 179 :  2.4416723262605333
cost at iteration 180 :  2.427445027509923
cost at iteration 181 :  2.413382099940702
cost at iteration 182 :  2.3994814241364444
cost at iteration 183 :  2.3857409134412375
cost at iteration 184 :  2.3721585133589196
cost at iteration 185 :  2.358732200965482
cost at iteration 186 :  2.3454599843335524
cost at iteration 187 :  2.3323399019696147
cost at iteration 188 :  2.3193700222628615
cost at iteration 189 :  2.3065484429460783
cost at iteration 190 :  2.293873290567446
cost at iteration 191 :  2.2813427199739977
cost at iteration 192 :  2.2689549138054415
cost at iteration 193 :  2.256708081998858
cost at iteration 194 :  2.244600461303798
cost at iteration 195 :  2.2326303148069733
cost at iteration 196 :  2.2207959314676313
cost at iteration 197 :  2.2090956256614747
cost at iteration 198 :  2.1975277367350303
cost at iteration 199 :  2.1860906285684303
cost at iteration 200 :  2.1747826891474222
cost at iteration 201 :  2.1636023301441427
cost at iteration 202 :  2.1525479865062165
cost at iteration 203 :  2.141618116054566
cost at iteration 204 :  2.130811199088807
cost at iteration 205 :  2.1201257380012137
cost at iteration 206 :  2.109560256897829
cost at iteration 207 :  2.099113301227706
cost at iteration 208 :  2.0887834374192007
cost at iteration 209 :  2.078569252523587
cost at iteration 210 :  2.0684693538658676
cost at iteration 211 :  2.0584823687022507
cost at iteration 212 :  2.0486069438847796
cost at iteration 213 :  2.038841745532014
cost at iteration 214 :  2.029185458706767
cost at iteration 215 :  2.0196367870996013
cost at iteration 216 :  2.010194452718936
cost at iteration 217 :  2.0008571955867405
cost at iteration 218 :  1.9916237734406657
cost at iteration 219 :  1.9824929614414906
cost at iteration 220 :  1.9734635518862256
cost at iteration 221 :  1.9645343539271067
cost at iteration 222 :  1.9557041932955332
cost at iteration 223 :  1.9469719120314786
cost at iteration 224 :  1.9383363682180776
cost at iteration 225 :  1.929796435721272
cost at iteration 226 :  1.9213510039341817
cost at iteration 227 :  1.9129989775265157
cost at iteration 228 :  1.9047392761987394
cost at iteration 229 :  1.8965708344405479
cost at iteration 230 :  1.8884926012941636
cost at iteration 231 :  1.8805035401220322
cost at iteration 232 :  1.8726026283785673
cost at iteration 233 :  1.8647888573865248
cost at iteration 234 :  1.8570612321171416
cost at iteration 235 :  1.8494187709746175
cost at iteration 236 :  1.841860505584373
cost at iteration 237 :  1.8343854805853856
cost at iteration 238 :  1.8269927534260753
cost at iteration 239 :  1.8196813941641619
cost at iteration 240 :  1.8124504852699752
cost at iteration 241 :  1.8052991214334695
cost at iteration 242 :  1.7982264093745732
cost at iteration 243 :  1.7912314676572485
cost at iteration 244 :  1.7843134265065905
cost at iteration 245 :  1.7774714276292956
cost at iteration 246 :  1.7707046240374826
cost at iteration 247 :  1.7640121798754804
cost at iteration 248 :  1.757393270249858
cost at iteration 249 :  1.7508470810623302
cost at iteration 250 :  1.7443728088456039
cost at iteration 251 :  1.737969660602488
cost at iteration 252 :  1.731636853647192
cost at iteration 253 :  1.725373615450051
cost at iteration 254 :  1.7191791834845536
cost at iteration 255 :  1.7130528050772313
cost at iteration 256 :  1.7069937372601294
cost at iteration 257 :  1.70100124662571
cost at iteration 258 :  1.6950746091846793
cost at iteration 259 :  1.6892131102256227
cost at iteration 260 :  1.683416044177652
cost at iteration 261 :  1.6776827144752418
cost at iteration 262 :  1.6720124334249806
cost at iteration 263 :  1.6664045220754165
cost at iteration 264 :  1.6608583100883452
cost at iteration 265 :  1.6553731356126518
cost at iteration 266 :  1.6499483451603318
cost at iteration 267 :  1.6445832934844677
cost at iteration 268 :  1.639277343459288
cost at iteration 269 :  1.6340298659624126
cost at iteration 270 :  1.6288402397588018
cost at iteration 271 :  1.6237078513869072
cost at iteration 272 :  1.6186320950464848
cost at iteration 273 :  1.6136123724886295
cost at iteration 274 :  1.6086480929071705
cost at iteration 275 :  1.6037386728322989
cost at iteration 276 :  1.5988835360257612
cost at iteration 277 :  1.5940821133777066
cost at iteration 278 :  1.5893338428053423
cost at iteration 279 :  1.5846381691534075
cost at iteration 280 :  1.5799945440958447
cost at iteration 281 :  1.5754024260395718
cost at iteration 282 :  1.5708612800294617
cost at iteration 283 :  1.5663705776550487
cost at iteration 284 :  1.56192979695875
cost at iteration 285 :  1.557538422345481
cost at iteration 286 :  1.5531959444936536
cost at iteration 287 :  1.548901860268023
cost at iteration 288 :  1.5446556726333833
cost at iteration 289 :  1.5404568905700111
cost at iteration 290 :  1.5363050289903868
cost at iteration 291 :  1.5321996086571132
cost at iteration 292 :  1.5281401561023675
cost at iteration 293 :  1.5241262035484744
cost at iteration 294 :  1.5201572888296888
cost at iteration 295 :  1.5162329553153822
cost at iteration 296 :  1.5123527518343562
cost at iteration 297 :  1.508516232600388
cost at iteration 298 :  1.5047229571387146
cost at iteration 299 :  1.5009724902141657
cost at iteration 300 :  1.4972644017598626
cost at iteration 301 :  1.493598266807315
cost at iteration 302 :  1.4899736654177438
cost at iteration 303 :  1.4863901826139896
cost at iteration 304 :  1.4828474083142
cost at iteration 305 :  1.4793449372656573
cost at iteration 306 :  1.475882368980472
cost at iteration 307 :  1.472459307671675
cost at iteration 308 :  1.4690753621905175
cost at iteration 309 :  1.4657301459648182
cost at iteration 310 :  1.4624232769379895
cost at iteration 311 :  1.459154377509365
cost at iteration 312 :  1.4559230744750364
cost at iteration 313 :  1.4527289989697263
cost at iteration 314 :  1.4495717864098823
cost at iteration 315 :  1.4464510764369465
cost at iteration 316 :  1.4433665128621622
cost at iteration 317 :  1.4403177436117856
cost at iteration 318 :  1.4373044206732672
cost at iteration 319 :  1.4343262000422319
cost at iteration 320 :  1.431382741670177
cost at iteration 321 :  1.428473709413187
cost at iteration 322 :  1.4255987709810212
cost at iteration 323 :  1.4227575978872453
cost at iteration 324 :  1.4199498654001614
cost at iteration 325 :  1.417175252494063
cost at iteration 326 :  1.4144334418017905
cost at iteration 327 :  1.4117241195673917
cost at iteration 328 :  1.409046975599941
cost at iteration 329 :  1.4064017032277107
cost at iteration 330 :  1.4037879992532847
cost at iteration 331 :  1.4012055639090877
cost at iteration 332 :  1.3986541008136941
cost at iteration 333 :  1.396133316928895
cost at iteration 334 :  1.3936429225169267
cost at iteration 335 :  1.3911826310988555
cost at iteration 336 :  1.3887521594132843
cost at iteration 337 :  1.3863512273758267
cost at iteration 338 :  1.3839795580386887
cost at iteration 339 :  1.3816368775516255
cost at iteration 340 :  1.3793229151226805
cost at iteration 341 :  1.3770374029799521
cost at iteration 342 :  1.3747800763339746
cost at iteration 343 :  1.372550673339937
cost at iteration 344 :  1.3703489350613844
cost at iteration 345 :  1.3681746054339314
cost at iteration 346 :  1.3660274312292202
cost at iteration 347 :  1.3639071620201206
cost at iteration 348 :  1.3618135501456547
cost at iteration 349 :  1.3597463506771514
cost at iteration 350 :  1.3577053213840993
cost at iteration 351 :  1.3556902227012104
cost at iteration 352 :  1.3537008176953111
cost at iteration 353 :  1.3517368720332115
cost at iteration 354 :  1.3497981539496846
cost at iteration 355 :  1.3478844342161382
cost at iteration 356 :  1.3459954861093983
cost at iteration 357 :  1.344131085381208
cost at iteration 358 :  1.342291010228171
cost at iteration 359 :  1.3404750412617914
cost at iteration 360 :  1.3386829614792632
cost at iteration 361 :  1.3369145562345277
cost at iteration 362 :  1.3351696132096162
cost at iteration 363 :  1.333447922386621
cost at iteration 364 :  1.3317492760198009
cost at iteration 365 :  1.3300734686084696
cost at iteration 366 :  1.3284202968695196
cost at iteration 367 :  1.3267895597111736
cost at iteration 368 :  1.325181058206404
cost at iteration 369 :  1.3235945955673651
cost at iteration 370 :  1.322029977119296
cost at iteration 371 :  1.3204870102757758
cost at iteration 372 :  1.3189655045134994
cost at iteration 373 :  1.3174652713478447
cost at iteration 374 :  1.3159861243086943
cost at iteration 375 :  1.314527878916308
cost at iteration 376 :  1.3130903526578683
cost at iteration 377 :  1.3116733649644474
cost at iteration 378 :  1.310276737187427
cost at iteration 379 :  1.3089002925764317
cost at iteration 380 :  1.3075438562566144
cost at iteration 381 :  1.306207255206912
cost at iteration 382 :  1.3048903182377662
cost at iteration 383 :  1.3035928759701794
cost at iteration 384 :  1.3023147608141397
cost at iteration 385 :  1.3010558069479445
cost at iteration 386 :  1.2998158502972734
cost at iteration 387 :  1.2985947285150612
cost at iteration 388 :  1.29739228096135
cost at iteration 389 :  1.2962083486832132
cost at iteration 390 :  1.2950427743954231
cost at iteration 391 :  1.293895402461041
cost at iteration 392 :  1.2927660788722422
cost at iteration 393 :  1.2916546512315608
cost at iteration 394 :  1.2905609687333333
cost at iteration 395 :  1.2894848821451101
cost at iteration 396 :  1.2884262437899847
cost at iteration 397 :  1.2873849075281925
cost at iteration 398 :  1.2863607287398187
cost at iteration 399 :  1.2853535643072802
cost at iteration 400 :  1.2843632725979808
cost at iteration 401 :  1.2833897134475774
cost at iteration 402 :  1.2824327481429731
cost at iteration 403 :  1.2814922394058883
cost at iteration 404 :  1.2805680513765503
cost at iteration 405 :  1.2796600495974786
cost at iteration 406 :  1.2787681009975265
cost at iteration 407 :  1.2778920738762698
cost at iteration 408 :  1.2770318378884116
cost at iteration 409 :  1.2761872640284655
cost at iteration 410 :  1.2753582246155069
cost at iteration 411 :  1.2745445932784343
cost at iteration 412 :  1.273746244941086
cost at iteration 413 :  1.2729630558076461
cost at iteration 414 :  1.272194903348229
cost at iteration 415 :  1.271441666284778
cost at iteration 416 :  1.2707032245769816
cost at iteration 417 :  1.2699794594083245
cost at iteration 418 :  1.2692702531725586
cost at iteration 419 :  1.268575489460036
cost at iteration 420 :  1.2678950530445068
cost at iteration 421 :  1.2672288298697183
cost at iteration 422 :  1.266576707036594
cost at iteration 423 :  1.2659385727904116
cost at iteration 424 :  1.2653143165078875
cost at iteration 425 :  1.264703828684697
cost at iteration 426 :  1.2641070009232351
cost at iteration 427 :  1.2635237259201735
cost at iteration 428 :  1.2629538974546575
cost at iteration 429 :  1.2623974103758737
cost at iteration 430 :  1.2618541605918403
cost at iteration 431 :  1.2613240450573697
cost at iteration 432 :  1.2608069617625788
cost at iteration 433 :  1.260302809721812
cost at iteration 434 :  1.259811488962035
cost at iteration 435 :  1.259332900511973
cost at iteration 436 :  1.258866946391148
cost at iteration 437 :  1.2584135295989214
cost at iteration 438 :  1.2579725541039437
cost at iteration 439 :  1.2575439248334177
cost at iteration 440 :  1.2571275476629469
cost at iteration 441 :  1.2567233294058464
cost at iteration 442 :  1.256331177803315
cost at iteration 443 :  1.255951001514014
cost at iteration 444 :  1.2555827101045347
cost at iteration 445 :  1.2552262140391546
cost at iteration 446 :  1.2548814246703741
cost at iteration 447 :  1.254548254229384
cost at iteration 448 :  1.254226615816268
cost at iteration 449 :  1.2539164233909819
cost at iteration 450 :  1.2536175917640013
cost at iteration 451 :  1.2533300365868727
cost at iteration 452 :  1.2530536743438068
cost at iteration 453 :  1.25278842234207
cost at iteration 454 :  1.2525341987036143
cost at iteration 455 :  1.252290922356204
cost at iteration 456 :  1.2520585130247126
cost at iteration 457 :  1.251836891222852
cost at iteration 458 :  1.2516259782445036
cost at iteration 459 :  1.2514256961555208
cost at iteration 460 :  1.251235967785526
cost at iteration 461 :  1.2510567167196875
cost at iteration 462 :  1.2508878672909112
cost at iteration 463 :  1.2507293445715988
cost at iteration 464 :  1.2505810743660208
cost at iteration 465 :  1.2504429832025268
cost at iteration 466 :  1.2503149983259212
cost at iteration 467 :  1.2501970476897117
cost at iteration 468 :  1.250089059948902
cost at iteration 469 :  1.2499909644523188
cost at iteration 470 :  1.249902691235647
cost at iteration 471 :  1.2498241710136928
cost at iteration 472 :  1.2497553351738397
cost at iteration 473 :  1.249696115768505
cost at iteration 474 :  1.249646445508458
cost at iteration 475 :  1.2496062577556994
cost at iteration 476 :  1.2495754865168993
cost at iteration 477 :  1.2495540664364273
cost at iteration 478 :  1.2495419327897916
cost at iteration 479 :  1.2495390214770865
scale =  5e-05
New cost during line search:  1.242616032021647
cost at iteration 480 :  1.242616032021647
scale =  5e-05
New cost during line search:  1.2426156656099823
cost at iteration 481 :  1.2426156656099823
scale =  5e-05
New cost during line search:  1.242617568717261
scale =  2.5e-05
New cost during line search:  1.2426165396036053
scale =  1.25e-05
New cost during line search:  1.2426160832181463
scale =  6.25e-06
New cost during line search:  1.2426158695670388
scale =  3.125e-06
New cost during line search:  1.24261576637678
scale =  1.5625e-06
New cost during line search:  1.2426157156904207
scale =  7.8125e-07
New cost during line search:  1.2426156905745076
scale =  3.90625e-07
New cost during line search:  1.242615678073288
scale =  1.953125e-07
New cost during line search:  1.2426156718368953
scale =  9.765625e-08
New cost during line search:  1.2426156687222547
scale =  4.8828125e-08
New cost during line search:  1.2426156671658437
scale =  2.44140625e-08
New cost during line search:  1.242615666387826
scale =  1.220703125e-08
New cost during line search:  1.2426156659988907
scale =  6.103515625e-09
New cost during line search:  1.2426156658044298
scale =  3.0517578125e-09
New cost during line search:  1.2426156657071838
scale =  1.52587890625e-09
New cost during line search:  1.2426156656586058
scale =  7.62939453125e-10
New cost during line search:  1.242615665634289
scale =  3.814697265625e-10
New cost during line search:  1.242615665622116
scale =  1.9073486328125e-10
New cost during line search:  1.2426156656160388
scale =  9.5367431640625e-11
New cost during line search:  1.2426156656130096
scale =  4.76837158203125e-11
New cost during line search:  1.2426156656114715
scale =  2.384185791015625e-11
New cost during line search:  1.2426156656107425
scale =  1.1920928955078126e-11
New cost during line search:  1.2426156656103344
scale =  5.960464477539063e-12
New cost during line search:  1.2426156656101475
scale =  2.9802322387695314e-12
New cost during line search:  1.2426156656100888
scale =  1.4901161193847657e-12
New cost during line search:  1.2426156656099985
scale =  7.450580596923828e-13
No more descent can be found
=================Mesh Quality Check=================
max quality: 1.4420554859725245  min quality: 0.5166402075400832  mean quality: 0.8593656411107816
===================================================
check_mesh_quality()
gfX = GridFunction(VEC)

b = pose_deformation_equation(gamma=1)

gfset = GridFunction(VEC)
gfset.Set((0,0))
gfsetOld = GridFunction(VEC)

scene = Draw(gfset,mesh,"gfset")
SetVisualization (deformation=True)

optimize_shape(linear = True)

mesh.SetDeformation(gfset)
check_mesh_quality()
=================Mesh Quality Check=================
max quality: 1.4420554859725245  min quality: 0.5166402075400832  mean quality: 0.8593656411107816
===================================================
Using gamma = 1
Cost at initial design 13.962388102754325
cost at iteration 0 :  13.962388102754325
cost at iteration 1 :  13.825772613698412
cost at iteration 2 :  13.708401589313848
cost at iteration 3 :  13.592734781168076
cost at iteration 4 :  13.478731015689526
cost at iteration 5 :  13.366350749245594
cost at iteration 6 :  13.255555977256009
cost at iteration 7 :  13.146310149336143
cost at iteration 8 :  13.038578090036056
cost at iteration 9 :  12.932325924773831
cost at iteration 10 :  12.827521010591616
cost at iteration 11 :  12.724131871390405
cost at iteration 12 :  12.62212813732685
cost at iteration 13 :  12.521480488076158
cost at iteration 14 :  12.422160599690764
cost at iteration 15 :  12.324141094801222
cost at iteration 16 :  12.227395495926604
cost at iteration 17 :  12.131898181677869
cost at iteration 18 :  12.03762434565384
cost at iteration 19 :  11.944549957844472
cost at iteration 20 :  11.85265172836919
cost at iteration 21 :  11.761907073390157
cost at iteration 22 :  11.672294083053641
cost at iteration 23 :  11.583791491321508
cost at iteration 24 :  11.496378647565159
cost at iteration 25 :  11.410035489805589
cost at iteration 26 :  11.324742519486739
cost at iteration 27 :  11.240480777683652
cost at iteration 28 :  11.15723182264898
cost at iteration 29 :  11.074977708610167
cost at iteration 30 :  10.993700965735849
cost at iteration 31 :  10.913384581196372
cost at iteration 32 :  10.834011981246448
cost at iteration 33 :  10.755567014265564
cost at iteration 34 :  10.678033934694788
cost at iteration 35 :  10.60139738781274
cost at iteration 36 :  10.525642395299371
cost at iteration 37 :  10.4507543415359
cost at iteration 38 :  10.376718960597987
cost at iteration 39 :  10.30352232389744
cost at iteration 40 :  10.231150828433448
cost at iteration 41 :  10.159591185616444
cost at iteration 42 :  10.088830410629775
cost at iteration 43 :  10.018855812297367
cost at iteration 44 :  9.94965498342675
cost at iteration 45 :  9.881215791600304
cost at iteration 46 :  9.813526370387274
cost at iteration 47 :  9.746575110953774
cost at iteration 48 :  9.680350654046608
cost at iteration 49 :  9.614841882330133
cost at iteration 50 :  9.550037913055997
cost at iteration 51 :  9.485928091047798
cost at iteration 52 :  9.42250198198223
cost at iteration 53 :  9.359749365951204
cost at iteration 54 :  9.29766023128913
cost at iteration 55 :  9.236224768651516
cost at iteration 56 :  9.175433365331248
cost at iteration 57 :  9.115276599799587
cost at iteration 58 :  9.055745236461473
cost at iteration 59 :  8.996830220611923
cost at iteration 60 :  8.938522673585503
cost at iteration 61 :  8.880813888087435
cost at iteration 62 :  8.82369532369782
cost at iteration 63 :  8.76715860254074
cost at iteration 64 :  8.711195505109323
cost at iteration 65 :  8.655797966239959
cost at iteration 66 :  8.600958071228272
cost at iteration 67 :  8.546668052079939
cost at iteration 68 :  8.492920283890067
cost at iteration 69 :  8.439707281345619
cost at iteration 70 :  8.387021695344725
cost at iteration 71 :  8.33485630972766
cost at iteration 72 :  8.283204038114569
cost at iteration 73 :  8.232057920845817
cost at iteration 74 :  8.181411122019073
cost at iteration 75 :  8.131256926620479
cost at iteration 76 :  8.081588737744742
cost at iteration 77 :  8.032400073901128
cost at iteration 78 :  7.983684566401613
cost at iteration 79 :  7.9354359568272175
cost at iteration 80 :  7.887648094570333
cost at iteration 81 :  7.840314934449432
cost at iteration 82 :  7.7934305343928365
cost at iteration 83 :  7.746989053190266
cost at iteration 84 :  7.700984748307666
cost at iteration 85 :  7.655411973764373
cost at iteration 86 :  7.610265178069866
cost at iteration 87 :  7.565538902217125
cost at iteration 88 :  7.521227777731821
cost at iteration 89 :  7.4773265247739875
cost at iteration 90 :  7.433829950291335
cost at iteration 91 :  7.390732946221673
cost at iteration 92 :  7.348030487743153
cost at iteration 93 :  7.30571763157005
cost at iteration 94 :  7.263789514293292
cost at iteration 95 :  7.222241350763824
cost at iteration 96 :  7.181068432517166
cost at iteration 97 :  7.140266126237953
cost at iteration 98 :  7.099829872263646
cost at iteration 99 :  7.0597551831250165
cost at iteration 100 :  7.020037642123312
cost at iteration 101 :  6.980672901942279
cost at iteration 102 :  6.9416566832942195
cost at iteration 103 :  6.902984773598968
cost at iteration 104 :  6.864653025694407
cost at iteration 105 :  6.8266573565784245
cost at iteration 106 :  6.788993746180299
cost at iteration 107 :  6.751658236161375
cost at iteration 108 :  6.7146469287437895
cost at iteration 109 :  6.677955985566658
cost at iteration 110 :  6.641581626568691
cost at iteration 111 :  6.605520128896418
cost at iteration 112 :  6.569767825837985
cost at iteration 113 :  6.534321105780406
cost at iteration 114 :  6.499176411191092
cost at iteration 115 :  6.4643302376220655
cost at iteration 116 :  6.42977913273647
cost at iteration 117 :  6.395519695356891
cost at iteration 118 :  6.361548574534466
cost at iteration 119 :  6.327862468638919
cost at iteration 120 :  6.294458124468221
cost at iteration 121 :  6.261332336377701
cost at iteration 122 :  6.228481945428376
cost at iteration 123 :  6.195903838553306
cost at iteration 124 :  6.163594947741946
cost at iteration 125 :  6.131552249242174
cost at iteration 126 :  6.099772762778906
cost at iteration 127 :  6.068253550789795
cost at iteration 128 :  6.03699171767664
cost at iteration 129 :  6.005984409072565
cost at iteration 130 :  5.975228811124798
cost at iteration 131 :  5.944722149792106
cost at iteration 132 :  5.914461690156875
cost at iteration 133 :  5.88444473575141
cost at iteration 134 :  5.854668627898181
cost at iteration 135 :  5.825130745063481
cost at iteration 136 :  5.79582850222414
cost at iteration 137 :  5.766759350247391
cost at iteration 138 :  5.737920775283032
cost at iteration 139 :  5.70931029816814
cost at iteration 140 :  5.680925473843277
cost at iteration 141 :  5.652763890780821
cost at iteration 142 :  5.624823170424442
cost at iteration 143 :  5.59710096663976
cost at iteration 144 :  5.569594965175941
cost at iteration 145 :  5.542302883137909
cost at iteration 146 :  5.515222468468557
cost at iteration 147 :  5.488351499441873
cost at iteration 148 :  5.461687784165084
cost at iteration 149 :  5.435229160090971
cost at iteration 150 :  5.408973493539715
cost at iteration 151 :  5.3829186792295065
cost at iteration 152 :  5.357062639816603
cost at iteration 153 :  5.331403325444119
cost at iteration 154 :  5.305938713299228
cost at iteration 155 :  5.280666807179087
cost at iteration 156 :  5.255585637064755
cost at iteration 157 :  5.230693258703373
cost at iteration 158 :  5.205987753198138
cost at iteration 159 :  5.181467226605926
cost at iteration 160 :  5.157129809542473
cost at iteration 161 :  5.132973656795099
cost at iteration 162 :  5.1089969469424465
cost at iteration 163 :  5.085197881981353
cost at iteration 164 :  5.0615746869606255
cost at iteration 165 :  5.038125609621544
cost at iteration 166 :  5.014848920045193
cost at iteration 167 :  4.991742910305761
cost at iteration 168 :  4.968805894130814
cost at iteration 169 :  4.946036206567214
cost at iteration 170 :  4.923432203653479
cost at iteration 171 :  4.90099226209791
cost at iteration 172 :  4.878714778962587
cost at iteration 173 :  4.8565981713529744
cost at iteration 174 :  4.834640876113276
cost at iteration 175 :  4.812841349527101
cost at iteration 176 :  4.791198067023474
cost at iteration 177 :  4.769709522888272
cost at iteration 178 :  4.748374229980345
cost at iteration 179 :  4.727190719453457
cost at iteration 180 :  4.706157540482098
cost at iteration 181 :  4.685273259993144
cost at iteration 182 :  4.664536462401425
cost at iteration 183 :  4.643945749350507
cost at iteration 184 :  4.6234997394575155
cost at iteration 185 :  4.603197068062952
cost at iteration 186 :  4.583036386984256
cost at iteration 187 :  4.563016364274103
cost at iteration 188 :  4.543135683982712
cost at iteration 189 :  4.523393045924095
cost at iteration 190 :  4.50378716544673
cost at iteration 191 :  4.484316773207632
cost at iteration 192 :  4.464980614950813
cost at iteration 193 :  4.445777451289177
cost at iteration 194 :  4.426706057490039
cost at iteration 195 :  4.407765223264819
cost at iteration 196 :  4.388953752561654
cost at iteration 197 :  4.370270463361901
cost at iteration 198 :  4.351714187479962
cost at iteration 199 :  4.333283770366416
cost at iteration 200 :  4.3149780709145515
cost at iteration 201 :  4.296795961270009
cost at iteration 202 :  4.27873632664377
cost at iteration 203 :  4.260798065127959
cost at iteration 204 :  4.242980087514956
cost at iteration 205 :  4.225281317119616
cost at iteration 206 :  4.207700689603904
cost at iteration 207 :  4.190237152804947
cost at iteration 208 :  4.1728896665655855
cost at iteration 209 :  4.155657202567927
cost at iteration 210 :  4.138538744169364
cost at iteration 211 :  4.1215332862414815
cost at iteration 212 :  4.10463983501145
cost at iteration 213 :  4.087857407906072
cost at iteration 214 :  4.071185033398243
cost at iteration 215 :  4.054621750856027
cost at iteration 216 :  4.0381666103940494
cost at iteration 217 :  4.021818672727305
cost at iteration 218 :  4.005577009027234
cost at iteration 219 :  3.989440700780418
cost at iteration 220 :  3.973408839648983
cost at iteration 221 :  3.957480527333594
cost at iteration 222 :  3.941654875438743
cost at iteration 223 :  3.925931005339714
cost at iteration 224 :  3.9103080480520895
cost at iteration 225 :  3.894785144102904
cost at iteration 226 :  3.8793614434042225
cost at iteration 227 :  3.864036105128496
cost at iteration 228 :  3.848808297585798
cost at iteration 229 :  3.8336771981030195
cost at iteration 230 :  3.8186419929052002
cost at iteration 231 :  3.8037018769982414
cost at iteration 232 :  3.788856054053679
cost at iteration 233 :  3.7741037362954715
cost at iteration 234 :  3.759444144388007
cost at iteration 235 :  3.7448765073262873
cost at iteration 236 :  3.7304000623275106
cost at iteration 237 :  3.7160140547244733
cost at iteration 238 :  3.7017177378606054
cost at iteration 239 :  3.6875103729864156
cost at iteration 240 :  3.6733912291577493
cost at iteration 241 :  3.659359583135402
cost at iteration 242 :  3.6454147192863853
cost at iteration 243 :  3.631555929486598
cost at iteration 244 :  3.617782513024946
cost at iteration 245 :  3.6040937765090435
cost at iteration 246 :  3.590489033772125
cost at iteration 247 :  3.5769676057815403
cost at iteration 248 :  3.5635288205483513
cost at iteration 249 :  3.5501720130387353
cost at iteration 250 :  3.5368965250861537
cost at iteration 251 :  3.5237017053051995
cost at iteration 252 :  3.510586909006592
cost at iteration 253 :  3.4975514981134235
cost at iteration 254 :  3.4845948410786276
cost at iteration 255 :  3.4717163128037187
cost at iteration 256 :  3.458915294558712
cost at iteration 257 :  3.4461911739031206
cost at iteration 258 :  3.4335433446081187
cost at iteration 259 :  3.420971206580122
cost at iteration 260 :  3.4084741657849307
cost at iteration 261 :  3.396051634173521
cost at iteration 262 :  3.383703029608659
cost at iteration 263 :  3.3714277757923736
cost at iteration 264 :  3.359225302194955
cost at iteration 265 :  3.347095043984623
cost at iteration 266 :  3.3350364419580947
cost at iteration 267 :  3.32304894247264
cost at iteration 268 :  3.3111319973786593
cost at iteration 269 :  3.299285063953235
cost at iteration 270 :  3.287507604834982
cost at iteration 271 :  3.2757990879595726
cost at iteration 272 :  3.2641589864959597
cost at iteration 273 :  3.2525867787840115
cost at iteration 274 :  3.2410819482725657
cost at iteration 275 :  3.2296439834587516
cost at iteration 276 :  3.218272377827663
cost at iteration 277 :  3.206966629793401
cost at iteration 278 :  3.1957262426406157
cost at iteration 279 :  3.1845507244668863
cost at iteration 280 :  3.1734395881261572
cost at iteration 281 :  3.162392351172511
cost at iteration 282 :  3.1514085358052246
cost at iteration 283 :  3.1404876688140906
cost at iteration 284 :  3.129629281526015
cost at iteration 285 :  3.1188329097518332
cost at iteration 286 :  3.108098093734223
cost at iteration 287 :  3.0974243780961825
cost at iteration 288 :  3.086811311790229
cost at iteration 289 :  3.0762584480482973
cost at iteration 290 :  3.065765344332298
cost at iteration 291 :  3.0553315622854633
cost at iteration 292 :  3.044956667684129
cost at iteration 293 :  3.034640230390382
cost at iteration 294 :  3.024381824305271
cost at iteration 295 :  3.014181027322641
cost at iteration 296 :  3.004037421283487
cost at iteration 297 :  2.993950591931224
cost at iteration 298 :  2.9839201288671626
cost at iteration 299 :  2.9739456255068797
cost at iteration 300 :  2.9640266790370684
cost at iteration 301 :  2.95416289037284
cost at iteration 302 :  2.9443538641159943
cost at iteration 303 :  2.93459920851317
cost at iteration 304 :  2.9248985354152612
cost at iteration 305 :  2.9152514602369415
cost at iteration 306 :  2.905657601916772
cost at iteration 307 :  2.8961165828780073
cost at iteration 308 :  2.8866280289897066
cost at iteration 309 :  2.877191569528519
cost at iteration 310 :  2.867806837140803
cost at iteration 311 :  2.8584734678055432
cost at iteration 312 :  2.849191100797221
cost at iteration 313 :  2.839959378649823
cost at iteration 314 :  2.830777947120737
cost at iteration 315 :  2.8216464551555402
cost at iteration 316 :  2.81256455485293
cost at iteration 317 :  2.803531901430348
cost at iteration 318 :  2.794548153189873
cost at iteration 319 :  2.785612971484631
cost at iteration 320 :  2.7767260206856648
cost at iteration 321 :  2.7678869681491616
cost at iteration 322 :  2.7590954841840993
cost at iteration 323 :  2.750351242020376
cost at iteration 324 :  2.7416539177772665
cost at iteration 325 :  2.733003190432423
cost at iteration 326 :  2.7243987417909237
cost at iteration 327 :  2.715840256455236
cost at iteration 328 :  2.7073274217950845
cost at iteration 329 :  2.6988599279179004
cost at iteration 330 :  2.690437467639729
cost at iteration 331 :  2.6820597364562824
cost at iteration 332 :  2.673726432514631
cost at iteration 333 :  2.6654372565848696
cost at iteration 334 :  2.6571919120325873
cost at iteration 335 :  2.6489901047912983
cost at iteration 336 :  2.6408315433355103
cost at iteration 337 :  2.6327159386538064
cost at iteration 338 :  2.624643004222619
cost at iteration 339 :  2.6166124559800505
cost at iteration 340 :  2.608624012300176
cost at iteration 341 :  2.6006773939674965
cost at iteration 342 :  2.5927723241520804
cost at iteration 343 :  2.584908528384302
cost at iteration 344 :  2.5770857345308076
cost at iteration 345 :  2.56930367276997
cost at iteration 346 :  2.5615620755681605
cost at iteration 347 :  2.5538606776560253
cost at iteration 348 :  2.546199216005103
cost at iteration 349 :  2.5385774298049624
cost at iteration 350 :  2.5309950604402665
cost at iteration 351 :  2.5234518514682605
cost at iteration 352 :  2.515947548596616
cost at iteration 353 :  2.5084818996615112
cost at iteration 354 :  2.501054654605828
cost at iteration 355 :  2.4936655654577473
cost at iteration 356 :  2.4863143863096298
cost at iteration 357 :  2.4790008732970406
cost at iteration 358 :  2.471724784578061
cost at iteration 359 :  2.4644858803129277
cost at iteration 360 :  2.4572839226437084
cost at iteration 361 :  2.450118675674617
cost at iteration 362 :  2.4429899054519395
cost at iteration 363 :  2.435897379944866
cost at iteration 364 :  2.4288408690261813
cost at iteration 365 :  2.421820144453088
cost at iteration 366 :  2.414834979848614
cost at iteration 367 :  2.407885150682844
cost at iteration 368 :  2.4009704342548415
cost at iteration 369 :  2.3940906096741497
cost at iteration 370 :  2.387245457843142
cost at iteration 371 :  2.3804347614392123
cost at iteration 372 :  2.3736583048972406
cost at iteration 373 :  2.3669158743922702
cost at iteration 374 :  2.360207257822471
cost at iteration 375 :  2.353532244792174
cost at iteration 376 :  2.346890626595135
cost at iteration 377 :  2.3402821961981233
cost at iteration 378 :  2.3337067482244263
cost at iteration 379 :  2.3271640789378374
cost at iteration 380 :  2.3206539862266475
cost at iteration 381 :  2.3141762695878505
cost at iteration 382 :  2.307730730111632
cost at iteration 383 :  2.301317170465829
cost at iteration 384 :  2.294935394880756
cost at iteration 385 :  2.2885852091342054
cost at iteration 386 :  2.28226642053642
cost at iteration 387 :  2.2759788379154244
cost at iteration 388 :  2.269722271602483
cost at iteration 389 :  2.2634965334176997
cost at iteration 390 :  2.2573014366557596
cost at iteration 391 :  2.2511367960719033
cost at iteration 392 :  2.245002427867939
cost at iteration 393 :  2.2388981496785156
cost at iteration 394 :  2.2328237805576
cost at iteration 395 :  2.226779140964867
cost at iteration 396 :  2.22076405275258
cost at iteration 397 :  2.2147783391522573
cost at iteration 398 :  2.2088218247617903
cost at iteration 399 :  2.2028943355326214
cost at iteration 400 :  2.1969956987567536
cost at iteration 401 :  2.1911257430545645
cost at iteration 402 :  2.185284298362089
cost at iteration 403 :  2.179471195918811
cost at iteration 404 :  2.1736862682554667
cost at iteration 405 :  2.167929349182026
cost at iteration 406 :  2.162200273775836
cost at iteration 407 :  2.1564988783697827
cost at iteration 408 :  2.150825000540728
cost at iteration 409 :  2.1451784790978556
cost at iteration 410 :  2.1395591540714687
cost at iteration 411 :  2.133966866701642
cost at iteration 412 :  2.1284014594269918
cost at iteration 413 :  2.122862775873796
cost at iteration 414 :  2.1173506608450197
cost at iteration 415 :  2.1118649603095836
cost at iteration 416 :  2.106405521391615
cost at iteration 417 :  2.1009721923599995
cost at iteration 418 :  2.0955648226178654
cost at iteration 419 :  2.0901832626922348
cost at iteration 420 :  2.084827364223976
cost at iteration 421 :  2.07949697995747
cost at iteration 422 :  2.074191963730797
cost at iteration 423 :  2.0689121704657145
cost at iteration 424 :  2.063657456158043
cost at iteration 425 :  2.0584276778677353
cost at iteration 426 :  2.053222693709529
cost at iteration 427 :  2.0480423628433977
cost at iteration 428 :  2.0428865454651546
cost at iteration 429 :  2.037755102797126
cost at iteration 430 :  2.032647897079085
cost at iteration 431 :  2.0275647915591453
cost at iteration 432 :  2.0225056504846646
cost at iteration 433 :  2.017470339093534
cost at iteration 434 :  2.0124587236052056
cost at iteration 435 :  2.0074706712121477
cost at iteration 436 :  2.0025060500710534
cost at iteration 437 :  1.9975647292945102
cost at iteration 438 :  1.9926465789423775
cost at iteration 439 :  1.9877514700136314
cost at iteration 440 :  1.982879274437885
cost at iteration 441 :  1.9780298650674404
cost at iteration 442 :  1.9732031156690253
cost at iteration 443 :  1.9683989009158784
cost at iteration 444 :  1.9636170963797466
cost at iteration 445 :  1.9588575785232047
cost at iteration 446 :  1.9541202246916343
cost at iteration 447 :  1.9494049131058178
cost at iteration 448 :  1.9447115228541376
cost at iteration 449 :  1.9400399338851813
cost at iteration 450 :  1.9353900270002058
cost at iteration 451 :  1.930761683845777
cost at iteration 452 :  1.9261547869065354
cost at iteration 453 :  1.921569219497869
cost at iteration 454 :  1.9170048657589447
cost at iteration 455 :  1.9124616106453542
cost at iteration 456 :  1.9079393399224565
cost at iteration 457 :  1.9034379401581238
cost at iteration 458 :  1.898957298716071
cost at iteration 459 :  1.8944973037489519
cost at iteration 460 :  1.8900578441918245
cost at iteration 461 :  1.88563880975517
cost at iteration 462 :  1.8812400909186653
cost at iteration 463 :  1.8768615789243626
cost at iteration 464 :  1.8725031657704054
cost at iteration 465 :  1.8681647442045735
cost at iteration 466 :  1.8638462077179465
cost at iteration 467 :  1.8595474505386596
cost at iteration 468 :  1.8552683676256985
cost at iteration 469 :  1.8510088546627164
cost at iteration 470 :  1.846768808051967
cost at iteration 471 :  1.8425481249083493
cost at iteration 472 :  1.8383467030533396
cost at iteration 473 :  1.8341644410092341
cost at iteration 474 :  1.8300012379930357
cost at iteration 475 :  1.8258569939110172
cost at iteration 476 :  1.8217316093526965
cost at iteration 477 :  1.8176249855853683
cost at iteration 478 :  1.81353702454827
cost at iteration 479 :  1.8094676288472926
cost at iteration 480 :  1.8054167017491978
cost at iteration 481 :  1.8013841471763428
cost at iteration 482 :  1.7973698697012255
cost at iteration 483 :  1.7933737745411453
cost at iteration 484 :  1.7893957675528274
cost at iteration 485 :  1.785435755227312
cost at iteration 486 :  1.781493644684605
cost at iteration 487 :  1.7775693436687017
cost at iteration 488 :  1.7736627605423319
cost at iteration 489 :  1.7697738042820614
cost at iteration 490 :  1.76590238447314
cost at iteration 491 :  1.7620484113046821
cost at iteration 492 :  1.7582117955646677
cost at iteration 493 :  1.7543924486351927
cost at iteration 494 :  1.7505902824875663
cost at iteration 495 :  1.7468052096775428
cost at iteration 496 :  1.7430371433407568
cost at iteration 497 :  1.7392859971877888
cost at iteration 498 :  1.735551685499818
cost at iteration 499 :  1.7318341231238346
cost at iteration 500 :  1.7281332254681958
cost at iteration 501 :  1.7244489084981307
cost at iteration 502 :  1.7207810887312072
cost at iteration 503 :  1.7171296832329783
cost at iteration 504 :  1.7134946096127064
cost at iteration 505 :  1.7098757860188099
cost at iteration 506 :  1.7062731311347472
cost at iteration 507 :  1.7026865641747577
cost at iteration 508 :  1.6991160048795564
cost at iteration 509 :  1.6955613735122526
cost at iteration 510 :  1.6920225908542317
cost at iteration 511 :  1.688499578200949
cost at iteration 512 :  1.6849922573579899
cost at iteration 513 :  1.6815005506370084
cost at iteration 514 :  1.678024380851817
cost at iteration 515 :  1.6745636713143124
cost at iteration 516 :  1.6711183458306498
cost at iteration 517 :  1.667688328697402
cost at iteration 518 :  1.6642735446976848
cost at iteration 519 :  1.6608739190973945
cost at iteration 520 :  1.6574893776413469
cost at iteration 521 :  1.6541198465497196
cost at iteration 522 :  1.6507652525142187
cost at iteration 523 :  1.6474255226944374
cost at iteration 524 :  1.6441005847143233
cost at iteration 525 :  1.6407903666584605
cost at iteration 526 :  1.6374947970686191
cost at iteration 527 :  1.6342138049401824
cost at iteration 528 :  1.6309473197185778
cost at iteration 529 :  1.6276952712959678
cost at iteration 530 :  1.6244575900077578
cost at iteration 531 :  1.6212342066290346
cost at iteration 532 :  1.6180250523715307
cost at iteration 533 :  1.6148300588799223
cost at iteration 534 :  1.61164915822883
cost at iteration 535 :  1.6084822829193206
cost at iteration 536 :  1.6053293658757424
cost at iteration 537 :  1.6021903404425513
cost at iteration 538 :  1.5990651403810356
cost at iteration 539 :  1.595953699866223
cost at iteration 540 :  1.5928559534836677
cost at iteration 541 :  1.5897718362264412
cost at iteration 542 :  1.5867012834919412
cost at iteration 543 :  1.5836442310790193
cost at iteration 544 :  1.5806006151847636
cost at iteration 545 :  1.5775703724015537
cost at iteration 546 :  1.5745534397141807
cost at iteration 547 :  1.5715497544968733
cost at iteration 548 :  1.5685592545102658
cost at iteration 549 :  1.565581877898658
cost at iteration 550 :  1.562617563187081
cost at iteration 551 :  1.5596662492784263
cost at iteration 552 :  1.5567278754507556
cost at iteration 553 :  1.5538023813542918
cost at iteration 554 :  1.5508897070089045
cost at iteration 555 :  1.5479897928011683
cost at iteration 556 :  1.5451025794818032
cost at iteration 557 :  1.5422280081628288
cost at iteration 558 :  1.5393660203150048
cost at iteration 559 :  1.5365165577651203
cost at iteration 560 :  1.5336795626935251
cost at iteration 561 :  1.530854977631256
cost at iteration 562 :  1.528042745457633
cost at iteration 563 :  1.5252428093977621
cost at iteration 564 :  1.522455113019884
cost at iteration 565 :  1.519679600232826
cost at iteration 566 :  1.516916215283718
cost at iteration 567 :  1.514164902755255
cost at iteration 568 :  1.5114256075635313
cost at iteration 569 :  1.5086982749553746
cost at iteration 570 :  1.5059828505060466
cost at iteration 571 :  1.5032792801168928
cost at iteration 572 :  1.500587510012929
cost at iteration 573 :  1.4979074867404847
cost at iteration 574 :  1.4952391571649495
cost at iteration 575 :  1.4925824684683713
cost at iteration 576 :  1.489937368147271
cost at iteration 577 :  1.4873038040103368
cost at iteration 578 :  1.4846817241761814
cost at iteration 579 :  1.4820710770710444
cost at iteration 580 :  1.4794718114267855
cost at iteration 581 :  1.4768838762784446
cost at iteration 582 :  1.474307220962263
cost at iteration 583 :  1.4717417951134668
cost at iteration 584 :  1.4691875486641068
cost at iteration 585 :  1.4666444318409806
cost at iteration 586 :  1.4641123951635793
cost at iteration 587 :  1.461591389441834
cost at iteration 588 :  1.4590813657743456
cost at iteration 589 :  1.4565822755460305
cost at iteration 590 :  1.45409407042634
cost at iteration 591 :  1.4516167023670545
cost at iteration 592 :  1.4491501236004325
cost at iteration 593 :  1.4466942866372101
cost at iteration 594 :  1.4442491442645735
cost at iteration 595 :  1.4418146495442976
cost at iteration 596 :  1.4393907558107133
cost at iteration 597 :  1.436977416668975
cost at iteration 598 :  1.4345745859929717
cost at iteration 599 :  1.4321822179235317
=================Mesh Quality Check=================
max quality: 1.1832106186704068  min quality: 0.5059467464043369  mean quality: 0.8036611401632984
===================================================

non linear with p#

check_mesh_quality()
c = pose_deformation_equation(p=2)

gfX = GridFunction(VEC)
gfX.Set((x *(1-x)*y*(1-y),0))

gfsetOld = GridFunction(VEC)
gfset = GridFunction(VEC)
gfset.Set((0,0))

scene = Draw(gfset,mesh,"gfset")
SetVisualization (deformation=True)


optimize_shape(linear=False)
mesh.SetDeformation(gfset)
check_mesh_quality()
=================Mesh Quality Check=================
max quality: 1.1832106186704068  min quality: 0.5059467464043369  mean quality: 0.8036611401632984
===================================================
Using p = 2
Cost at initial design 13.962388102754339
cost at iteration 0 :  13.962388102754339
Newton iteration  0
err =  37.410916722694864
Newton iteration  1
err =  1.8869214813442775e-13
cost at iteration 1 :  13.735607075299352
Newton iteration  0
err =  0.5046122613881028
Newton iteration  1
err =  4.4254487315379895e-14
cost at iteration 2 :  13.526734084141825
Newton iteration  0
err =  0.4939055405164836
Newton iteration  1
err =  4.363456293335879e-14
cost at iteration 3 :  13.3225169520256
Newton iteration  0
err =  0.48353678713810977
Newton iteration  1
err =  4.3915394551139703e-14
cost at iteration 4 :  13.122814774368264
Newton iteration  0
err =  0.473492208258463
Newton iteration  1
err =  4.218799185213006e-14
cost at iteration 5 :  12.9274922406713
Newton iteration  0
err =  0.46375867035191304
Newton iteration  1
err =  4.150894460032659e-14
cost at iteration 6 :  12.736419356284364
Newton iteration  0
err =  0.4543236642196425
Newton iteration  1
err =  4.182856003489478e-14
cost at iteration 7 :  12.54947118072345
Newton iteration  0
err =  0.4451752719111788
Newton iteration  1
err =  4.132104602605138e-14
cost at iteration 8 :  12.366527581410036
Newton iteration  0
err =  0.4363021355767377
Newton iteration  1
err =  4.0760320972295e-14
cost at iteration 9 :  12.187473001786444
Newton iteration  0
err =  0.42769342812216765
Newton iteration  1
err =  3.829579513147409e-14
cost at iteration 10 :  12.0121962428425
Newton iteration  0
err =  0.4193388255528737
Newton iteration  1
err =  3.9107369256031866e-14
cost at iteration 11 :  11.840590257157793
Newton iteration  0
err =  0.4112284808975276
Newton iteration  1
err =  3.829863816750436e-14
cost at iteration 12 :  11.672551954633821
Newton iteration  0
err =  0.4033529996112698
Newton iteration  1
err =  3.8760547806858766e-14
cost at iteration 13 :  11.50798201915008
Newton iteration  0
err =  0.3957034163676085
Newton iteration  1
err =  3.735585281847013e-14
cost at iteration 14 :  11.34678473543316
Newton iteration  0
err =  0.3882711731502286
Newton iteration  1
err =  3.679358962415122e-14
cost at iteration 15 :  11.188867825481541
Newton iteration  0
err =  0.3810480985658445
Newton iteration  1
err =  3.6355290237819265e-14
cost at iteration 16 :  11.034142293936283
Newton iteration  0
err =  0.3740263883041259
Newton iteration  1
err =  3.6423255470679594e-14
cost at iteration 17 :  10.882522281828033
Newton iteration  0
err =  0.3671985866726246
Newton iteration  1
err =  3.5168145324956277e-14
cost at iteration 18 :  10.733924928177329
Newton iteration  0
err =  0.3605575691446813
Newton iteration  1
err =  3.618354054421274e-14
cost at iteration 19 :  10.588270238956433
Newton iteration  0
err =  0.3540965258578649
Newton iteration  1
err =  3.466003726418036e-14
cost at iteration 20 :  10.445480962958058
Newton iteration  0
err =  0.3478089460072846
Newton iteration  1
err =  3.417259683553717e-14
cost at iteration 21 :  10.305482474149086
Newton iteration  0
err =  0.3416886030808486
Newton iteration  1
err =  3.363225970136936e-14
cost at iteration 22 :  10.1682026601132
Newton iteration  0
err =  0.3357295408871752
Newton iteration  1
err =  3.3657680969529875e-14
cost at iteration 23 :  10.03357181621658
Newton iteration  0
err =  0.32992606033011274
Newton iteration  1
err =  3.306700354588007e-14
cost at iteration 24 :  9.901522545154027
Newton iteration  0
err =  0.3242727068864402
Newton iteration  1
err =  3.3252305797305016e-14
cost at iteration 25 :  9.77198966155669
Newton iteration  0
err =  0.3187642587465984
Newton iteration  1
err =  3.226146710550975e-14
cost at iteration 26 :  9.644910101363886
Newton iteration  0
err =  0.3133957155809668
Newton iteration  1
err =  3.1820854747058744e-14
cost at iteration 27 :  9.520222835680451
Newton iteration  0
err =  0.3081622878946944
Newton iteration  1
err =  3.0166350286928206e-14
cost at iteration 28 :  9.397868788861103
Newton iteration  0
err =  0.30305938694005086
Newton iteration  1
err =  3.069483372461462e-14
cost at iteration 29 :  9.277790760578783
Newton iteration  0
err =  0.29808261515309536
Newton iteration  1
err =  3.1562027199591346e-14
cost at iteration 30 :  9.159933351650135
Newton iteration  0
err =  0.2932277570861597
Newton iteration  1
err =  3.1170655957994946e-14
cost at iteration 31 :  9.04424289340649
Newton iteration  0
err =  0.28849077080878344
Newton iteration  1
err =  3.098536667357062e-14
cost at iteration 32 :  8.930667380411426
Newton iteration  0
err =  0.28386777975081473
Newton iteration  1
err =  3.0822062436365796e-14
cost at iteration 33 :  8.819156406339548
Newton iteration  0
err =  0.27935506496283785
Newton iteration  1
err =  3.05640745024032e-14
cost at iteration 34 :  8.709661102842102
Newton iteration  0
err =  0.2749490577729587
Newton iteration  1
err =  2.9802068828511664e-14
cost at iteration 35 :  8.602134081235981
Newton iteration  0
err =  0.27064633281544587
Newton iteration  1
err =  2.9897951777283047e-14
cost at iteration 36 :  8.496529376863588
Newton iteration  0
err =  0.26644360141381834
Newton iteration  1
err =  2.979879115117831e-14
cost at iteration 37 :  8.392802395979752
Newton iteration  0
err =  0.2623377052972754
Newton iteration  1
err =  2.974981838288552e-14
cost at iteration 38 :  8.290909865030997
Newton iteration  0
err =  0.2583256106333532
Newton iteration  1
err =  2.916544314229156e-14
cost at iteration 39 :  8.190809782200082
Newton iteration  0
err =  0.2544044023596154
Newton iteration  1
err =  2.897880130158268e-14
cost at iteration 40 :  8.092461371098164
Newton iteration  0
err =  0.250571278798772
Newton iteration  1
err =  2.9149188223265397e-14
cost at iteration 41 :  7.995825036490887
Newton iteration  0
err =  0.24682354654094163
Newton iteration  1
err =  2.8388198077377854e-14
cost at iteration 42 :  7.900862321955695
Newton iteration  0
err =  0.243158615580995
Newton iteration  1
err =  2.927281966432458e-14
cost at iteration 43 :  7.8075358693684205
Newton iteration  0
err =  0.23957399469454663
Newton iteration  1
err =  2.8065210185213697e-14
cost at iteration 44 :  7.715809380129594
Newton iteration  0
err =  0.23606728704329613
Newton iteration  1
err =  2.7923198285731088e-14
cost at iteration 45 :  7.62564757803949
Newton iteration  0
err =  0.2326361859949322
Newton iteration  1
err =  2.770632745826125e-14
cost at iteration 46 :  7.5370161737411125
Newton iteration  0
err =  0.22927847114741998
Newton iteration  1
err =  2.7424721947950033e-14
cost at iteration 47 :  7.449881830653638
Newton iteration  0
err =  0.22599200454819124
Newton iteration  1
err =  2.6656721630787922e-14
cost at iteration 48 :  7.364212132320758
Newton iteration  0
err =  0.222774727095234
Newton iteration  1
err =  2.781901017996575e-14
cost at iteration 49 :  7.279975551107125
Newton iteration  0
err =  0.21962465511355989
Newton iteration  1
err =  2.5895490588667782e-14
cost at iteration 50 :  7.197141418175964
Newton iteration  0
err =  0.21653987709637457
Newton iteration  1
err =  2.697390467807852e-14
cost at iteration 51 :  7.115679894686386
Newton iteration  0
err =  0.2135185506024936
Newton iteration  1
err =  2.6024578110089555e-14
cost at iteration 52 :  7.0355619441525
Newton iteration  0
err =  0.21055889930284774
Newton iteration  1
err =  2.5894851503522644e-14
cost at iteration 53 :  6.956759305908962
Newton iteration  0
err =  0.20765921016731087
Newton iteration  1
err =  2.6162624102862953e-14
cost at iteration 54 :  6.879244469630437
Newton iteration  0
err =  0.20481783078554475
Newton iteration  1
err =  2.472819737989273e-14
cost at iteration 55 :  6.802990650856605
Newton iteration  0
err =  0.2020331668146619
Newton iteration  1
err =  2.4954669047009874e-14
cost at iteration 56 :  6.727971767475329
Newton iteration  0
err =  0.19930367954724418
Newton iteration  1
err =  2.5392850327199726e-14
cost at iteration 57 :  6.654162417119806
Newton iteration  0
err =  0.19662788359360964
Newton iteration  1
err =  2.4883511309256957e-14
cost at iteration 58 :  6.58153785543873
Newton iteration  0
err =  0.19400434467260577
Newton iteration  1
err =  2.495309677577878e-14
cost at iteration 59 :  6.510073975198968
Newton iteration  0
err =  0.19143167750582743
Newton iteration  1
err =  2.417306722744785e-14
cost at iteration 60 :  6.439747286183503
Newton iteration  0
err =  0.18890854380886424
Newton iteration  1
err =  2.395625253251977e-14
cost at iteration 61 :  6.370534895849263
Newton iteration  0
err =  0.1864336503764679
Newton iteration  1
err =  2.4179326833341218e-14
cost at iteration 62 :  6.302414490711011
Newton iteration  0
err =  0.18400574725527524
Newton iteration  1
err =  2.344351739420432e-14
cost at iteration 63 :  6.235364318419066
Newton iteration  0
err =  0.18162362600094725
Newton iteration  1
err =  2.2853774338649303e-14
cost at iteration 64 :  6.169363170500706
Newton iteration  0
err =  0.17928611801448188
Newton iteration  1
err =  2.3180929241673046e-14
cost at iteration 65 :  6.104390365736333
Newton iteration  0
err =  0.1769920929551006
Newton iteration  1
err =  2.29884149808342e-14
cost at iteration 66 :  6.040425734142831
Newton iteration  0
err =  0.1747404572242098
Newton iteration  1
err =  2.2407926063751413e-14
cost at iteration 67 :  5.977449601538489
Newton iteration  0
err =  0.1725301525187306
Newton iteration  1
err =  2.206037414049852e-14
cost at iteration 68 :  5.915442774664087
Newton iteration  0
err =  0.17036015444878166
Newton iteration  1
err =  2.2592987813551075e-14
cost at iteration 69 :  5.854386526837393
Newton iteration  0
err =  0.16822947121793563
Newton iteration  1
err =  2.2332218979478733e-14
cost at iteration 70 :  5.7942625841178
Newton iteration  0
err =  0.16613714236173774
Newton iteration  1
err =  2.3162700265990732e-14
cost at iteration 71 :  5.735053111960767
Newton iteration  0
err =  0.16408223754294976
Newton iteration  1
err =  2.0775037205834536e-14
cost at iteration 72 :  5.676740702340764
Newton iteration  0
err =  0.16206385539922852
Newton iteration  1
err =  2.0964340905126647e-14
cost at iteration 73 :  5.619308361324378
Newton iteration  0
err =  0.16008112244253891
Newton iteration  1
err =  2.0825835946906057e-14
cost at iteration 74 :  5.562739497074374
Newton iteration  0
err =  0.15813319200569742
Newton iteration  1
err =  2.0811534547971275e-14
cost at iteration 75 :  5.5070179082677875
Newton iteration  0
err =  0.15621924323612166
Newton iteration  1
err =  2.0992040522753272e-14
cost at iteration 76 :  5.452127772910877
Newton iteration  0
err =  0.15433848013196255
Newton iteration  1
err =  2.0233013255686562e-14
cost at iteration 77 :  5.3980536375354795
Newton iteration  0
err =  0.15249013062120986
Newton iteration  1
err =  2.1085335798391568e-14
cost at iteration 78 :  5.344780406760929
Newton iteration  0
err =  0.15067344567927646
Newton iteration  1
err =  2.031073308090293e-14
cost at iteration 79 :  5.292293333208101
Newton iteration  0
err =  0.14888769848541278
Newton iteration  1
err =  2.01752988900952e-14
cost at iteration 80 :  5.240578007750145
Newton iteration  0
err =  0.14713218361408853
Newton iteration  1
err =  2.0531484264300172e-14
cost at iteration 81 :  5.189620350088669
Newton iteration  0
err =  0.14540621626164266
Newton iteration  1
err =  2.0550330125789433e-14
cost at iteration 82 :  5.139406599640823
Newton iteration  0
err =  0.14370913150480172
Newton iteration  1
err =  1.9492766330431437e-14
cost at iteration 83 :  5.089923306726742
Newton iteration  0
err =  0.14204028359084603
Newton iteration  1
err =  1.988326561839075e-14
cost at iteration 84 :  5.041157324045196
Newton iteration  0
err =  0.14039904525719915
Newton iteration  1
err =  2.0092245142869654e-14
cost at iteration 85 :  4.993095798426067
Newton iteration  0
err =  0.1387848070794409
Newton iteration  1
err =  1.9595886835325316e-14
cost at iteration 86 :  4.945726162850187
Newton iteration  0
err =  0.13719697684644117
Newton iteration  1
err =  1.9221262179063385e-14
cost at iteration 87 :  4.899036128725188
Newton iteration  0
err =  0.13563497896062468
Newton iteration  1
err =  1.9633910333095527e-14
cost at iteration 88 :  4.853013678409138
Newton iteration  0
err =  0.1340982538638475
Newton iteration  1
err =  2.0035138301752925e-14
cost at iteration 89 :  4.80764705797135
Newton iteration  0
err =  0.13258625748560676
Newton iteration  1
err =  1.89933144326837e-14
cost at iteration 90 :  4.762924770182638
Newton iteration  0
err =  0.13109846071423814
Newton iteration  1
err =  1.9084821237066645e-14
cost at iteration 91 :  4.7188355677261935
Newton iteration  0
err =  0.12963434888950617
Newton iteration  1
err =  1.908089189428801e-14
cost at iteration 92 :  4.6753684466205785
Newton iteration  0
err =  0.1281934213149521
Newton iteration  1
err =  1.9183504421377437e-14
cost at iteration 93 :  4.632512639848179
Newton iteration  0
err =  0.12677519079062152
Newton iteration  1
err =  1.9424698294404588e-14
cost at iteration 94 :  4.590257611180182
Newton iteration  0
err =  0.12537918316337374
Newton iteration  1
err =  1.8321507594471017e-14
cost at iteration 95 :  4.548593049192574
Newton iteration  0
err =  0.12400493689564202
Newton iteration  1
err =  1.8534113542856063e-14
cost at iteration 96 :  4.507508861465317
Newton iteration  0
err =  0.12265200265086788
Newton iteration  1
err =  1.8902169458220728e-14
cost at iteration 97 :  4.466995168958002
Newton iteration  0
err =  0.12131994289514317
Newton iteration  1
err =  1.8254690088613713e-14
cost at iteration 98 :  4.427042300557031
Newton iteration  0
err =  0.12000833151442203
Newton iteration  1
err =  1.8367255808737127e-14
cost at iteration 99 :  4.3876407877865615
Newton iteration  0
err =  0.11871675344618546
Newton iteration  1
err =  1.7992305711889887e-14
cost at iteration 100 :  4.348781359678971
Newton iteration  0
err =  0.11744480432582347
Newton iteration  1
err =  1.884271770104458e-14
cost at iteration 101 :  4.310454937798219
Newton iteration  0
err =  0.11619209014621458
Newton iteration  1
err =  1.7319665579270483e-14
cost at iteration 102 :  4.2726526314112485
Newton iteration  0
err =  0.11495822692982222
Newton iteration  1
err =  1.774913269200734e-14
cost at iteration 103 :  4.235365732802734
Newton iteration  0
err =  0.1137428404140657
Newton iteration  1
err =  1.7829726332219462e-14
cost at iteration 104 :  4.1985857127274455
Newton iteration  0
err =  0.11254556574790923
Newton iteration  1
err =  1.7405712286700443e-14
cost at iteration 105 :  4.162304215995716
Newton iteration  0
err =  0.11136604719987281
Newton iteration  1
err =  1.756172201034928e-14
cost at iteration 106 :  4.126513057188112
Newton iteration  0
err =  0.1102039378767274
Newton iteration  1
err =  1.717652900010677e-14
cost at iteration 107 :  4.09120421649457
Newton iteration  0
err =  0.109058899453365
Newton iteration  1
err =  1.7852769246321188e-14
cost at iteration 108 :  4.0563698356731175
Newton iteration  0
err =  0.10793060191130138
Newton iteration  1
err =  1.746576537520936e-14
cost at iteration 109 :  4.0220022141255996
Newton iteration  0
err =  0.10681872328794062
Newton iteration  1
err =  1.742867705322332e-14
cost at iteration 110 :  3.9880938050851227
Newton iteration  0
err =  0.10572294943417872
Newton iteration  1
err =  1.676228308777403e-14
cost at iteration 111 :  3.954637211912643
Newton iteration  0
err =  0.10464297378118169
Newton iteration  1
err =  1.669486191878093e-14
cost at iteration 112 :  3.9216251844979313
Newton iteration  0
err =  0.10357849711504515
Newton iteration  1
err =  1.689193527853468e-14
cost at iteration 113 :  3.889050615763077
Newton iteration  0
err =  0.10252922736021403
Newton iteration  1
err =  1.7766500066388993e-14
cost at iteration 114 :  3.856906538263273
Newton iteration  0
err =  0.10149487936977013
Newton iteration  1
err =  1.6725253937344247e-14
cost at iteration 115 :  3.8251861208834375
Newton iteration  0
err =  0.10047517472378296
Newton iteration  1
err =  1.6238910010123773e-14
cost at iteration 116 :  3.793882665626482
Newton iteration  0
err =  0.09946984153432696
Newton iteration  1
err =  1.6395255617308234e-14
cost at iteration 117 :  3.7629896044906963
Newton iteration  0
err =  0.09847861425743942
Newton iteration  1
err =  1.5727642125647482e-14
cost at iteration 118 :  3.732500496433299
Newton iteration  0
err =  0.09750123351158986
Newton iteration  1
err =  1.568034726830412e-14
cost at iteration 119 :  3.7024090244177703
Newton iteration  0
err =  0.09653744590232963
Newton iteration  1
err =  1.6037676393262092e-14
cost at iteration 120 :  3.672708992541458
Newton iteration  0
err =  0.09558700385309694
Newton iteration  1
err =  1.6033422571830485e-14
cost at iteration 121 :  3.643394323241782
Newton iteration  0
err =  0.0946496654415284
Newton iteration  1
err =  1.6751463785325726e-14
cost at iteration 122 :  3.614459054578219
Newton iteration  0
err =  0.09372519424160372
Newton iteration  1
err =  1.5615148123479615e-14
cost at iteration 123 :  3.5858973375876624
Newton iteration  0
err =  0.09281335917098991
Newton iteration  1
err =  1.6064218745597007e-14
cost at iteration 124 :  3.557703433710705
Newton iteration  0
err =  0.0919139343430389
Newton iteration  1
err =  1.60652103742555e-14
cost at iteration 125 :  3.5298717122874064
Newton iteration  0
err =  0.09102669892492495
Newton iteration  1
err =  1.5693592792509927e-14
cost at iteration 126 :  3.502396648119054
Newton iteration  0
err =  0.09015143699881792
Newton iteration  1
err =  1.5608244602409452e-14
cost at iteration 127 :  3.475272819095457
Newton iteration  0
err =  0.08928793742898121
Newton iteration  1
err =  1.578886090981116e-14
cost at iteration 128 :  3.4484949038845554
Newton iteration  0
err =  0.08843599373253572
Newton iteration  1
err =  1.5148270305905005e-14
cost at iteration 129 :  3.4220576796829736
Newton iteration  0
err =  0.08759540395469989
Newton iteration  1
err =  1.5070675741909975e-14
cost at iteration 130 :  3.3959560200255274
Newton iteration  0
err =  0.08676597054772177
Newton iteration  1
err =  1.5696967849700388e-14
cost at iteration 131 :  3.3701848926520497
Newton iteration  0
err =  0.08594750025410273
Newton iteration  1
err =  1.500713706640934e-14
cost at iteration 132 :  3.3447393574297646
Newton iteration  0
err =  0.08513980399342592
Newton iteration  1
err =  1.516075667593954e-14
cost at iteration 133 :  3.319614564329166
Newton iteration  0
err =  0.08434269675235208
Newton iteration  1
err =  1.509716673750858e-14
cost at iteration 134 :  3.2948057514526643
Newton iteration  0
err =  0.0835559974789482
Newton iteration  1
err =  1.515976880432024e-14
cost at iteration 135 :  3.2703082431133517
Newton iteration  0
err =  0.08277952897943444
Newton iteration  1
err =  1.5225933734585967e-14
cost at iteration 136 :  3.2461174479632677
Newton iteration  0
err =  0.08201311781874558
Newton iteration  1
err =  1.4631550580765774e-14
cost at iteration 137 :  3.2222288571694198
Newton iteration  0
err =  0.08125659422386022
Newton iteration  1
err =  1.442773747625835e-14
cost at iteration 138 :  3.198638042635913
Newton iteration  0
err =  0.08050979199041783
Newton iteration  1
err =  1.3838707331107502e-14
cost at iteration 139 :  3.175340655271209
Newton iteration  0
err =  0.07977254839172224
Newton iteration  1
err =  1.4623982618863782e-14
cost at iteration 140 :  3.1523324232990286
Newton iteration  0
err =  0.07904470409115459
Newton iteration  1
err =  1.453715757656939e-14
cost at iteration 141 :  3.1296091506115546
Newton iteration  0
err =  0.07832610305660725
Newton iteration  1
err =  1.4349225516809257e-14
cost at iteration 142 :  3.107166715164028
Newton iteration  0
err =  0.07761659247786462
Newton iteration  1
err =  1.431283348533617e-14
cost at iteration 143 :  3.0850010674090296
Newton iteration  0
err =  0.07691602268648064
Newton iteration  1
err =  1.3941155891119381e-14
cost at iteration 144 :  3.063108228769916
Newton iteration  0
err =  0.07622424707793228
Newton iteration  1
err =  1.4540268373360353e-14
cost at iteration 145 :  3.041484290151681
Newton iteration  0
err =  0.07554112203586151
Newton iteration  1
err =  1.37695543053199e-14
cost at iteration 146 :  3.0201254104887747
Newton iteration  0
err =  0.07486650685916968
Newton iteration  1
err =  1.4098390029979167e-14
cost at iteration 147 :  2.9990278153282937
Newton iteration  0
err =  0.07420026369051555
Newton iteration  1
err =  1.352876507212318e-14
cost at iteration 148 :  2.9781877954480436
Newton iteration  0
err =  0.07354225744750793
Newton iteration  1
err =  1.3943105355031476e-14
cost at iteration 149 :  2.9576017055079973
Newton iteration  0
err =  0.07289235575546647
Newton iteration  1
err =  1.3870780890945365e-14
cost at iteration 150 :  2.9372659627345326
Newton iteration  0
err =  0.0722504288823758
Newton iteration  1
err =  1.443055284328891e-14
cost at iteration 151 :  2.917177045636685
Newton iteration  0
err =  0.07161634967566578
Newton iteration  1
err =  1.3677476582190203e-14
cost at iteration 152 :  2.897331492752833
Newton iteration  0
err =  0.07098999350057675
Newton iteration  1
err =  1.3224128636860656e-14
cost at iteration 153 :  2.8777259014281062
Newton iteration  0
err =  0.07037123818071718
Newton iteration  1
err =  1.3239281864153681e-14
cost at iteration 154 :  2.8583569266202615
Newton iteration  0
err =  0.06975996393977388
Newton iteration  1
err =  1.3366927330370637e-14
cost at iteration 155 :  2.839221279734652
Newton iteration  0
err =  0.06915605334540897
Newton iteration  1
err =  1.3027927935681592e-14
cost at iteration 156 :  2.8203157274863218
Newton iteration  0
err =  0.06855939125407003
Newton iteration  1
err =  1.3285580588119147e-14
cost at iteration 157 :  2.8016370907894563
Newton iteration  0
err =  0.06796986475810478
Newton iteration  1
err =  1.2929030979614705e-14
cost at iteration 158 :  2.783182243672615
Newton iteration  0
err =  0.06738736313372283
Newton iteration  1
err =  1.3232730305129686e-14
cost at iteration 159 :  2.7649481122195065
Newton iteration  0
err =  0.066811777790548
Newton iteration  1
err =  1.3319236468666535e-14
cost at iteration 160 :  2.746931673534717
Newton iteration  0
err =  0.06624300222282171
Newton iteration  1
err =  1.2468109900809142e-14
cost at iteration 161 :  2.7291299547333683
Newton iteration  0
err =  0.06568093196158645
Newton iteration  1
err =  1.2489869650503688e-14
cost at iteration 162 :  2.711540031954288
Newton iteration  0
err =  0.0651254645283338
Newton iteration  1
err =  1.208742928932281e-14
cost at iteration 163 :  2.694159029395877
Newton iteration  0
err =  0.06457649938984061
Newton iteration  1
err =  1.2145135797743977e-14
cost at iteration 164 :  2.6769841183745475
Newton iteration  0
err =  0.06403393791446017
Newton iteration  1
err =  1.2379827987839742e-14
cost at iteration 165 :  2.660012516404277
Newton iteration  0
err =  0.06349768332898342
Newton iteration  1
err =  1.1496994039170966e-14
cost at iteration 166 :  2.6432414862978044
Newton iteration  0
err =  0.06296764067752605
Newton iteration  1
err =  1.1795401191088332e-14
cost at iteration 167 :  2.6266683352878166
Newton iteration  0
err =  0.06244371678062058
Newton iteration  1
err =  1.1963292218033105e-14
cost at iteration 168 :  2.610290414168519
Newton iteration  0
err =  0.061925820195903813
Newton iteration  1
err =  1.19608416839995e-14
cost at iteration 169 :  2.5941051164563564
Newton iteration  0
err =  0.06141386117996197
Newton iteration  1
err =  1.127006141057721e-14
cost at iteration 170 :  2.578109877569801
Newton iteration  0
err =  0.060907751650451314
Newton iteration  1
err =  1.1685477546714063e-14
cost at iteration 171 :  2.562302174027734
Newton iteration  0
err =  0.060407405150140196
Newton iteration  1
err =  1.1870829775855695e-14
cost at iteration 172 :  2.5466795226655083
Newton iteration  0
err =  0.05991273681108933
Newton iteration  1
err =  1.1201846772749181e-14
cost at iteration 173 :  2.5312394798687916
Newton iteration  0
err =  0.05942366332028936
Newton iteration  1
err =  1.1982155855904087e-14
cost at iteration 174 :  2.515979640824251
Newton iteration  0
err =  0.058940102885943535
Newton iteration  1
err =  1.1523888841153917e-14
cost at iteration 175 :  2.5008976387870554
Newton iteration  0
err =  0.05846197520441013
Newton iteration  1
err =  1.1056230714241185e-14
cost at iteration 176 :  2.48599114436449
Newton iteration  0
err =  0.05798920142857383
Newton iteration  1
err =  1.1365928150840587e-14
cost at iteration 177 :  2.471257864815464
Newton iteration  0
err =  0.0575217041364129
Newton iteration  1
err =  1.1287123139500858e-14
cost at iteration 178 :  2.4566955433651088
Newton iteration  0
err =  0.05705940730062341
Newton iteration  1
err =  1.1209242436060726e-14
cost at iteration 179 :  2.4423019585348325
Newton iteration  0
err =  0.056602236259090166
Newton iteration  1
err =  1.1447131095083275e-14
cost at iteration 180 :  2.428074923486592
Newton iteration  0
err =  0.05615011768587242
Newton iteration  1
err =  1.1720226700249744e-14
cost at iteration 181 :  2.4140122853817636
Newton iteration  0
err =  0.055702979563093016
Newton iteration  1
err =  1.150235688536667e-14
cost at iteration 182 :  2.400111924753508
Newton iteration  0
err =  0.05526075115335434
Newton iteration  1
err =  1.1764797866550177e-14
cost at iteration 183 :  2.3863717548932826
Newton iteration  0
err =  0.05482336297315851
Newton iteration  1
err =  1.0966104965182652e-14
cost at iteration 184 :  2.372789721249856
Newton iteration  0
err =  0.054390746766426724
Newton iteration  1
err =  1.3886464952224779e-14
cost at iteration 185 :  2.359363800841905
Newton iteration  0
err =  0.053962835479271824
Newton iteration  1
err =  1.1007439335482281e-14
cost at iteration 186 :  2.346092001682705
Newton iteration  0
err =  0.0535395632349775
Newton iteration  1
err =  1.0910820269684659e-14
cost at iteration 187 :  2.3329723622172556
Newton iteration  0
err =  0.05312086530969652
Newton iteration  1
err =  1.1712142769008469e-14
cost at iteration 188 :  2.320002950771466
Newton iteration  0
err =  0.05270667810877874
Newton iteration  1
err =  1.0992817642057991e-14
cost at iteration 189 :  2.3071818650128186
Newton iteration  0
err =  0.052296939143755224
Newton iteration  1
err =  1.0568932073922597e-14
cost at iteration 190 :  2.2945072314226214
Newton iteration  0
err =  0.05189158700946137
Newton iteration  1
err =  1.07806256627934e-14
cost at iteration 191 :  2.281977204779218
Newton iteration  0
err =  0.051490561362246716
Newton iteration  1
err =  1.0602168922705796e-14
cost at iteration 192 :  2.269589967652175
Newton iteration  0
err =  0.05109380289845302
Newton iteration  1
err =  1.092523542843782e-14
cost at iteration 193 :  2.2573437299068657
Newton iteration  0
err =  0.05070125333329303
Newton iteration  1
err =  1.1203151189071521e-14
cost at iteration 194 :  2.2452367282195502
Newton iteration  0
err =  0.050312855380289456
Newton iteration  1
err =  1.1161432876548931e-14
cost at iteration 195 :  2.233267225602631
Newton iteration  0
err =  0.04992855273148209
Newton iteration  1
err =  1.079607730955441e-14
cost at iteration 196 :  2.2214335109393466
Newton iteration  0
err =  0.04954829003755455
Newton iteration  1
err =  1.0625111162586442e-14
cost at iteration 197 :  2.2097338985286035
Newton iteration  0
err =  0.049172012889007605
Newton iteration  1
err =  1.0826366877414574e-14
cost at iteration 198 :  2.1981667276386285
Newton iteration  0
err =  0.048799667797245996
Newton iteration  1
err =  1.0645821841132533e-14
cost at iteration 199 :  2.186730362070373
Newton iteration  0
err =  0.048431202176540865
Newton iteration  1
err =  1.0287885543556046e-14
cost at iteration 200 :  2.175423189729347
Newton iteration  0
err =  0.04806656432609454
Newton iteration  1
err =  1.0135520412293773e-14
cost at iteration 201 :  2.164243622206382
Newton iteration  0
err =  0.04770570341269443
Newton iteration  1
err =  1.0571375978387038e-14
cost at iteration 202 :  2.1531900943670617
Newton iteration  0
err =  0.04734856945355085
Newton iteration  1
err =  1.055526034513479e-14
cost at iteration 203 :  2.142261063949316
Newton iteration  0
err =  0.04699511329996871
Newton iteration  1
err =  1.052770201069458e-14
cost at iteration 204 :  2.1314550111690376
Newton iteration  0
err =  0.04664528662084675
Newton iteration  1
err =  1.027684725761626e-14
cost at iteration 205 :  2.120770438333968
Newton iteration  0
err =  0.04629904188682957
Newton iteration  1
err =  1.0186749844140971e-14
cost at iteration 206 :  2.110205869465021
Newton iteration  0
err =  0.045956332354955345
Newton iteration  1
err =  1.0220591737986611e-14
cost at iteration 207 :  2.0997598499254235
Newton iteration  0
err =  0.04561711205321229
Newton iteration  1
err =  1.0480088207610617e-14
cost at iteration 208 :  2.0894309460571536
Newton iteration  0
err =  0.04528133576600432
Newton iteration  1
err =  1.0308991655980906e-14
cost at iteration 209 :  2.079217744824544
Newton iteration  0
err =  0.04494895901937785
Newton iteration  1
err =  1.0055049983706952e-14
cost at iteration 210 :  2.0691188534652327
Newton iteration  0
err =  0.0446199380670377
Newton iteration  1
err =  1.0267110724762877e-14
cost at iteration 211 :  2.059132899147489
Newton iteration  0
err =  0.04429422987630728
Newton iteration  1
err =  1.04209815132137e-14
cost at iteration 212 :  2.0492585286351055
Newton iteration  0
err =  0.043971792114818345
Newton iteration  1
err =  9.799951514141768e-15
cost at iteration 213 :  2.0394944079579895
Newton iteration  0
err =  0.043652583136896204
Newton iteration  1
err =  1.0173706033227935e-14
cost at iteration 214 :  2.0298392220899695
Newton iteration  0
err =  0.04333656197090666
Newton iteration  1
err =  9.923679075918152e-15
cost at iteration 215 :  2.020291674632348
Newton iteration  0
err =  0.043023688306315984
Newton iteration  1
err =  1.0310929777184594e-14
cost at iteration 216 :  2.010850487504017
Newton iteration  0
err =  0.042713922481594425
Newton iteration  1
err =  9.86972975455831e-15
cost at iteration 217 :  2.001514400637221
Newton iteration  0
err =  0.042407225471618556
Newton iteration  1
err =  1.0228614191292706e-14
cost at iteration 218 :  1.9922821716796657
Newton iteration  0
err =  0.04210355887614734
Newton iteration  1
err =  9.571072201271871e-15
cost at iteration 219 :  1.9831525757021429
Newton iteration  0
err =  0.04180288490822253
Newton iteration  1
err =  9.961330197063449e-15
cost at iteration 220 :  1.974124404911329
Newton iteration  0
err =  0.04150516638236573
Newton iteration  1
err =  9.264800212373807e-15
cost at iteration 221 :  1.9651964683693186
Newton iteration  0
err =  0.04121036670397088
Newton iteration  1
err =  9.611506619708045e-15
cost at iteration 222 :  1.9563675917170427
Newton iteration  0
err =  0.04091844985799478
Newton iteration  1
err =  9.758634754387624e-15
cost at iteration 223 :  1.9476366169041894
Newton iteration  0
err =  0.040629380398669875
Newton iteration  1
err =  9.581598069818273e-15
cost at iteration 224 :  1.939002401923543
Newton iteration  0
err =  0.04034312343856659
Newton iteration  1
err =  9.84477477369681e-15
cost at iteration 225 :  1.9304638205506306
Newton iteration  0
err =  0.040059644638886754
Newton iteration  1
err =  9.56323526775078e-15
cost at iteration 226 :  1.9220197620883508
Newton iteration  0
err =  0.03977891019917632
Newton iteration  1
err =  9.406645109163497e-15
cost at iteration 227 :  1.9136691311161722
Newton iteration  0
err =  0.03950088684755654
Newton iteration  1
err =  9.550944053726653e-15
cost at iteration 228 :  1.9054108472444122
Newton iteration  0
err =  0.0392255418311759
Newton iteration  1
err =  9.509808278021026e-15
cost at iteration 229 :  1.8972438448727744
Newton iteration  0
err =  0.03895284290677868
Newton iteration  1
err =  9.543664633182348e-15
cost at iteration 230 :  1.8891670729536754
Newton iteration  0
err =  0.03868275833162445
Newton iteration  1
err =  9.392355158156096e-15
cost at iteration 231 :  1.8811794947598897
Newton iteration  0
err =  0.038415256854325106
Newton iteration  1
err =  9.19469262035993e-15
cost at iteration 232 :  1.8732800876563072
Newton iteration  0
err =  0.038150307706107354
Newton iteration  1
err =  9.353927558094818e-15
cost at iteration 233 :  1.8654678428763722
Newton iteration  0
err =  0.03788788059219876
Newton iteration  1
err =  9.364163939199157e-15
cost at iteration 234 :  1.857741765302336
Newton iteration  0
err =  0.037627945683314824
Newton iteration  1
err =  9.000083050109269e-15
cost at iteration 235 :  1.850100873249579
Newton iteration  0
err =  0.03737047360740218
Newton iteration  1
err =  9.161811975728316e-15
cost at iteration 236 :  1.8425441982549575
Newton iteration  0
err =  0.037115435441485496
Newton iteration  1
err =  9.054763167763328e-15
cost at iteration 237 :  1.83507078486923
Newton iteration  0
err =  0.03686280270382594
Newton iteration  1
err =  9.091233634866697e-15
cost at iteration 238 :  1.8276796904528105
Newton iteration  0
err =  0.0366125473459175
Newton iteration  1
err =  9.358502731421925e-15
cost at iteration 239 :  1.820369984975784
Newton iteration  0
err =  0.036364641745107965
Newton iteration  1
err =  9.02674909724339e-15
cost at iteration 240 :  1.8131407508211574
Newton iteration  0
err =  0.03611905869692779
Newton iteration  1
err =  9.243068807294554e-15
cost at iteration 241 :  1.8059910825919052
Newton iteration  0
err =  0.03587577140784971
Newton iteration  1
err =  9.120282639027886e-15
cost at iteration 242 :  1.798920086921396
Newton iteration  0
err =  0.0356347534880395
Newton iteration  1
err =  8.907400695484304e-15
cost at iteration 243 :  1.791926882287321
Newton iteration  0
err =  0.035395978944427274
Newton iteration  1
err =  9.116842223798918e-15
cost at iteration 244 :  1.7850105988289253
Newton iteration  0
err =  0.035159422173700335
Newton iteration  1
err =  9.776157539729783e-15
cost at iteration 245 :  1.7781703781674647
Newton iteration  0
err =  0.03492505795557075
Newton iteration  1
err =  9.040293483995416e-15
cost at iteration 246 :  1.7714053732299668
Newton iteration  0
err =  0.03469286144619063
Newton iteration  1
err =  8.854989880527818e-15
cost at iteration 247 :  1.7647147480762109
Newton iteration  0
err =  0.034462808171624636
Newton iteration  1
err =  8.698419961648744e-15
cost at iteration 248 :  1.75809767772856
Newton iteration  0
err =  0.03423487402141801
Newton iteration  1
err =  8.821326787631025e-15
cost at iteration 249 :  1.7515533480049628
Newton iteration  0
err =  0.03400903524246598
Newton iteration  1
err =  8.79791680765444e-15
cost at iteration 250 :  1.7450809553549127
Newton iteration  0
err =  0.03378526843285553
Newton iteration  1
err =  8.754063350856767e-15
cost at iteration 251 :  1.738679706698288
Newton iteration  0
err =  0.03356355053579524
Newton iteration  1
err =  9.600195281082473e-15
cost at iteration 252 :  1.73234881926702
Newton iteration  0
err =  0.033343858833751294
Newton iteration  1
err =  8.62249460571416e-15
cost at iteration 253 :  1.7260875204496053
Newton iteration  0
err =  0.03312617094281825
Newton iteration  1
err =  8.913762178379869e-15
cost at iteration 254 :  1.7198950476381134
Newton iteration  0
err =  0.03291046480679048
Newton iteration  1
err =  8.2545569559591e-15
cost at iteration 255 :  1.7137706480782255
Newton iteration  0
err =  0.03269671869188567
Newton iteration  1
err =  9.309215277754185e-15
cost at iteration 256 :  1.707713578721636
Newton iteration  0
err =  0.03248491118110761
Newton iteration  1
err =  8.742127662734057e-15
cost at iteration 257 :  1.7017231060809876
Newton iteration  0
err =  0.03227502116901751
Newton iteration  1
err =  8.869286219743448e-15
cost at iteration 258 :  1.6957985060875984
Newton iteration  0
err =  0.03206702785637755
Newton iteration  1
err =  8.875557449390921e-15
cost at iteration 259 :  1.689939063951304
Newton iteration  0
err =  0.03186091074517047
Newton iteration  1
err =  8.533144650340584e-15
cost at iteration 260 :  1.684144074022993
Newton iteration  0
err =  0.03165664963344586
Newton iteration  1
err =  8.558013309963941e-15
cost at iteration 261 :  1.678412839659345
Newton iteration  0
err =  0.03145422461037338
Newton iteration  1
err =  8.676350916376211e-15
cost at iteration 262 :  1.6727446730899125
Newton iteration  0
err =  0.03125361605145948
Newton iteration  1
err =  8.390347964326866e-15
cost at iteration 263 :  1.667138895286503
Newton iteration  0
err =  0.031054804613713714
Newton iteration  1
err =  8.5360704219807e-15
cost at iteration 264 :  1.6615948358349548
Newton iteration  0
err =  0.03085777123106561
Newton iteration  1
err =  8.265887461544966e-15
cost at iteration 265 :  1.656111832808618
Newton iteration  0
err =  0.030662497109662193
Newton iteration  1
err =  8.331200950328902e-15
cost at iteration 266 :  1.6506892326446538
Newton iteration  0
err =  0.030468963723464165
Newton iteration  1
err =  9.095647199825274e-15
cost at iteration 267 :  1.6453263900217099
Newton iteration  0
err =  0.030277152809807112
Newton iteration  1
err =  8.14405257501335e-15
cost at iteration 268 :  1.640022667740313
Newton iteration  0
err =  0.03008704636505082
Newton iteration  1
err =  8.296304556449326e-15
cost at iteration 269 :  1.6347774366049042
Newton iteration  0
err =  0.029898626640313064
Newton iteration  1
err =  8.458986810645183e-15
cost at iteration 270 :  1.6295900753078367
Newton iteration  0
err =  0.029711876137422703
Newton iteration  1
err =  8.351842636506941e-15
cost at iteration 271 :  1.6244599703155504
Newton iteration  0
err =  0.029526777604571237
Newton iteration  1
err =  8.139835503796986e-15
cost at iteration 272 :  1.6193865157564487
Newton iteration  0
err =  0.029343314032563468
Newton iteration  1
err =  7.885385115291997e-15
cost at iteration 273 :  1.6143691133107332
Newton iteration  0
err =  0.029161468650673324
Newton iteration  1
err =  8.33638732154356e-15
cost at iteration 274 :  1.6094071721021148
Newton iteration  0
err =  0.02898122492285875
Newton iteration  1
err =  8.269826229598571e-15
cost at iteration 275 :  1.604500108591133
Newton iteration  0
err =  0.028802566543947545
Newton iteration  1
err =  8.058224588326022e-15
cost at iteration 276 :  1.5996473464705097
Newton iteration  0
err =  0.02862547743589727
Newton iteration  1
err =  9.045948371734885e-15
cost at iteration 277 :  1.5948483165619856
Newton iteration  0
err =  0.02844994174405817
Newton iteration  1
err =  8.176254215301458e-15
cost at iteration 278 :  1.5901024567150657
Newton iteration  0
err =  0.02827594383370081
Newton iteration  1
err =  8.240728358040853e-15
cost at iteration 279 :  1.5854092117072285
Newton iteration  0
err =  0.028103468286348143
Newton iteration  1
err =  8.391086223953496e-15
cost at iteration 280 :  1.580768033145884
Newton iteration  0
err =  0.027932499896412267
Newton iteration  1
err =  8.557067579853645e-15
cost at iteration 281 :  1.5761783793719932
Newton iteration  0
err =  0.02776302366763554
Newton iteration  1
err =  7.737330875101282e-15
cost at iteration 282 :  1.571639715365117
Newton iteration  0
err =  0.027595024809915766
Newton iteration  1
err =  8.125209069618717e-15
cost at iteration 283 :  1.5671515126500843
Newton iteration  0
err =  0.027428488735891616
Newton iteration  1
err =  8.026627691704212e-15
cost at iteration 284 :  1.5627132492051914
Newton iteration  0
err =  0.027263401057694883
Newton iteration  1
err =  7.78343667748406e-15
cost at iteration 285 :  1.55832440937178
Newton iteration  0
err =  0.027099747583913612
Newton iteration  1
err =  7.982513934726638e-15
cost at iteration 286 :  1.5539844837655796
Newton iteration  0
err =  0.026937514316254164
Newton iteration  1
err =  7.973795284697961e-15
cost at iteration 287 :  1.549692969189044
Newton iteration  0
err =  0.02677668744670462
Newton iteration  1
err =  7.48007636062614e-15
cost at iteration 288 :  1.5454493685452957
Newton iteration  0
err =  0.026617253354342545
Newton iteration  1
err =  7.500205247459698e-15
cost at iteration 289 :  1.5412531907537554
Newton iteration  0
err =  0.02645919860240618
Newton iteration  1
err =  7.650441019754136e-15
cost at iteration 290 :  1.5371039506663529
Newton iteration  0
err =  0.02630250993554459
Newton iteration  1
err =  7.583058943236667e-15
cost at iteration 291 :  1.5330011689861016
Newton iteration  0
err =  0.026147174276699287
Newton iteration  1
err =  7.664425478913346e-15
cost at iteration 292 :  1.528944372185966
Newton iteration  0
err =  0.025993178724475798
Newton iteration  1
err =  7.575901710533628e-15
cost at iteration 293 :  1.52493309242979
Newton iteration  0
err =  0.025840510550363294
Newton iteration  1
err =  8.050068489905822e-15
cost at iteration 294 :  1.5209668674939323
Newton iteration  0
err =  0.025689157195950975
Newton iteration  1
err =  7.882893167940142e-15
cost at iteration 295 :  1.5170452406905635
Newton iteration  0
err =  0.025539106270275315
Newton iteration  1
err =  7.366113700053655e-15
cost at iteration 296 :  1.5131677607918528
Newton iteration  0
err =  0.025390345547250854
Newton iteration  1
err =  7.540186047190117e-15
cost at iteration 297 :  1.5093339819555065
Newton iteration  0
err =  0.025242862963100743
Newton iteration  1
err =  7.20292447477048e-15
cost at iteration 298 :  1.505543463651462
Newton iteration  0
err =  0.02509664661367184
Newton iteration  1
err =  7.435562457103207e-15
cost at iteration 299 :  1.5017957705896965
Newton iteration  0
err =  0.024951684752153428
Newton iteration  1
err =  7.584314153186356e-15
cost at iteration 300 :  1.4980904726491875
Newton iteration  0
err =  0.024807965786468485
Newton iteration  1
err =  7.623211761175861e-15
cost at iteration 301 :  1.4944271448079767
Newton iteration  0
err =  0.024665478276945887
Newton iteration  1
err =  7.534415310744736e-15
cost at iteration 302 :  1.4908053670742984
Newton iteration  0
err =  0.024524210933878448
Newton iteration  1
err =  6.948741615770957e-15
cost at iteration 303 :  1.487224724418807
Newton iteration  0
err =  0.024384152615336496
Newton iteration  1
err =  7.647145393369557e-15
cost at iteration 304 :  1.4836848067078714
Newton iteration  0
err =  0.024245292324699415
Newton iteration  1
err =  7.06546624158858e-15
cost at iteration 305 :  1.4801852086377938
Newton iteration  0
err =  0.02410761920848713
Newton iteration  1
err =  7.53735802647261e-15
cost at iteration 306 :  1.4767255296702095
Newton iteration  0
err =  0.02397112255423721
Newton iteration  1
err =  7.283318303577721e-15
cost at iteration 307 :  1.4733053739682442
Newton iteration  0
err =  0.023835791788056202
Newton iteration  1
err =  7.333658832930132e-15
cost at iteration 308 :  1.4699243503340491
Newton iteration  0
err =  0.02370161647281716
Newton iteration  1
err =  7.160977070905204e-15
cost at iteration 309 :  1.4665820721467613
Newton iteration  0
err =  0.023568586305822124
Newton iteration  1
err =  7.180923128721638e-15
cost at iteration 310 :  1.463278157301896
Newton iteration  0
err =  0.023436691116734555
Newton iteration  1
err =  7.452080377693612e-15
cost at iteration 311 :  1.4600122281512942
Newton iteration  0
err =  0.023305920865706663
Newton iteration  1
err =  7.469386492665296e-15
cost at iteration 312 :  1.456783911444226
Newton iteration  0
err =  0.023176265641149874
Newton iteration  1
err =  7.05563947776917e-15
cost at iteration 313 :  1.453592838269358
Newton iteration  0
err =  0.023047715657897542
Newton iteration  1
err =  6.9702950709158405e-15
cost at iteration 314 :  1.4504386439974586
Newton iteration  0
err =  0.022920261255288895
Newton iteration  1
err =  7.067074135898416e-15
cost at iteration 315 :  1.4473209682250379
Newton iteration  0
err =  0.02279389289509277
Newton iteration  1
err =  7.359202477254275e-15
cost at iteration 316 :  1.4442394547189723
Newton iteration  0
err =  0.022668601159807355
Newton iteration  1
err =  6.9680686706595326e-15
cost at iteration 317 :  1.4411937513617123
Newton iteration  0
err =  0.022544376750721073
Newton iteration  1
err =  7.818626527741644e-15
cost at iteration 318 :  1.4381835100975193
Newton iteration  0
err =  0.02242121048603246
Newton iteration  1
err =  6.868513216054526e-15
cost at iteration 319 :  1.4352083868794991
Newton iteration  0
err =  0.02229909329921659
Newton iteration  1
err =  6.9981385995630444e-15
cost at iteration 320 :  1.4322680416171023
Newton iteration  0
err =  0.02217801623707167
Newton iteration  1
err =  6.809975678551943e-15
cost at iteration 321 :  1.4293621381249835
Newton iteration  0
err =  0.022057970458135755
Newton iteration  1
err =  7.039104936825374e-15
cost at iteration 322 :  1.4264903440721564
Newton iteration  0
err =  0.02193894723085272
Newton iteration  1
err =  6.874515387267434e-15
cost at iteration 323 :  1.4236523309318958
Newton iteration  0
err =  0.021820937931965635
Newton iteration  1
err =  6.8161884695552965e-15
cost at iteration 324 :  1.4208477739328556
Newton iteration  0
err =  0.02170393404488064
Newton iteration  1
err =  6.946049796171718e-15
cost at iteration 325 :  1.4180763520103388
Newton iteration  0
err =  0.021587927157880753
Newton iteration  1
err =  7.020565052423586e-15
cost at iteration 326 :  1.4153377477586393
Newton iteration  0
err =  0.02147290896275567
Newton iteration  1
err =  6.906553854623312e-15
cost at iteration 327 :  1.4126316473839353
Newton iteration  0
err =  0.02135887125293732
Newton iteration  1
err =  7.086819471574583e-15
cost at iteration 328 :  1.4099577406579167
Newton iteration  0
err =  0.021245805922218336
Newton iteration  1
err =  6.789760753033439e-15
cost at iteration 329 :  1.4073157208721752
Newton iteration  0
err =  0.021133704963023296
Newton iteration  1
err =  6.626435611997565e-15
cost at iteration 330 :  1.4047052847930914
Newton iteration  0
err =  0.021022560464935443
Newton iteration  1
err =  7.079174804108771e-15
cost at iteration 331 :  1.402126132617442
Newton iteration  0
err =  0.020912364613301416
Newton iteration  1
err =  7.596051467262828e-15
cost at iteration 332 :  1.3995779679288523
Newton iteration  0
err =  0.02080310968764293
Newton iteration  1
err =  6.485422564876489e-15
cost at iteration 333 :  1.3970604976545442
Newton iteration  0
err =  0.02069478806022344
Newton iteration  1
err =  6.2259718238701865e-15
cost at iteration 334 :  1.394573432022856
Newton iteration  0
err =  0.02058739219473534
Newton iteration  1
err =  6.4775913418277475e-15
cost at iteration 335 :  1.3921164845214715
Newton iteration  0
err =  0.020480914644782924
Newton iteration  1
err =  6.7253109615051366e-15
cost at iteration 336 :  1.3896893718562897
Newton iteration  0
err =  0.0203753480525228
Newton iteration  1
err =  6.418206625657877e-15
cost at iteration 337 :  1.3872918139105188
Newton iteration  0
err =  0.020270685147378193
Newton iteration  1
err =  6.531838168622758e-15
cost at iteration 338 :  1.3849235337048131
Newton iteration  0
err =  0.020166918744588175
Newton iteration  1
err =  6.591120509306196e-15
cost at iteration 339 :  1.3825842573577827
Newton iteration  0
err =  0.020064041743884434
Newton iteration  1
err =  6.7771590611222375e-15
cost at iteration 340 :  1.3802737140468042
Newton iteration  0
err =  0.019962047128375904
Newton iteration  1
err =  6.417868904569248e-15
cost at iteration 341 :  1.3779916359700686
Newton iteration  0
err =  0.019860927963027058
Newton iteration  1
err =  6.6640964111746985e-15
cost at iteration 342 :  1.3757377583084465
Newton iteration  0
err =  0.019760677393508196
Newton iteration  1
err =  6.40250117841675e-15
cost at iteration 343 :  1.3735118191882918
Newton iteration  0
err =  0.019661288644921077
Newton iteration  1
err =  6.814382063579104e-15
cost at iteration 344 :  1.3713135596446613
Newton iteration  0
err =  0.01956275502068379
Newton iteration  1
err =  6.185956616481845e-15
cost at iteration 345 :  1.3691427235852867
Newton iteration  0
err =  0.01946506990107182
Newton iteration  1
err =  6.267419111187413e-15
cost at iteration 346 :  1.3669990577544866
Newton iteration  0
err =  0.01936822674232329
Newton iteration  1
err =  6.534867497082318e-15
cost at iteration 347 :  1.364882311698321
Newton iteration  0
err =  0.01927221907513693
Newton iteration  1
err =  6.756303440668222e-15
cost at iteration 348 :  1.3627922377295518
Newton iteration  0
err =  0.019177040503908484
Newton iteration  1
err =  6.33259031543098e-15
cost at iteration 349 :  1.36072859089365
Newton iteration  0
err =  0.019082684705179116
Newton iteration  1
err =  6.229324677804126e-15
cost at iteration 350 :  1.3586911289349992
Newton iteration  0
err =  0.018989145426779777
Newton iteration  1
err =  6.4555105148366385e-15
cost at iteration 351 :  1.3566796122634308
Newton iteration  0
err =  0.018896416486651316
Newton iteration  1
err =  6.2273718513285985e-15
cost at iteration 352 :  1.3546938039217367
Newton iteration  0
err =  0.018804491771654424
Newton iteration  1
err =  6.3108726943394105e-15
cost at iteration 353 :  1.3527334695529398
Newton iteration  0
err =  0.018713365236733608
Newton iteration  1
err =  7.233120106997587e-15
cost at iteration 354 :  1.3507983773687198
Newton iteration  0
err =  0.018623030903518654
Newton iteration  1
err =  6.404668438865247e-15
cost at iteration 355 :  1.3488882981177963
Newton iteration  0
err =  0.018533482859504503
Newton iteration  1
err =  5.969444456880012e-15
cost at iteration 356 :  1.3470030050548205
Newton iteration  0
err =  0.018444715257035078
Newton iteration  1
err =  6.122444158604868e-15
cost at iteration 357 :  1.3451422739099563
Newton iteration  0
err =  0.018356722312132482
Newton iteration  1
err =  5.892235518465987e-15
cost at iteration 358 :  1.3433058828586582
Newton iteration  0
err =  0.018269498303551086
Newton iteration  1
err =  6.173646213943639e-15
cost at iteration 359 :  1.3414936124918218
Newton iteration  0
err =  0.018183037571826315
Newton iteration  1
err =  5.810765128175706e-15
cost at iteration 360 :  1.339705245786437
Newton iteration  0
err =  0.01809733451835703
Newton iteration  1
err =  5.989111736222585e-15
cost at iteration 361 :  1.337940568076825
Newton iteration  0
err =  0.018012383604240397
Newton iteration  1
err =  5.945514451287301e-15
cost at iteration 362 :  1.3361993670259684
Newton iteration  0
err =  0.017928179349425
Newton iteration  1
err =  6.0940408806413e-15
cost at iteration 363 :  1.3344814325972343
Newton iteration  0
err =  0.01784471633191367
Newton iteration  1
err =  5.851655961315726e-15
cost at iteration 364 :  1.332786557026774
Newton iteration  0
err =  0.017761989186568128
Newton iteration  1
err =  6.001268321584067e-15
cost at iteration 365 :  1.3311145347960986
Newton iteration  0
err =  0.017679992604465528
Newton iteration  1
err =  6.369560087098946e-15
cost at iteration 366 :  1.3294651626050498
Newton iteration  0
err =  0.01759872133172245
Newton iteration  1
err =  6.168926911786663e-15
cost at iteration 367 :  1.3278382393450783
Newton iteration  0
err =  0.017518170168895984
Newton iteration  1
err =  5.90289735754596e-15
cost at iteration 368 :  1.3262335660729636
Newton iteration  0
err =  0.017438333969929377
Newton iteration  1
err =  6.174426911674652e-15
cost at iteration 369 :  1.324650945985037
Newton iteration  0
err =  0.017359207641296143
Newton iteration  1
err =  5.9948213860156785e-15
cost at iteration 370 :  1.3230901843914091
Newton iteration  0
err =  0.017280786141140293
Newton iteration  1
err =  5.79212216110854e-15
cost at iteration 371 :  1.321551088690667
Newton iteration  0
err =  0.017203064478691313
Newton iteration  1
err =  5.874699895966492e-15
cost at iteration 372 :  1.3200334683454054
Newton iteration  0
err =  0.0171260377128738
Newton iteration  1
err =  6.0502263532975194e-15
cost at iteration 373 :  1.3185371348568984
Newton iteration  0
err =  0.017049700952271974
Newton iteration  1
err =  6.223870960019374e-15
cost at iteration 374 :  1.3170619017417504
Newton iteration  0
err =  0.01697404935348647
Newton iteration  1
err =  6.043434207164425e-15
cost at iteration 375 :  1.3156075845072315
Newton iteration  0
err =  0.016899078121071037
Newton iteration  1
err =  5.943256304764971e-15
cost at iteration 376 :  1.314174000628201
Newton iteration  0
err =  0.016824782506145792
Newton iteration  1
err =  6.844394568230098e-15
cost at iteration 377 :  1.312760969523343
Newton iteration  0
err =  0.016751157806269276
Newton iteration  1
err =  6.8495683193937266e-15
cost at iteration 378 :  1.3113683125328477
Newton iteration  0
err =  0.01667819936392945
Newton iteration  1
err =  5.923953893807216e-15
cost at iteration 379 :  1.3099958528950681
Newton iteration  0
err =  0.016605902566455295
Newton iteration  1
err =  6.011100033495453e-15
cost at iteration 380 :  1.3086434157245443
Newton iteration  0
err =  0.01653426284487813
Newton iteration  1
err =  5.8669752003181395e-15
cost at iteration 381 :  1.307310827989854
Newton iteration  0
err =  0.016463275673377385
Newton iteration  1
err =  5.812744711672109e-15
cost at iteration 382 :  1.3059979184918404
Newton iteration  0
err =  0.01639293656853612
Newton iteration  1
err =  6.09228836665353e-15
cost at iteration 383 :  1.3047045178421026
Newton iteration  0
err =  0.016323241088515208
Newton iteration  1
err =  5.837020910162792e-15
cost at iteration 384 :  1.303430458441755
Newton iteration  0
err =  0.01625418483256933
Newton iteration  1
err =  5.9253014304586185e-15
cost at iteration 385 :  1.3021755744606078
Newton iteration  0
err =  0.016185763440107907
Newton iteration  1
err =  5.864652114231764e-15
cost at iteration 386 :  1.3009397018163489
Newton iteration  0
err =  0.016117972590316394
Newton iteration  1
err =  6.177141125724994e-15
cost at iteration 387 :  1.2997226781543867
Newton iteration  0
err =  0.016050808001115413
Newton iteration  1
err =  6.179931120366633e-15
cost at iteration 388 :  1.2985243428274889
Newton iteration  0
err =  0.01598426542889171
Newton iteration  1
err =  5.848886668157695e-15
cost at iteration 389 :  1.2973445368761696
Newton iteration  0
err =  0.015918340667485306
Newton iteration  1
err =  5.981263418983771e-15
cost at iteration 390 :  1.2961831030087758
Newton iteration  0
err =  0.015853029547899314
Newton iteration  1
err =  6.107898601837273e-15
cost at iteration 391 :  1.2950398855825656
Newton iteration  0
err =  0.015788327937311185
Newton iteration  1
err =  6.319386942874652e-15
cost at iteration 392 :  1.2939147305842142
Newton iteration  0
err =  0.01572423173879727
Newton iteration  1
err =  5.707167852131487e-15
cost at iteration 393 :  1.292807485611282
Newton iteration  0
err =  0.01566073689047814
Newton iteration  1
err =  5.822553018982528e-15
cost at iteration 394 :  1.2917179998534898
Newton iteration  0
err =  0.01559783936491258
Newton iteration  1
err =  5.570841773342408e-15
cost at iteration 395 :  1.290646124074247
Newton iteration  0
err =  0.015535535168737662
Newton iteration  1
err =  5.529300688326916e-15
cost at iteration 396 :  1.2895917105927512
Newton iteration  0
err =  0.015473820341799003
Newton iteration  1
err =  5.883189559447527e-15
cost at iteration 397 :  1.288554613265949
Newton iteration  0
err =  0.015412690956789921
Newton iteration  1
err =  5.8028068681042725e-15
cost at iteration 398 :  1.2875346874709503
Newton iteration  0
err =  0.015352143118500092
Newton iteration  1
err =  5.995074323107259e-15
cost at iteration 399 :  1.2865317900875928
Newton iteration  0
err =  0.01529217296345702
Newton iteration  1
err =  5.992024095095047e-15
cost at iteration 400 :  1.2855457794811842
Newton iteration  0
err =  0.01523277665917158
Newton iteration  1
err =  5.690705783984125e-15
cost at iteration 401 :  1.2845765154855913
Newton iteration  0
err =  0.015173950403808362
Newton iteration  1
err =  6.134458827159431e-15
cost at iteration 402 :  1.2836238593863976
Newton iteration  0
err =  0.015115690425328728
Newton iteration  1
err =  6.042845223990916e-15
cost at iteration 403 :  1.282687673904337
Newton iteration  0
err =  0.015057992981438884
Newton iteration  1
err =  5.897063683198683e-15
cost at iteration 404 :  1.2817678231791796
Newton iteration  0
err =  0.015000854358606215
Newton iteration  1
err =  5.966105664034391e-15
cost at iteration 405 :  1.2808641727531263
Newton iteration  0
err =  0.0149442708718234
Newton iteration  1
err =  5.659846629868067e-15
cost at iteration 406 :  1.2799765895554072
Newton iteration  0
err =  0.01488823886407971
Newton iteration  1
err =  5.842821371365607e-15
cost at iteration 407 :  1.2791049418861868
Newton iteration  0
err =  0.014832754705723657
Newton iteration  1
err =  5.5748200987261026e-15
cost at iteration 408 :  1.278249099401141
Newton iteration  0
err =  0.014777814794148999
Newton iteration  1
err =  5.690878976044175e-15
cost at iteration 409 :  1.277408933096219
Newton iteration  0
err =  0.014723415553248423
Newton iteration  1
err =  5.8441405195437625e-15
cost at iteration 410 :  1.2765843152922904
Newton iteration  0
err =  0.01466955343293476
Newton iteration  1
err =  5.45332165413058e-15
cost at iteration 411 :  1.2757751196203786
Newton iteration  0
err =  0.014616224908662075
Newton iteration  1
err =  5.5056489564586646e-15
cost at iteration 412 :  1.2749812210068245
Newton iteration  0
err =  0.014563426481126194
Newton iteration  1
err =  5.722646853019058e-15
cost at iteration 413 :  1.2742024956588198
Newton iteration  0
err =  0.01451115467551712
Newton iteration  1
err =  5.6274982105043215e-15
cost at iteration 414 :  1.2734388210497436
Newton iteration  0
err =  0.014459406041409398
Newton iteration  1
err =  6.202802125666374e-15
cost at iteration 415 :  1.2726900759052038
Newton iteration  0
err =  0.014408177152097254
Newton iteration  1
err =  5.770546981147802e-15
cost at iteration 416 :  1.2719561401888462
Newton iteration  0
err =  0.014357464604424672
Newton iteration  1
err =  5.793726147004763e-15
cost at iteration 417 :  1.2712368950887096
Newton iteration  0
err =  0.014307265017978453
Newton iteration  1
err =  5.948249531102655e-15
cost at iteration 418 :  1.2705322230031904
Newton iteration  0
err =  0.014257575035155875
Newton iteration  1
err =  5.549556453513066e-15
cost at iteration 419 :  1.269842007527918
Newton iteration  0
err =  0.014208391320224072
Newton iteration  1
err =  5.879297618825031e-15
cost at iteration 420 :  1.269166133441925
Newton iteration  0
err =  0.014159710559761859
Newton iteration  1
err =  5.7452734608892574e-15
cost at iteration 421 :  1.268504486695106
Newton iteration  0
err =  0.014111529461078945
Newton iteration  1
err =  6.3896739378319985e-15
cost at iteration 422 :  1.2678569543944442
Newton iteration  0
err =  0.014063844753122319
Newton iteration  1
err =  6.171689722003069e-15
cost at iteration 423 :  1.267223424791868
Newton iteration  0
err =  0.014016653185020015
Newton iteration  1
err =  5.366509694857744e-15
cost at iteration 424 :  1.266603787270908
Newton iteration  0
err =  0.013969951526594223
Newton iteration  1
err =  5.762464790329768e-15
cost at iteration 425 :  1.2659979323346957
Newton iteration  0
err =  0.013923736567252171
Newton iteration  1
err =  6.186034317632932e-15
cost at iteration 426 :  1.2654057515931374
Newton iteration  0
err =  0.013878005116458519
Newton iteration  1
err =  5.8697645829827826e-15
cost at iteration 427 :  1.2648271377511338
Newton iteration  0
err =  0.013832754002271762
Newton iteration  1
err =  5.336173758415826e-15
cost at iteration 428 :  1.264261984595809
Newton iteration  0
err =  0.01378798007244223
Newton iteration  1
err =  5.503303333346419e-15
cost at iteration 429 :  1.2637101869853793
Newton iteration  0
err =  0.013743680192794048
Newton iteration  1
err =  5.5886762406567975e-15
cost at iteration 430 :  1.263171640836775
Newton iteration  0
err =  0.013699851247552672
Newton iteration  1
err =  5.451590079248483e-15
cost at iteration 431 :  1.262646243114162
Newton iteration  0
err =  0.013656490139066701
Newton iteration  1
err =  6.1455988648969106e-15
cost at iteration 432 :  1.2621338918174283
Newton iteration  0
err =  0.013613593787267599
Newton iteration  1
err =  5.738012305029323e-15
cost at iteration 433 :  1.261634485970965
Newton iteration  0
err =  0.013571159129389066
Newton iteration  1
err =  5.778859492840211e-15
cost at iteration 434 :  1.261147925612145
Newton iteration  0
err =  0.013529183119917122
Newton iteration  1
err =  5.923486417032145e-15
cost at iteration 435 :  1.260674111780469
Newton iteration  0
err =  0.013487662730018233
Newton iteration  1
err =  5.559852202221101e-15
cost at iteration 436 :  1.2602129465064973
Newton iteration  0
err =  0.013446594947409438
Newton iteration  1
err =  5.6562147892789476e-15
cost at iteration 437 :  1.259764332801013
Newton iteration  0
err =  0.013405976776103916
Newton iteration  1
err =  5.703763191652325e-15
cost at iteration 438 :  1.259328174644364
Newton iteration  0
err =  0.01336580523616169
Newton iteration  1
err =  5.817800866376618e-15
cost at iteration 439 :  1.258904376975956
Newton iteration  0
err =  0.013326077363218803
Newton iteration  1
err =  6.3299934867846884e-15
cost at iteration 440 :  1.2584928456836433
Newton iteration  0
err =  0.013286790208554361
Newton iteration  1
err =  5.5787674553857555e-15
cost at iteration 441 :  1.2580934875935779
Newton iteration  0
err =  0.013247940838603353
Newton iteration  1
err =  5.875665028954993e-15
cost at iteration 442 :  1.257706210459827
Newton iteration  0
err =  0.013209526335043256
Newton iteration  1
err =  6.308463397658376e-15
cost at iteration 443 :  1.2573309229547052
Newton iteration  0
err =  0.013171543793982148
Newton iteration  1
err =  5.359965503124107e-15
cost at iteration 444 :  1.2569675346583193
Newton iteration  0
err =  0.013133990326247156
Newton iteration  1
err =  5.288284530563216e-15
cost at iteration 445 :  1.2566159560489947
Newton iteration  0
err =  0.013096863057138083
Newton iteration  1
err =  5.414375694774223e-15
cost at iteration 446 :  1.2562760984936177
Newton iteration  0
err =  0.013060159125939493
Newton iteration  1
err =  5.366710623576179e-15
cost at iteration 447 :  1.2559478742379209
Newton iteration  0
err =  0.01302387568592076
Newton iteration  1
err =  5.828671748546536e-15
cost at iteration 448 :  1.2556311963971858
Newton iteration  0
err =  0.012988009903911482
Newton iteration  1
err =  5.6671927280742446e-15
cost at iteration 449 :  1.2553259789464075
Newton iteration  0
err =  0.012952558960567578
Newton iteration  1
err =  5.560721988244578e-15
cost at iteration 450 :  1.2550321367116926
Newton iteration  0
err =  0.012917520049776685
Newton iteration  1
err =  5.534973803749219e-15
cost at iteration 451 :  1.254749585360689
Newton iteration  0
err =  0.012882890378465092
Newton iteration  1
err =  5.1384203831663886e-15
cost at iteration 452 :  1.2544782413936248
Newton iteration  0
err =  0.012848667166694545
Newton iteration  1
err =  5.523569549800466e-15
cost at iteration 453 :  1.2542180221344141
Newton iteration  0
err =  0.012814847647228254
Newton iteration  1
err =  5.545068107144237e-15
cost at iteration 454 :  1.2539688457217764
Newton iteration  0
err =  0.012781429065740818
Newton iteration  1
err =  5.625002827229345e-15
cost at iteration 455 :  1.2537306311006549
Newton iteration  0
err =  0.012748408680023575
Newton iteration  1
err =  5.382164647065205e-15
cost at iteration 456 :  1.2535032980134027
Newton iteration  0
err =  0.012715783760563001
Newton iteration  1
err =  5.626101902185502e-15
cost at iteration 457 :  1.253286766991362
Newton iteration  0
err =  0.012683551589885943
Newton iteration  1
err =  5.297123799678775e-15
cost at iteration 458 :  1.2530809593465972
Newton iteration  0
err =  0.012651709462571898
Newton iteration  1
err =  5.247514321301302e-15
cost at iteration 459 :  1.2528857971632268
Newton iteration  0
err =  0.012620254685198413
Newton iteration  1
err =  5.1864529907606165e-15
cost at iteration 460 :  1.2527012032895979
Newton iteration  0
err =  0.012589184575982121
Newton iteration  1
err =  5.929417633523389e-15
cost at iteration 461 :  1.2525271013297758
Newton iteration  0
err =  0.012558496465086446
Newton iteration  1
err =  5.350340246053722e-15
cost at iteration 462 :  1.2523634156358718
Newton iteration  0
err =  0.012528187693897547
Newton iteration  1
err =  5.184389037328501e-15
cost at iteration 463 :  1.2522100712998747
Newton iteration  0
err =  0.012498255615436701
Newton iteration  1
err =  5.319012884650349e-15
cost at iteration 464 :  1.2520669941458866
Newton iteration  0
err =  0.01246869759407648
Newton iteration  1
err =  5.560518466604486e-15
cost at iteration 465 :  1.2519341107224065
Newton iteration  0
err =  0.012439511005175144
Newton iteration  1
err =  6.059595639270101e-15
cost at iteration 466 :  1.2518113482945947
Newton iteration  0
err =  0.012410693235454653
Newton iteration  1
err =  5.520391388566282e-15
cost at iteration 467 :  1.251698634836755
Newton iteration  0
err =  0.012382241682422003
Newton iteration  1
err =  5.3865678087956365e-15
cost at iteration 468 :  1.2515958990247538
Newton iteration  0
err =  0.01235415375482515
Newton iteration  1
err =  5.2858908212564884e-15
cost at iteration 469 :  1.2515030702288537
Newton iteration  0
err =  0.01232642687193126
Newton iteration  1
err =  6.327437380811941e-15
cost at iteration 470 :  1.2514200785061478
Newton iteration  0
err =  0.01229905846382115
Newton iteration  1
err =  5.499855409901974e-15
cost at iteration 471 :  1.2513468545934314
Newton iteration  0
err =  0.012272045971458015
Newton iteration  1
err =  5.175984348766833e-15
cost at iteration 472 :  1.25128332990008
Newton iteration  0
err =  0.012245386846256773
Newton iteration  1
err =  6.083251342207983e-15
cost at iteration 473 :  1.2512294365010628
Newton iteration  0
err =  0.012219078550117295
Newton iteration  1
err =  5.80007850161522e-15
cost at iteration 474 :  1.2511851071297697
Newton iteration  0
err =  0.012193118555581094
Newton iteration  1
err =  5.45350643293988e-15
cost at iteration 475 :  1.2511502751713044
Newton iteration  0
err =  0.01216750434524448
Newton iteration  1
err =  5.961584381602905e-15
cost at iteration 476 :  1.2511248746554224
Newton iteration  0
err =  0.012142233412599776
Newton iteration  1
err =  5.273484203131633e-15
cost at iteration 477 :  1.2511088402501607
Newton iteration  0
err =  0.012117303260865664
Newton iteration  1
err =  5.1859008661746714e-15
cost at iteration 478 :  1.2511021072547472
Newton iteration  0
err =  0.012092711403952024
Newton iteration  1
err =  5.212613895612797e-15
scale =  5e-05
New cost during line search:  1.24416530401301
cost at iteration 479 :  1.24416530401301
Newton iteration  0
err =  0.006034500361259126
Newton iteration  1
err =  5.1850160002518524e-15
scale =  5e-05
New cost during line search:  1.2441630501393042
cost at iteration 480 :  1.2441630501393042
Newton iteration  0
err =  0.006028510531511427
Newton iteration  1
err =  5.699703000782396e-15
scale =  5e-05
New cost during line search:  1.2441630851978775
scale =  2.5e-05
New cost during line search:  1.2441629892777155
cost at iteration 481 :  1.2441629892777155
Newton iteration  0
err =  0.003011349663225307
Newton iteration  1
err =  5.3918315175923764e-15
scale =  5e-05
New cost during line search:  1.244164165163216
scale =  2.5e-05
New cost during line search:  1.2441634989728694
scale =  1.25e-05
New cost during line search:  1.2441632245647072
scale =  6.25e-06
New cost during line search:  1.244163102031241
scale =  3.125e-06
New cost during line search:  1.2441630444319927
scale =  1.5625e-06
New cost during line search:  1.2441630165492428
scale =  7.8125e-07
New cost during line search:  1.2441630028370567
scale =  3.90625e-07
New cost during line search:  1.244162996038274
scale =  1.953125e-07
New cost during line search:  1.2441629926532174
scale =  9.765625e-08
New cost during line search:  1.2441629909642733
scale =  4.8828125e-08
New cost during line search:  1.2441629901206601
scale =  2.44140625e-08
New cost during line search:  1.24416298969912
scale =  1.220703125e-08
New cost during line search:  1.2441629894884019
scale =  6.103515625e-09
New cost during line search:  1.2441629893830144
scale =  3.0517578125e-09
New cost during line search:  1.2441629893303696
scale =  1.52587890625e-09
New cost during line search:  1.244162989304043
scale =  7.62939453125e-10
New cost during line search:  1.2441629892908286
scale =  3.814697265625e-10
New cost during line search:  1.2441629892842996
scale =  1.9073486328125e-10
New cost during line search:  1.2441629892809771
scale =  9.5367431640625e-11
New cost during line search:  1.2441629892793518
scale =  4.76837158203125e-11
New cost during line search:  1.2441629892785286
scale =  2.384185791015625e-11
New cost during line search:  1.2441629892780874
scale =  1.1920928955078126e-11
New cost during line search:  1.2441629892778983
scale =  5.960464477539063e-12
New cost during line search:  1.2441629892777857
scale =  2.9802322387695314e-12
New cost during line search:  1.2441629892777466
scale =  1.4901161193847657e-12
New cost during line search:  1.2441629892777295
scale =  7.450580596923828e-13
No more descent can be found
=================Mesh Quality Check=================
max quality: 1.441807871875525  min quality: 0.5167710698008765  mean quality: 0.8594784827643153
===================================================
check_mesh_quality()
c = pose_deformation_equation(p=5)

gfX = GridFunction(VEC)
gfX.Set((x *(1-x)*y*(1-y),0))

gfset = GridFunction(VEC)
gfsetOld = GridFunction(VEC)
gfset.Set((0,0))

scene = Draw(gfset,mesh,"gfset")
SetVisualization (deformation=True)

optimize_shape(linear=False)
mesh.SetDeformation(gfset)
check_mesh_quality()
=================Mesh Quality Check=================
max quality: 1.441807871875525  min quality: 0.5167710698008765  mean quality: 0.8594784827643153
===================================================
Using p = 5
Cost at initial design 13.962388102754325
cost at iteration 0 :  13.962388102754325
Newton iteration  0
err =  69.38374121332043
Newton iteration  1
err =  95525973.66348527
Newton iteration  2
err =  46534454.951519154
Newton iteration  3
err =  22668761.32833453
Newton iteration  4
err =  11042844.290232865
Newton iteration  5
err =  5379403.322999564
Newton iteration  6
err =  2620518.71335221
Newton iteration  7
err =  1276557.6244199935
Newton iteration  8
err =  621861.3742348586
Newton iteration  9
err =  302933.1079457248
Newton iteration  10
err =  147570.61879367885
Newton iteration  11
err =  71887.44606638905
Newton iteration  12
err =  35019.19889199272
Newton iteration  13
err =  17059.226994771438
Newton iteration  14
err =  8310.219017450478
Newton iteration  15
err =  4048.233680696883
Newton iteration  16
err =  1972.0533601981629
Newton iteration  17
err =  960.6643252697774
Newton iteration  18
err =  467.97698908025376
Newton iteration  19
err =  227.9697209864447
Newton iteration  20
err =  111.05296347167638
Newton iteration  21
err =  54.098079001127424
Newton iteration  22
err =  26.35274210207727
Newton iteration  23
err =  12.837171975807319
Newton iteration  24
err =  6.256405109294945
Newton iteration  25
err =  3.055553881474144
Newton iteration  26
err =  1.5032775702643142
Newton iteration  27
err =  0.7539379142719076
Newton iteration  28
err =  0.39516650093499783
Newton iteration  29
err =  0.1962113835173666
Newton iteration  30
err =  0.07616340191361548
Newton iteration  31
err =  0.0161952102187628
Newton iteration  32
err =  0.0008527401283532236
Newton iteration  33
err =  2.4418088889929785e-06
Newton iteration  34
err =  2.0062203907033034e-11
cost at iteration 1 :  13.934752179253156
Newton iteration  0
err =  0.1653834014651807
Newton iteration  1
err =  3.979479133693186e-05
Newton iteration  2
err =  7.957871044945422e-10
cost at iteration 2 :  13.908869188124667
Newton iteration  0
err =  0.16007838567877328
Newton iteration  1
err =  3.870229556976139e-05
Newton iteration  2
err =  7.584259436602396e-10
cost at iteration 3 :  13.883178024567632
Newton iteration  0
err =  0.15494615756371
Newton iteration  1
err =  3.766084712112133e-05
Newton iteration  2
err =  7.234984402115556e-10
cost at iteration 4 :  13.857668858567187
Newton iteration  0
err =  0.1499811640330529
Newton iteration  1
err =  3.666793450350137e-05
Newton iteration  2
err =  6.908339011113341e-10
cost at iteration 5 :  13.832332481890504
Newton iteration  0
err =  0.14517803283706313
Newton iteration  1
err =  3.572121220495211e-05
Newton iteration  2
err =  6.60263167729356e-10
cost at iteration 6 :  13.807160268538782
Newton iteration  0
err =  0.14053156671882644
Newton iteration  1
err =  3.4818486783809284e-05
Newton iteration  2
err =  6.316453715935145e-10
cost at iteration 7 :  13.78214413771727
Newton iteration  0
err =  0.13603673776060105
Newton iteration  1
err =  3.395770417769166e-05
Newton iteration  2
err =  6.048447608968938e-10
cost at iteration 8 :  13.757276519160792
Newton iteration  0
err =  0.13168868191436386
Newton iteration  1
err =  3.3136938089904465e-05
Newton iteration  2
err =  5.797341565124422e-10
cost at iteration 9 :  13.732550320666476
Newton iteration  0
err =  0.12748269371029505
Newton iteration  1
err =  3.235437933179608e-05
Newton iteration  2
err =  5.56202405922325e-10
cost at iteration 10 :  13.707958897692258
Newton iteration  0
err =  0.12341422113652271
Newton iteration  1
err =  3.1608326096894855e-05
Newton iteration  2
err =  5.341329099250298e-10
cost at iteration 11 :  13.683496024889928
Newton iteration  0
err =  0.11947886068575622
Newton iteration  1
err =  3.089717506403407e-05
Newton iteration  2
err =  5.134375212418207e-10
cost at iteration 12 :  13.659155869449927
Newton iteration  0
err =  0.11567235256220998
Newton iteration  1
err =  3.0219413283596676e-05
Newton iteration  2
err =  4.940149573669752e-10
cost at iteration 13 :  13.634932966142163
Newton iteration  0
err =  0.11199057604318337
Newton iteration  1
err =  2.9573610736510203e-05
Newton iteration  2
err =  4.757887483157533e-10
cost at iteration 14 :  13.61082219394503
Newton iteration  0
err =  0.1084295449911515
Newton iteration  1
err =  2.895841360895716e-05
Newton iteration  2
err =  4.58674688496704e-10
cost at iteration 15 :  13.586818754161266
Newton iteration  0
err =  0.10498540351029195
Newton iteration  1
err =  2.8372538117719427e-05
Newton iteration  2
err =  4.426029538964038e-10
cost at iteration 16 :  13.562918149926984
Newton iteration  0
err =  0.10165442174235445
Newton iteration  1
err =  2.781476491709504e-05
Newton iteration  2
err =  4.2750353825139356e-10
cost at iteration 17 :  13.539116167023822
Newton iteration  0
err =  0.09843299179885555
Newton iteration  1
err =  2.7283934032305238e-05
Newton iteration  2
err =  4.1331902806538765e-10
cost at iteration 18 :  13.515408855912371
Newton iteration  0
err =  0.09531762382206924
Newton iteration  1
err =  2.6778940249599286e-05
Newton iteration  2
err =  3.9998511048991195e-10
cost at iteration 19 :  13.491792514908903
Newton iteration  0
err =  0.09230494217315402
Newton iteration  1
err =  2.629872896108289e-05
Newton iteration  2
err =  3.8745411602311227e-10
cost at iteration 20 :  13.468263674432137
Newton iteration  0
err =  0.08939168174115018
Newton iteration  1
err =  2.5842292424781e-05
Newton iteration  2
err =  3.756707515319509e-10
cost at iteration 21 :  13.444819082252364
Newton iteration  0
err =  0.0865746843694662
Newton iteration  1
err =  2.540866636047218e-05
Newton iteration  2
err =  3.6459457775500474e-10
cost at iteration 22 :  13.42145568967985
Newton iteration  0
err =  0.08385089539543973
Newton iteration  1
err =  2.499692694405459e-05
Newton iteration  2
err =  3.5418250913762986e-10
cost at iteration 23 :  13.398170638630528
Newton iteration  0
err =  0.08121736029902628
Newton iteration  1
err =  2.4606188084728796e-05
Newton iteration  2
err =  3.443927347695955e-10
cost at iteration 24 :  13.374961249516515
Newton iteration  0
err =  0.07867122145659185
Newton iteration  1
err =  2.423559898811694e-05
Newton iteration  2
err =  3.35190857824499e-10
cost at iteration 25 :  13.351825009905221
Newton iteration  0
err =  0.07620971499641337
Newton iteration  1
err =  2.3884342009685574e-05
Newton iteration  2
err =  3.2654130299985124e-10
cost at iteration 26 :  13.328759563901423
Newton iteration  0
err =  0.07383016775139907
Newton iteration  1
err =  2.355163072972866e-05
Newton iteration  2
err =  3.184129010365336e-10
cost at iteration 27 :  13.305762702204017
Newton iteration  0
err =  0.07152999430656994
Newton iteration  1
err =  2.323670824209527e-05
Newton iteration  2
err =  3.107789920303499e-10
cost at iteration 28 :  13.282832352795541
Newton iteration  0
err =  0.06930669413638757
Newton iteration  1
err =  2.2938845657968435e-05
Newton iteration  2
err =  3.0361176175778246e-10
cost at iteration 29 :  13.259966572224714
Newton iteration  0
err =  0.06715784883000046
Newton iteration  1
err =  2.265734076610297e-05
Newton iteration  2
err =  2.9688810519160065e-10
cost at iteration 30 :  13.237163537442639
Newton iteration  0
err =  0.06508111939966152
Newton iteration  1
err =  2.2391516862928474e-05
Newton iteration  2
err =  2.9058246355070095e-10
cost at iteration 31 :  13.214421538158547
Newton iteration  0
err =  0.06307424367048843
Newton iteration  1
err =  2.2140721702169e-05
Newton iteration  2
err =  2.846745086947417e-10
cost at iteration 32 :  13.191738969682694
Newton iteration  0
err =  0.06113503374691637
Newton iteration  1
err =  2.1904326588968548e-05
Newton iteration  2
err =  2.7914525544891714e-10
cost at iteration 33 :  13.169114326223509
Newton iteration  0
err =  0.059261373554205375
Newton iteration  1
err =  2.168172555486305e-05
Newton iteration  2
err =  2.7397597230784513e-10
cost at iteration 34 :  13.146546194611872
Newton iteration  0
err =  0.057451216450950256
Newton iteration  1
err =  2.147233463200448e-05
Newton iteration  2
err =  2.691545187931426e-10
cost at iteration 35 :  13.124033248424368
Newton iteration  0
err =  0.05570258291042372
Newton iteration  1
err =  2.1275591218018307e-05
Newton iteration  2
err =  2.646617110668048e-10
cost at iteration 36 :  13.10157424248032
Newton iteration  0
err =  0.05401355826715598
Newton iteration  1
err =  2.109095345805131e-05
Newton iteration  2
err =  2.6048540761035893e-10
cost at iteration 37 :  13.079168007689233
Newton iteration  0
err =  0.05238229052713169
Newton iteration  1
err =  2.091789974569797e-05
Newton iteration  2
err =  2.5660795548213085e-10
cost at iteration 38 :  13.056813446225696
Newton iteration  0
err =  0.05080698823785021
Newton iteration  1
err =  2.075592820393548e-05
Newton iteration  2
err =  2.530263506901833e-10
cost at iteration 39 :  13.034509527011236
Newton iteration  0
err =  0.04928591841650187
Newton iteration  1
err =  2.060455623027167e-05
Newton iteration  2
err =  2.497234385687985e-10
cost at iteration 40 :  13.012255281484032
Newton iteration  0
err =  0.04781740453345619
Newton iteration  1
err =  2.0463320062343084e-05
Newton iteration  2
err =  2.466916017938375e-10
cost at iteration 41 :  12.990049799637204
Newton iteration  0
err =  0.04639982454906527
Newton iteration  1
err =  2.033177436745505e-05
Newton iteration  2
err =  2.439211066824284e-10
cost at iteration 42 :  12.967892226309438
Newton iteration  0
err =  0.045031609000942496
Newton iteration  1
err =  2.0209491822957732e-05
Newton iteration  2
err =  2.414044578616109e-10
cost at iteration 43 :  12.94578175771133
Newton iteration  0
err =  0.043711239140922356
Newton iteration  1
err =  2.0096062729260756e-05
Newton iteration  2
err =  2.3913171029803674e-10
cost at iteration 44 :  12.923717638173049
Newton iteration  0
err =  0.042437245117982594
Newton iteration  1
err =  1.9991094615979194e-05
Newton iteration  2
err =  2.370986581939119e-10
cost at iteration 45 :  12.90169915709876
Newton iteration  0
err =  0.04120820420733661
Newton iteration  1
err =  1.9894211849599213e-05
Newton iteration  2
err =  2.352992442704463e-10
cost at iteration 46 :  12.879725646114599
Newton iteration  0
err =  0.040022739082168385
Newton iteration  1
err =  1.980505523952205e-05
Newton iteration  2
err =  2.337256606303884e-10
cost at iteration 47 :  12.857796476398315
Newton iteration  0
err =  0.038879516128408585
Newton iteration  1
err =  1.972328164685503e-05
Newton iteration  2
err =  2.3237397221252913e-10
cost at iteration 48 :  12.835911056178656
Newton iteration  0
err =  0.03777724379966409
Newton iteration  1
err =  1.9648563580626895e-05
Newton iteration  2
err =  2.312400645822846e-10
cost at iteration 49 :  12.814068828393879
Newton iteration  0
err =  0.03671467101220227
Newton iteration  1
err =  1.9580588802687717e-05
Newton iteration  2
err =  2.3031968943383468e-10
cost at iteration 50 :  12.792269268499494
Newton iteration  0
err =  0.035690585579110964
Newton iteration  1
err =  1.951905991753995e-05
Newton iteration  2
err =  2.2960584151286877e-10
cost at iteration 51 :  12.770511882414597
Newton iteration  0
err =  0.034703812682531816
Newton iteration  1
err =  1.946369395871884e-05
Newton iteration  2
err =  2.2909982075148822e-10
cost at iteration 52 :  12.748796204600286
Newton iteration  0
err =  0.03375321338396836
Newton iteration  1
err =  1.9414221986048657e-05
Newton iteration  2
err =  2.2879579784684218e-10
cost at iteration 53 :  12.727121796258643
Newton iteration  0
err =  0.03283768317226119
Newton iteration  1
err =  1.9370388659518247e-05
Newton iteration  2
err =  2.2869239362203937e-10
cost at iteration 54 :  12.705488243647784
Newton iteration  0
err =  0.03195615054965868
Newton iteration  1
err =  1.933195183103412e-05
Newton iteration  2
err =  2.2878511469500165e-10
cost at iteration 55 :  12.68389515650325
Newton iteration  0
err =  0.031107575655835402
Newton iteration  1
err =  1.929868212607858e-05
Newton iteration  2
err =  2.2907324036470654e-10
cost at iteration 56 :  12.662342166560158
Newton iteration  0
err =  0.030290948931083875
Newton iteration  1
err =  1.9270362523484163e-05
Newton iteration  2
err =  2.2955564201463845e-10
cost at iteration 57 :  12.640828926169023
Newton iteration  0
err =  0.029505289818631777
Newton iteration  1
err =  1.9246787951226814e-05
Newton iteration  2
err =  2.3023049095460936e-10
cost at iteration 58 :  12.61935510700028
Newton iteration  0
err =  0.028749645508260254
Newton iteration  1
err =  1.922776486537214e-05
Newton iteration  2
err =  2.3109646983083572e-10
cost at iteration 59 :  12.597920398830784
Newton iteration  0
err =  0.028023089721468292
Newton iteration  1
err =  1.9213110858595057e-05
Newton iteration  2
err =  2.321511859820282e-10
cost at iteration 60 :  12.576524508407468
Newton iteration  0
err =  0.02732472154018269
Newton iteration  1
err =  1.9202654246036426e-05
Newton iteration  2
err =  2.3339566618893463e-10
cost at iteration 61 :  12.555167158383666
Newton iteration  0
err =  0.02665366428086308
Newton iteration  1
err =  1.9196233687256903e-05
Newton iteration  2
err =  2.348285220768461e-10
cost at iteration 62 :  12.533848086322804
Newton iteration  0
err =  0.026009064415096463
Newton iteration  1
err =  1.9193697788336936e-05
Newton iteration  2
err =  2.364489416968931e-10
cost at iteration 63 :  12.51256704376589
Newton iteration  0
err =  0.025390090539169886
Newton iteration  1
err =  1.9194904734163153e-05
Newton iteration  2
err =  2.3825846491152083e-10
cost at iteration 64 :  12.491323795357815
Newton iteration  0
err =  0.024795932394166
Newton iteration  1
err =  1.9199721918646004e-05
Newton iteration  2
err =  2.4025729305966295e-10
cost at iteration 65 :  12.470118118029832
Newton iteration  0
err =  0.02422579993877051
Newton iteration  1
err =  1.9208025591171783e-05
Newton iteration  2
err =  2.424435362353477e-10
cost at iteration 66 :  12.44894980023375
Newton iteration  0
err =  0.0236789224761036
Newton iteration  1
err =  1.9219700503408023e-05
Newton iteration  2
err =  2.448217031206039e-10
cost at iteration 67 :  12.42781864122484
Newton iteration  0
err =  0.02315454783665432
Newton iteration  1
err =  1.923463958612323e-05
Newton iteration  2
err =  2.4738942506964383e-10
cost at iteration 68 :  12.40672445039091
Newton iteration  0
err =  0.022651941618807476
Newton iteration  1
err =  1.925274361801653e-05
Newton iteration  2
err =  2.501497563538057e-10
cost at iteration 69 :  12.385667046623333
Newton iteration  0
err =  0.022170386487614553
Newton iteration  1
err =  1.9273920921481275e-05
Newton iteration  2
err =  2.5310335673416403e-10
cost at iteration 70 :  12.364646257728754
Newton iteration  0
err =  0.021709181533457703
Newton iteration  1
err =  1.929808705961451e-05
Newton iteration  2
err =  2.5625381777137316e-10
cost at iteration 71 :  12.343661919878256
Newton iteration  0
err =  0.0212676416899205
Newton iteration  1
err =  1.9325164555610053e-05
Newton iteration  2
err =  2.5959944349798417e-10
cost at iteration 72 :  12.32271387709122
Newton iteration  0
err =  0.020845097212369715
Newton iteration  1
err =  1.935508262371922e-05
Newton iteration  2
err =  2.631453977900381e-10
cost at iteration 73 :  12.301801980752408
Newton iteration  0
err =  0.020440893215118897
Newton iteration  1
err =  1.9387776896846784e-05
Newton iteration  2
err =  2.6689592041981143e-10
cost at iteration 74 :  12.280926089159916
Newton iteration  0
err =  0.020054389267094404
Newton iteration  1
err =  1.9423189205114426e-05
Newton iteration  2
err =  2.708494746066684e-10
cost at iteration 75 :  12.260086067101575
Newton iteration  0
err =  0.019684959044516368
Newton iteration  1
err =  1.9461267314617638e-05
Newton iteration  2
err =  2.7501388383764224e-10
cost at iteration 76 :  12.23928178545892
Newton iteration  0
err =  0.019331990037441543
Newton iteration  1
err =  1.9501964732748676e-05
Newton iteration  2
err =  2.793915292194201e-10
cost at iteration 77 :  12.21851312083586
Newton iteration  0
err =  0.018994883308412872
Newton iteration  1
err =  1.9545240486624467e-05
Newton iteration  2
err =  2.839843674446678e-10
cost at iteration 78 :  12.197779955211471
Newton iteration  0
err =  0.018673053299673954
Newton iteration  1
err =  1.9591058936352454e-05
Newton iteration  2
err =  2.887981910456866e-10
cost at iteration 79 :  12.17708217561512
Newton iteration  0
err =  0.018365927684582145
Newton iteration  1
err =  1.9639389591711877e-05
Newton iteration  2
err =  2.938406026635228e-10
cost at iteration 80 :  12.15641967382163
Newton iteration  0
err =  0.018072947259868247
Newton iteration  1
err =  1.9690206947247617e-05
Newton iteration  2
err =  2.991119047573253e-10
cost at iteration 81 :  12.135792346066657
Newton iteration  0
err =  0.017793565873736518
Newton iteration  1
err =  1.9743490326403185e-05
Newton iteration  2
err =  3.0462088217597053e-10
cost at iteration 82 :  12.115200092780002
Newton iteration  0
err =  0.017527250384397196
Newton iteration  1
err =  1.9799223733302733e-05
Newton iteration  2
err =  3.1037277445294874e-10
cost at iteration 83 :  12.094642818336116
Newton iteration  0
err =  0.01727348064523154
Newton iteration  1
err =  1.9857395733847835e-05
Newton iteration  2
err =  3.1637348724595416e-10
cost at iteration 84 :  12.074120430820338
Newton iteration  0
err =  0.01703174950953975
Newton iteration  1
err =  1.9917999327007666e-05
Newton iteration  2
err =  3.226325483948901e-10
cost at iteration 85 :  12.053632841810622
Newton iteration  0
err =  0.016801562851512346
Newton iteration  1
err =  1.9981031848223014e-05
Newton iteration  2
err =  3.2915526433688025e-10
cost at iteration 86 :  12.033179966172652
Newton iteration  0
err =  0.016582439596270215
Newton iteration  1
err =  2.004649486746529e-05
Newton iteration  2
err =  3.359501607842966e-10
cost at iteration 87 :  12.012761721868623
Newton iteration  0
err =  0.016373911755880523
Newton iteration  1
err =  2.011439411338072e-05
Newton iteration  2
err =  3.430267176814542e-10
cost at iteration 88 :  11.9923780297781
Newton iteration  0
err =  0.016175524464836742
Newton iteration  1
err =  2.0184739404896924e-05
Newton iteration  2
err =  3.5039366913150407e-10
cost at iteration 89 :  11.972028813530821
Newton iteration  0
err =  0.015986836011493054
Newton iteration  1
err =  2.0257544587512306e-05
Newton iteration  2
err =  3.58061673042036e-10
cost at iteration 90 :  11.951713999349382
Newton iteration  0
err =  0.015807417861197735
Newton iteration  1
err =  2.0332827494493125e-05
Newton iteration  2
err =  3.6604139661135556e-10
cost at iteration 91 :  11.931433515903233
Newton iteration  0
err =  0.015636854667044085
Newton iteration  1
err =  2.0410609908154145e-05
Newton iteration  2
err =  3.7434359164310153e-10
cost at iteration 92 :  11.911187294171327
Newton iteration  0
err =  0.015474744265819422
Newton iteration  1
err =  2.0490917542715592e-05
Newton iteration  2
err =  3.829805289931906e-10
cost at iteration 93 :  11.890975267314026
Newton iteration  0
err =  0.015320697655525444
Newton iteration  1
err =  2.0573780029864628e-05
Newton iteration  2
err =  3.919654526177815e-10
cost at iteration 94 :  11.870797370553184
Newton iteration  0
err =  0.015174338953247718
Newton iteration  1
err =  2.0659230928063215e-05
Newton iteration  2
err =  4.013134672759139e-10
cost at iteration 95 :  11.850653541060415
Newton iteration  0
err =  0.015035305331077577
Newton iteration  1
err =  2.074730772780248e-05
Newton iteration  2
err =  4.110379676218823e-10
cost at iteration 96 :  11.830543717852102
Newton iteration  0
err =  0.01490324692931117
Newton iteration  1
err =  2.0838051892650166e-05
Newton iteration  2
err =  4.21154220623152e-10
cost at iteration 97 :  11.81046784169165
Newton iteration  0
err =  0.014777826746352021
Newton iteration  1
err =  2.0931508883741437e-05
Newton iteration  2
err =  4.316813778380524e-10
cost at iteration 98 :  11.790425854998444
Newton iteration  0
err =  0.014658720504719235
Newton iteration  1
err =  2.1027728228847623e-05
Newton iteration  2
err =  4.4263528697141394e-10
cost at iteration 99 :  11.770417701762034
Newton iteration  0
err =  0.014545616494004859
Newton iteration  1
err =  2.1126763575502378e-05
Newton iteration  2
err =  4.5403539840906605e-10
cost at iteration 100 :  11.750443327463028
Newton iteration  0
err =  0.01443821539114155
Newton iteration  1
err =  2.122867277838644e-05
Newton iteration  2
err =  4.659037570571297e-10
cost at iteration 101 :  11.730502678998207
Newton iteration  0
err =  0.014336230058629842
Newton iteration  1
err =  2.133351799243295e-05
Newton iteration  2
err =  4.782611136332689e-10
cost at iteration 102 :  11.710595704611597
Newton iteration  0
err =  0.01423938532270385
Newton iteration  1
err =  2.144136577921816e-05
Newton iteration  2
err =  4.911307903137802e-10
cost at iteration 103 :  11.69072235382991
Newton iteration  0
err =  0.014147417732197971
Newton iteration  1
err =  2.1552287227326138e-05
Newton iteration  2
err =  5.045399853199991e-10
cost at iteration 104 :  11.670882577401686
Newton iteration  0
err =  0.014060075300214232
Newton iteration  1
err =  2.1666358098721194e-05
Newton iteration  2
err =  5.185109486099702e-10
cost at iteration 105 :  11.651076327241729
Newton iteration  0
err =  0.013977117230139353
Newton iteration  1
err =  2.1783658972485023e-05
Newton iteration  2
err =  5.330748779614502e-10
cost at iteration 106 :  11.631303556378807
Newton iteration  0
err =  0.013898313628141494
Newton iteration  1
err =  2.190427541819263e-05
Newton iteration  2
err =  5.48261561764901e-10
cost at iteration 107 :  11.611564218906677
Newton iteration  0
err =  0.013823445203732573
Newton iteration  1
err =  2.202829818198757e-05
Newton iteration  2
err =  5.64104563503875e-10
cost at iteration 108 :  11.59185826993931
Newton iteration  0
err =  0.013752302960664497
Newton iteration  1
err =  2.2155823386193663e-05
Newton iteration  2
err =  5.806367617766908e-10
cost at iteration 109 :  11.572185665568735
Newton iteration  0
err =  0.01368468787998681
Newton iteration  1
err =  2.2286952753832947e-05
Newton iteration  2
err =  5.9789720255865e-10
cost at iteration 110 :  11.552546362825959
Newton iteration  0
err =  0.01362041059691616
Newton iteration  1
err =  2.2421793839785554e-05
Newton iteration  2
err =  6.159236576046548e-10
cost at iteration 111 :  11.532940319645318
Newton iteration  0
err =  0.01355929107373845
Newton iteration  1
err =  2.256046029099803e-05
Newton iteration  2
err =  6.347601047436843e-10
cost at iteration 112 :  11.513367494830643
Newton iteration  0
err =  0.013501158270124457
Newton iteration  1
err =  2.2703072118057304e-05
Newton iteration  2
err =  6.544522667119729e-10
cost at iteration 113 :  11.49382784802461
Newton iteration  0
err =  0.013445849812561434
Newton iteration  1
err =  2.2849755989266716e-05
Newton iteration  2
err =  6.75048097427665e-10
cost at iteration 114 :  11.474321339680246
Newton iteration  0
err =  0.013393211664539233
Newton iteration  1
err =  2.300064553798688e-05
Newton iteration  2
err =  6.965990842230796e-10
cost at iteration 115 :  11.45484793103473
Newton iteration  0
err =  0.013343097798627279
Newton iteration  1
err =  2.3155881699313018e-05
Newton iteration  2
err =  7.191631784842471e-10
cost at iteration 116 :  11.43540758408521
Newton iteration  0
err =  0.013295369871895071
Newton iteration  1
err =  2.3315613053601007e-05
Newton iteration  2
err =  7.427980683812784e-10
cost at iteration 117 :  11.416000261566968
Newton iteration  0
err =  0.01324989690581398
Newton iteration  1
err =  2.3479996194372114e-05
Newton iteration  2
err =  7.67570677313356e-10
cost at iteration 118 :  11.396625926932856
Newton iteration  0
err =  0.013206554971345267
Newton iteration  1
err =  2.36491961197143e-05
Newton iteration  2
err =  7.935477015028336e-10
cost at iteration 119 :  11.377284544335392
Newton iteration  0
err =  0.013165226880569389
Newton iteration  1
err =  2.3823386628648956e-05
Newton iteration  2
err =  8.20804991172258e-10
cost at iteration 120 :  11.35797607860985
Newton iteration  0
err =  0.013125801885137164
Newton iteration  1
err =  2.4002750746670334e-05
Newton iteration  2
err =  8.494206907645162e-10
cost at iteration 121 :  11.338700495259182
Newton iteration  0
err =  0.013088175382484624
Newton iteration  1
err =  2.4187481157376705e-05
Newton iteration  2
err =  8.794786734006334e-10
cost at iteration 122 :  11.319457760441008
Newton iteration  0
err =  0.013052248630370458
Newton iteration  1
err =  2.4377780654375423e-05
Newton iteration  2
err =  9.110734558766703e-10
cost at iteration 123 :  11.300247840955231
Newton iteration  0
err =  0.01301792846980446
Newton iteration  1
err =  2.457386260109387e-05
Newton iteration  2
err =  9.443018561769567e-10
cost at iteration 124 :  11.281070704233523
Newton iteration  0
err =  0.012985127057254117
Newton iteration  1
err =  2.4775951402067472e-05
Newton iteration  2
err =  9.7926818329193e-10
cost at iteration 125 :  11.261926318330408
Newton iteration  0
err =  0.012953761605895824
Newton iteration  1
err =  2.498428297817992e-05
Newton iteration  2
err =  1.0160863631884724e-09
cost at iteration 126 :  11.24281465191515
Newton iteration  0
err =  0.012923754136357716
Newton iteration  1
err =  2.519910524733234e-05
Newton iteration  2
err =  1.0548784742779523e-09
cost at iteration 127 :  11.223735674265345
Newton iteration  0
err =  0.01289503123703417
Newton iteration  1
err =  2.5420678605761696e-05
Newton iteration  2
err =  1.0957739518233584e-09
cost at iteration 128 :  11.20468935526117
Newton iteration  0
err =  0.012867523833871042
Newton iteration  1
err =  2.564927639418488e-05
Newton iteration  2
err =  1.1389127165181045e-09
cost at iteration 129 :  11.185675665381524
Newton iteration  0
err =  0.012841166969829003
Newton iteration  1
err =  2.5885185364067842e-05
Newton iteration  2
err =  1.1844455041994396e-09
cost at iteration 130 :  11.1666945757007
Newton iteration  0
err =  0.012815899593699812
Newton iteration  1
err =  2.612870611265229e-05
Newton iteration  2
err =  1.232533474787202e-09
cost at iteration 131 :  11.147746057886511
Newton iteration  0
err =  0.012791664358429943
Newton iteration  1
err =  2.6380153505690426e-05
Newton iteration  2
err =  1.2833479819657505e-09
cost at iteration 132 :  11.12883008419928
Newton iteration  0
err =  0.012768407428599802
Newton iteration  1
err =  2.6639857055408187e-05
Newton iteration  2
err =  1.3370762977409443e-09
cost at iteration 133 :  11.10994662749188
Newton iteration  0
err =  0.012746078296952785
Newton iteration  1
err =  2.690816125670828e-05
Newton iteration  2
err =  1.393915233529478e-09
cost at iteration 134 :  11.091095661211272
Newton iteration  0
err =  0.012724629609825467
Newton iteration  1
err =  2.7185425871996256e-05
Newton iteration  2
err =  1.4540795181166342e-09
cost at iteration 135 :  11.072277159400391
Newton iteration  0
err =  0.012704017001136833
Newton iteration  1
err =  2.747202615075324e-05
Newton iteration  2
err =  1.5177930923282437e-09
cost at iteration 136 :  11.053491096701482
Newton iteration  0
err =  0.012684198934655378
Newton iteration  1
err =  2.7768352966197956e-05
Newton iteration  2
err =  1.5853022893425856e-09
cost at iteration 137 :  11.034737448360112
Newton iteration  0
err =  0.012665136554533088
Newton iteration  1
err =  2.8074812864442036e-05
Newton iteration  2
err =  1.6568643119494645e-09
cost at iteration 138 :  11.016016190230676
Newton iteration  0
err =  0.012646793543384932
Newton iteration  1
err =  2.8391828007110677e-05
Newton iteration  2
err =  1.732756715136318e-09
cost at iteration 139 :  10.997327298782277
Newton iteration  0
err =  0.0126291359880176
Newton iteration  1
err =  2.871983599469449e-05
Newton iteration  2
err =  1.8132740508754525e-09
cost at iteration 140 :  10.978670751106096
Newton iteration  0
err =  0.012612132252244815
Newton iteration  1
err =  2.9059289556744532e-05
Newton iteration  2
err =  1.898728731499999e-09
cost at iteration 141 :  10.960046524923301
Newton iteration  0
err =  0.012595752856544335
Newton iteration  1
err =  2.9410656097024283e-05
Newton iteration  2
err =  1.989451688034235e-09
cost at iteration 142 :  10.94145459859431
Newton iteration  0
err =  0.012579970364208846
Newton iteration  1
err =  2.9774417074828568e-05
Newton iteration  2
err =  2.0857948895311195e-09
cost at iteration 143 :  10.922894951128796
Newton iteration  0
err =  0.012564759273698973
Newton iteration  1
err =  3.015106722799608e-05
Newton iteration  2
err =  2.1881262140167715e-09
cost at iteration 144 :  10.904367562196974
Newton iteration  0
err =  0.012550095916704536
Newton iteration  1
err =  3.0541113628607755e-05
Newton iteration  2
err =  2.2968364254121554e-09
cost at iteration 145 :  10.885872412140944
Newton iteration  0
err =  0.01253595836177467
Newton iteration  1
err =  3.094507457994192e-05
Newton iteration  2
err =  2.412329711945419e-09
cost at iteration 146 :  10.867409481988647
Newton iteration  0
err =  0.012522326322899202
Newton iteration  1
err =  3.13634783654517e-05
Newton iteration  2
err =  2.535032036065279e-09
cost at iteration 147 :  10.848978753466945
Newton iteration  0
err =  0.012509181072888996
Newton iteration  1
err =  3.179686190408378e-05
Newton iteration  2
err =  2.665383852604303e-09
cost at iteration 148 :  10.830580209017246
Newton iteration  0
err =  0.012496505361047754
Newton iteration  1
err =  3.224576936246809e-05
Newton iteration  2
err =  2.803838911677908e-09
cost at iteration 149 :  10.812213831810764
Newton iteration  0
err =  0.012484283334820089
Newton iteration  1
err =  3.271075080353895e-05
Newton iteration  2
err =  2.9508666115106892e-09
cost at iteration 150 :  10.79387960576576
Newton iteration  0
err =  0.01247250046494682
Newton iteration  1
err =  3.3192361037198554e-05
Newton iteration  2
err =  3.1069454092879777e-09
cost at iteration 151 :  10.7755775155651
Newton iteration  0
err =  0.012461143473996093
Newton iteration  1
err =  3.369115883843271e-05
Newton iteration  2
err =  3.2725598046699425e-09
cost at iteration 152 :  10.757307546674973
Newton iteration  0
err =  0.012450200267533315
Newton iteration  1
err =  3.420770680369968e-05
Newton iteration  2
err =  3.4482065726890676e-09
cost at iteration 153 :  10.73906968536392
Newton iteration  0
err =  0.0124396598680036
Newton iteration  1
err =  3.47425722001425e-05
Newton iteration  2
err =  3.634385129403419e-09
cost at iteration 154 :  10.720863918723941
Newton iteration  0
err =  0.012429512350786064
Newton iteration  1
err =  3.529632925498369e-05
Newton iteration  2
err =  3.831602352332445e-09
cost at iteration 155 :  10.702690234690634
Newton iteration  0
err =  0.012419748782284053
Newton iteration  1
err =  3.5869563490169264e-05
Newton iteration  2
err =  4.04037551724632e-09
cost at iteration 156 :  10.68454862206549
Newton iteration  0
err =  0.012410361159985889
Newton iteration  1
err =  3.646287880226096e-05
Newton iteration  2
err =  4.261241050538261e-09
cost at iteration 157 :  10.66643907053829
Newton iteration  0
err =  0.012401342354375485
Newton iteration  1
err =  3.707690820468913e-05
Newton iteration  2
err =  4.494764212998257e-09
cost at iteration 158 :  10.648361570709627
Newton iteration  0
err =  0.012392686053052757
Newton iteration  1
err =  3.771232923082586e-05
Newton iteration  2
err =  4.741566314469282e-09
cost at iteration 159 :  10.630316114114846
Newton iteration  0
err =  0.012384386707272339
Newton iteration  1
err =  3.8369885186867116e-05
Newton iteration  2
err =  5.0023479379468866e-09
cost at iteration 160 :  10.612302693247567
Newton iteration  0
err =  0.012376439481937729
Newton iteration  1
err =  3.905041344994022e-05
Newton iteration  2
err =  5.2779387052899975e-09
cost at iteration 161 :  10.594321301583989
Newton iteration  0
err =  0.012368840209780777
Newton iteration  1
err =  3.975488198813815e-05
Newton iteration  2
err =  5.569352490039284e-09
cost at iteration 162 :  10.576371933607524
Newton iteration  0
err =  0.01236158535155964
Newton iteration  1
err =  4.048443505122463e-05
Newton iteration  2
err =  5.877858953409508e-09
cost at iteration 163 :  10.55845458483335
Newton iteration  0
err =  0.012354671964298906
Newton iteration  1
err =  4.124044847075189e-05
Newton iteration  2
err =  6.205074367692843e-09
cost at iteration 164 :  10.54056925183342
Newton iteration  0
err =  0.012348097680111816
Newton iteration  1
err =  4.2024594228211534e-05
Newton iteration  2
err =  6.553066371884739e-09
cost at iteration 165 :  10.522715932262177
Newton iteration  0
err =  0.01234186069892655
Newton iteration  1
err =  4.283891270554718e-05
Newton iteration  2
err =  6.9244586919814264e-09
cost at iteration 166 :  10.504894624882294
Newton iteration  0
err =  0.012335959798989264
Newton iteration  1
err =  4.3685889375221465e-05
Newton iteration  2
err =  7.322552687718486e-09
cost at iteration 167 :  10.48710532959084
Newton iteration  0
err =  0.012330394369104367
Newton iteration  1
err =  4.456853075168658e-05
Newton iteration  2
err =  7.751426945462715e-09
cost at iteration 168 :  10.469348047448188
Newton iteration  0
err =  0.012325164467236461
Newton iteration  1
err =  4.549043224489716e-05
Newton iteration  2
err =  8.216033633436336e-09
cost at iteration 169 :  10.451622780706623
Newton iteration  0
err =  0.012320270909423678
Newton iteration  1
err =  4.64558287751062e-05
Newton iteration  2
err =  8.72225710114321e-09
cost at iteration 170 :  10.433929532842695
Newton iteration  0
err =  0.012315715392350192
Newton iteration  1
err =  4.746961799161584e-05
Newton iteration  2
err =  9.276960180847907e-09
cost at iteration 171 :  10.416268308592462
Newton iteration  0
err =  0.012311500651797438
Newton iteration  1
err =  4.853734658012786e-05
Newton iteration  2
err =  9.887989606364224e-09
cost at iteration 172 :  10.398639113991496
Newton iteration  0
err =  0.012307630656274677
Newton iteration  1
err =  4.9665153013417095e-05
Newton iteration  2
err =  1.0564187102901804e-08
Newton iteration  3
err =  1.8233564040988024e-14
cost at iteration 173 :  10.381041956420185
Newton iteration  0
err =  0.012304110833451495
Newton iteration  1
err =  5.085966562464336e-05
Newton iteration  2
err =  1.1315390989351814e-08
Newton iteration  3
err =  1.7787299707904648e-14
cost at iteration 174 :  10.363476844657052
Newton iteration  0
err =  0.01230094832267872
Newton iteration  1
err =  5.212786276184161e-05
Newton iteration  2
err =  1.2152475052473384e-08
Newton iteration  3
err =  1.7349121556607614e-14
cost at iteration 175 :  10.345943788940067
Newton iteration  0
err =  0.012298152244176751
Newton iteration  1
err =  5.3476910883903163e-05
Newton iteration  2
err =  1.3087412704283389e-08
Newton iteration  3
err =  1.871455396751302e-14
cost at iteration 176 :  10.328442801039486
Newton iteration  0
err =  0.012295733972304414
Newton iteration  1
err =  5.4914004689581795e-05
Newton iteration  2
err =  1.4133413158873195e-08
Newton iteration  3
err =  1.8808958329644537e-14
cost at iteration 177 :  10.310973894341576
Newton iteration  0
err =  0.01229370739825941
Newton iteration  1
err =  5.644623835126511e-05
Newton iteration  2
err =  1.5305099215434312e-08
Newton iteration  3
err =  1.9326667953338592e-14
cost at iteration 178 :  10.293537083947292
Newton iteration  0
err =  0.012292089167024359
Newton iteration  1
err =  5.808053629417801e-05
Newton iteration  2
err =  1.661875188392611e-08
Newton iteration  3
err =  1.8109226494483065e-14
cost at iteration 179 :  10.276132386785072
Newton iteration  0
err =  0.012290898874954062
Newton iteration  1
err =  5.9823664946027574e-05
Newton iteration  2
err =  1.8092613091068748e-08
Newton iteration  3
err =  1.7564372682958154e-14
cost at iteration 180 :  10.258759821739352
Newton iteration  0
err =  0.012290159216890376
Newton iteration  1
err =  6.16823343295594e-05
Newton iteration  2
err =  1.9747261659996388e-08
Newton iteration  3
err =  1.8765461197360116e-14
cost at iteration 181 :  10.241419409796512
Newton iteration  0
err =  0.012289896076918064
Newton iteration  1
err =  6.366338323003658e-05
Newton iteration  2
err =  2.16060791505228e-08
Newton iteration  3
err =  1.7787466749043256e-14
cost at iteration 182 :  10.224111174207465
Newton iteration  0
err =  0.012290138561118276
Newton iteration  1
err =  6.577402774877794e-05
Newton iteration  2
err =  2.3695820893734014e-08
Newton iteration  3
err =  1.801377840474905e-14
cost at iteration 183 :  10.206835140667062
Newton iteration  0
err =  0.012290918976429973
Newton iteration  1
err =  6.802214420423649e-05
Newton iteration  2
err =  2.6047313122456472e-08
Newton iteration  3
err =  1.9448335228159856e-14
cost at iteration 184 :  10.189591337511139
Newton iteration  0
err =  0.012292272762426974
Newton iteration  1
err =  7.0416555552149e-05
Newton iteration  2
err =  2.8696285519762736e-08
Newton iteration  3
err =  1.884096207872271e-14
cost at iteration 185 :  10.172379795929428
Newton iteration  0
err =  0.01229423838576332
Newton iteration  1
err =  7.296729618248695e-05
Newton iteration  2
err =  3.168433802872885e-08
Newton iteration  3
err =  2.078065411271825e-14
cost at iteration 186 :  10.155200550194584
Newton iteration  0
err =  0.012296857206278145
Newton iteration  1
err =  7.56858413437247e-05
Newton iteration  2
err =  3.50600525264954e-08
Newton iteration  3
err =  2.0564839526370023e-14
cost at iteration 187 :  10.138053637906946
Newton iteration  0
err =  0.01230017332275999
Newton iteration  1
err =  7.858530212621059e-05
Newton iteration  2
err =  3.888028032501466e-08
Newton iteration  3
err =  2.2868968974132482e-14
cost at iteration 188 :  10.120939100253366
Newton iteration  0
err =  0.012304233403718912
Newton iteration  1
err =  8.168060273910043e-05
Newton iteration  2
err =  4.321172135047506e-08
Newton iteration  3
err =  2.5096560457613616e-14
cost at iteration 189 :  10.103856982279527
Newton iteration  0
err =  0.01230908650697172
Newton iteration  1
err =  8.498867215398047e-05
Newton iteration  2
err =  4.8133011040960096e-08
Newton iteration  3
err =  2.894508299023985e-14
cost at iteration 190 :  10.086807333176145
Newton iteration  0
err =  0.01231478389176927
Newton iteration  1
err =  8.852869677817834e-05
Newton iteration  2
err =  5.373772185714223e-08
Newton iteration  3
err =  3.384465764475792e-14
cost at iteration 191 :  10.069790206576535
Newton iteration  0
err =  0.0123213788309675
Newton iteration  1
err =  9.232249433997831e-05
Newton iteration  2
err =  6.013895852418721e-08
Newton iteration  3
err =  4.087117127122929e-14
cost at iteration 192 :  10.052805660865943
Newton iteration  0
err =  0.012328926440053645
Newton iteration  1
err =  9.639508034211186e-05
Newton iteration  2
err =  6.747654059348854e-08
Newton iteration  3
err =  5.0028235460533076e-14
cost at iteration 193 :  10.035853759502295
Newton iteration  0
err =  0.012337483556003531
Newton iteration  1
err =  0.0001007755028593423
Newton iteration  2
err =  7.592819852547094e-08
Newton iteration  3
err =  6.165859050605587e-14
cost at iteration 194 :  10.018934571347906
Newton iteration  0
err =  0.012347108721927739
Newton iteration  1
err =  0.00010549800964131747
Newton iteration  2
err =  8.572656006446195e-08
Newton iteration  3
err =  7.766143612890255e-14
cost at iteration 195 :  10.002048171014076
Newton iteration  0
err =  0.012357862361309683
Newton iteration  1
err =  0.0001106035672606993
Newton iteration  2
err =  9.718389217316569e-08
Newton iteration  3
err =  9.941112682673601e-14
cost at iteration 196 :  9.985194639220616
Newton iteration  0
err =  0.012369807247722811
Newton iteration  1
err =  0.00011614165387388719
Newton iteration  2
err =  1.1072622013189192e-07
Newton iteration  3
err =  1.291765036593663e-13
cost at iteration 197 :  9.968374063176789
Newton iteration  0
err =  0.012383009373436452
Newton iteration  1
err =  0.00012217207389195398
Newton iteration  2
err =  1.2693727999214811e-07
Newton iteration  3
err =  1.7158506037343033e-13
cost at iteration 198 :  9.951586536990716
Newton iteration  0
err =  0.012397539255241416
Newton iteration  1
err =  0.00012876628622800467
Newton iteration  2
err =  1.4661056483943413e-07
Newton iteration  3
err =  2.3486002309556344e-13
cost at iteration 199 :  9.934832162119742
Newton iteration  0
err =  0.012413473528612424
Newton iteration  1
err =  0.00013600741548291289
Newton iteration  2
err =  1.7080467304352815e-07
Newton iteration  3
err =  3.318844884080939e-13
cost at iteration 200 :  9.918111047873042
Newton iteration  0
err =  0.012430896288924097
Newton iteration  1
err =  0.00014398776886633165
Newton iteration  2
err =  2.008938217699363e-07
Newton iteration  3
err =  4.834922260035211e-13
cost at iteration 201 :  9.901423311976957
Newton iteration  0
err =  0.012449898936251718
Newton iteration  1
err =  0.0001528023314832948
Newton iteration  2
err =  2.3860148924241467e-07
Newton iteration  3
err =  7.242987490647969e-13
cost at iteration 202 :  9.88476908119633
Newton iteration  0
err =  0.01247057614923254
Newton iteration  1
err =  0.00016253626045895244
Newton iteration  2
err =  2.8599747434497306e-07
Newton iteration  3
err =  1.1084882808312597e-12
cost at iteration 203 :  9.868148491975601
Newton iteration  0
err =  0.012493013915133206
Newton iteration  1
err =  0.00017324352100803514
Newton iteration  2
err =  3.4541831917974555e-07
Newton iteration  3
err =  1.7247213047934312e-12
cost at iteration 204 :  9.85156169099608
Newton iteration  0
err =  0.012517263083462142
Newton iteration  1
err =  0.00018491197370533902
Newton iteration  2
err =  4.192219678381694e-07
Newton iteration  3
err =  2.7070359429132887e-12
cost at iteration 205 :  9.83500883543544
Newton iteration  0
err =  0.0125432883813894
Newton iteration  1
err =  0.00019740701464829262
Newton iteration  2
err =  5.091865962984497e-07
Newton iteration  3
err =  4.251363827784453e-12
cost at iteration 206 :  9.81849009252109
Newton iteration  0
err =  0.012570877816956915
Newton iteration  1
err =  0.00021038189387592853
Newton iteration  2
err =  6.151900786917023e-07
Newton iteration  3
err =  6.609292927535375e-12
cost at iteration 207 :  9.802005637667902
Newton iteration  0
err =  0.012599490802467373
Newton iteration  1
err =  0.00022314209699310388
Newton iteration  2
err =  7.326004776070764e-07
Newton iteration  3
err =  9.988610152632147e-12
cost at iteration 208 :  9.785555650020397
Newton iteration  0
err =  0.01262801715316255
Newton iteration  1
err =  0.00023446765222186924
Newton iteration  2
err =  8.479037836207625e-07
Newton iteration  3
err =  1.4218161692357797e-11
cost at iteration 209 :  9.769140303575066
Newton iteration  0
err =  0.012654422777052572
Newton iteration  1
err =  0.00024247016156277073
Newton iteration  2
err =  9.336824296693441e-07
Newton iteration  3
err =  1.8102876993505744e-11
cost at iteration 210 :  9.752759751352622
Newton iteration  0
err =  0.012675300323488287
Newton iteration  1
err =  0.000244770602460703
Newton iteration  2
err =  9.492701950481514e-07
Newton iteration  3
err =  1.9129926515377813e-11
cost at iteration 211 :  9.736414099903048
Newton iteration  0
err =  0.012685486689362702
Newton iteration  1
err =  0.0002396685808087788
Newton iteration  2
err =  8.605358063796654e-07
Newton iteration  3
err =  1.5352502899024855e-11
cost at iteration 212 :  9.720103373262262
Newton iteration  0
err =  0.012678212691960454
Newton iteration  1
err =  0.00022900847970583473
Newton iteration  2
err =  6.840324732515209e-07
Newton iteration  3
err =  8.805191704451343e-12
cost at iteration 213 :  9.703827471531778
Newton iteration  0
err =  0.012646541206614497
Newton iteration  1
err =  0.0002208957185917427
Newton iteration  2
err =  5.22966689346756e-07
Newton iteration  3
err =  4.1416630326860395e-12
cost at iteration 214 :  9.687586139223368
Newton iteration  0
err =  0.012586392844751335
Newton iteration  1
err =  0.00022450335903146442
Newton iteration  2
err =  5.211695561220227e-07
Newton iteration  3
err =  3.8258944850706094e-12
cost at iteration 215 :  9.671378963930431
Newton iteration  0
err =  0.012499633096018899
Newton iteration  1
err =  0.0002372786082359165
Newton iteration  2
err =  6.614704289140628e-07
Newton iteration  3
err =  7.14451032496481e-12
cost at iteration 216 :  9.655205414130458
Newton iteration  0
err =  0.012394274759087597
Newton iteration  1
err =  0.00024707393778918864
Newton iteration  2
err =  7.887869700017052e-07
Newton iteration  3
err =  1.08018893333042e-11
cost at iteration 217 :  9.639064900436525
Newton iteration  0
err =  0.012280998957710105
Newton iteration  1
err =  0.00024622226405694047
Newton iteration  2
err =  8.147828082753429e-07
Newton iteration  3
err =  1.1504768130178703e-11
cost at iteration 218 :  9.62295683269248
Newton iteration  0
err =  0.012168955736803157
Newton iteration  1
err =  0.00023513476937008805
Newton iteration  2
err =  7.456842229774209e-07
Newton iteration  3
err =  9.34197911162977e-12
cost at iteration 219 :  9.606880656389658
Newton iteration  0
err =  0.012063789570376734
Newton iteration  1
err =  0.00021784170102066448
Newton iteration  2
err =  6.284615085426893e-07
Newton iteration  3
err =  6.371559931310892e-12
cost at iteration 220 :  9.590835868498706
Newton iteration  0
err =  0.01196793615792913
Newton iteration  1
err =  0.00019822397766633868
Newton iteration  2
err =  5.039291793968795e-07
Newton iteration  3
err =  3.926784384183423e-12
cost at iteration 221 :  9.574822020125978
Newton iteration  0
err =  0.01188177468161272
Newton iteration  1
err =  0.0001787905223304271
Newton iteration  2
err =  3.934708505207827e-07
Newton iteration  3
err =  2.294102495784768e-12
cost at iteration 222 :  9.558838712805453
Newton iteration  0
err =  0.011804654631116875
Newton iteration  1
err =  0.00016080852190692008
Newton iteration  2
err =  3.0377791865066424e-07
Newton iteration  3
err =  1.3102969066548512e-12
cost at iteration 223 :  9.542885592522522
Newton iteration  0
err =  0.011735537866697104
Newton iteration  1
err =  0.0001447521705739963
Newton iteration  2
err =  2.341364977884601e-07
Newton iteration  3
err =  7.454128240984134e-13
cost at iteration 224 :  9.526962343458358
Newton iteration  0
err =  0.011673328150561187
Newton iteration  1
err =  0.00013067310191967093
Newton iteration  2
err =  1.8119767720502107e-07
Newton iteration  3
err =  4.274975110605572e-13
cost at iteration 225 :  9.51106868226787
Newton iteration  0
err =  0.011617014009844835
Newton iteration  1
err =  0.00011842963885865633
Newton iteration  2
err =  1.4126355874317424e-07
Newton iteration  3
err =  2.4905222740444714e-13
cost at iteration 226 :  9.495204353139398
Newton iteration  0
err =  0.011565716250459919
Newton iteration  1
err =  0.00010781010052778667
Newton iteration  2
err =  1.1113371406889459e-07
Newton iteration  3
err =  1.4827984752100391e-13
cost at iteration 227 :  9.479369123642961
Newton iteration  0
err =  0.01151869208482935
Newton iteration  1
err =  9.859280053177322e-05
Newton iteration  2
err =  8.829261205829106e-08
Newton iteration  3
err =  9.040825947731074e-14
cost at iteration 228 :  9.463562781282617
Newton iteration  0
err =  0.011475322082215923
Newton iteration  1
err =  9.05721292484673e-05
Newton iteration  2
err =  7.084949711119817e-08
Newton iteration  3
err =  5.723896100731357e-14
cost at iteration 229 :  9.447785130644588
Newton iteration  0
err =  0.011435092112948192
Newton iteration  1
err =  8.356771945999547e-05
Newton iteration  2
err =  5.741230685377768e-08
Newton iteration  3
err =  3.877100067040091e-14
cost at iteration 230 :  9.432035991037212
Newton iteration  0
err =  0.011397575470893687
Newton iteration  1
err =  7.742572790745611e-05
Newton iteration  2
err =  4.6964240445235616e-08
Newton iteration  3
err =  2.772227345326458e-14
cost at iteration 231 :  9.41631519453251
Newton iteration  0
err =  0.011362417083487081
Newton iteration  1
err =  7.201678665348485e-05
Newton iteration  2
err =  3.876312563224858e-08
Newton iteration  3
err =  2.2218246478313102e-14
cost at iteration 232 :  9.40062258433328
Newton iteration  0
err =  0.011329320268357642
Newton iteration  1
err =  6.723281973122311e-05
Newton iteration  2
err =  3.226534187028741e-08
Newton iteration  3
err =  1.9198398906464284e-14
cost at iteration 233 :  9.38495801340529
Newton iteration  0
err =  0.01129803590576554
Newton iteration  1
err =  6.298372621363269e-05
Newton iteration  2
err =  2.707030943748325e-08
Newton iteration  3
err =  1.9041089393383047e-14
cost at iteration 234 :  9.36932134332474
Newton iteration  0
err =  0.01126835369522867
Newton iteration  1
err =  5.919434424469781e-05
Newton iteration  2
err =  2.2880660244994022e-08
Newton iteration  3
err =  1.7097884672059588e-14
cost at iteration 235 :  9.35371244330148
Newton iteration  0
err =  0.01124009513181097
Newton iteration  1
err =  5.580183218689942e-05
Newton iteration  2
err =  1.9473834896623463e-08
Newton iteration  3
err =  1.677509532865484e-14
cost at iteration 236 :  9.338131189347598
Newton iteration  0
err =  0.011213107869280954
Newton iteration  1
err =  5.275347774490113e-05
Newton iteration  2
err =  1.6681841626400464e-08
Newton iteration  3
err =  1.6964672898322622e-14
cost at iteration 237 :  9.322577463564683
Newton iteration  0
err =  0.011187261189272615
Newton iteration  1
err =  5.0004894921011636e-05
Newton iteration  2
err =  1.4376785249343987e-08
Newton iteration  3
err =  1.6444674769673342e-14
cost at iteration 238 :  9.307051153529985
Newton iteration  0
err =  0.011162442348189949
Newton iteration  1
err =  4.7518552680031076e-05
Newton iteration  2
err =  1.2460468455832e-08
Newton iteration  3
err =  1.6081874153370927e-14
cost at iteration 239 :  9.291552151765647
Newton iteration  0
err =  0.011138553619880219
Newton iteration  1
err =  4.52625788364844e-05
Newton iteration  2
err =  1.0856870145464966e-08
Newton iteration  3
err =  1.6909940691669037e-14
cost at iteration 240 :  9.27608035527692
Newton iteration  0
err =  0.011115509890539213
Newton iteration  1
err =  4.320978850616844e-05
Newton iteration  2
err =  9.506662684896346e-09
cost at iteration 241 :  9.260635665148042
Newton iteration  0
err =  0.01109323669328479
Newton iteration  1
err =  4.133689418979085e-05
Newton iteration  2
err =  8.363190255496531e-09
cost at iteration 242 :  9.245217986188868
Newton iteration  0
err =  0.011071668593876668
Newton iteration  1
err =  3.9623862108218774e-05
Newton iteration  2
err =  7.389490675793974e-09
cost at iteration 243 :  9.229827226622453
Newton iteration  0
err =  0.011050747858202132
Newton iteration  1
err =  3.8053386170797945e-05
Newton iteration  2
err =  6.556072433715143e-09
cost at iteration 244 :  9.214463297808825
Newton iteration  0
err =  0.011030423346911684
Newton iteration  1
err =  3.6610456584234745e-05
Newton iteration  2
err =  5.839248769735352e-09
cost at iteration 245 :  9.199126114000364
Newton iteration  0
err =  0.011010649593856916
Newton iteration  1
err =  3.528200476058788e-05
Newton iteration  2
err =  5.219869923687483e-09
cost at iteration 246 :  9.183815592122667
Newton iteration  0
err =  0.010991386034311827
Newton iteration  1
err =  3.405660986798309e-05
Newton iteration  2
err =  4.682362556963818e-09
cost at iteration 247 :  9.168531651579954
Newton iteration  0
err =  0.010972596355459385
Newton iteration  1
err =  3.292425537989784e-05
Newton iteration  2
err =  4.21398565053313e-09
cost at iteration 248 :  9.153274214078976
Newton iteration  0
err =  0.01095424794744804
Newton iteration  1
err =  3.1876126267885305e-05
Newton iteration  2
err =  3.804259953751312e-09
cost at iteration 249 :  9.138043203472176
Newton iteration  0
err =  0.010936311437394968
Newton iteration  1
err =  3.090443940851102e-05
Newton iteration  2
err =  3.4445183710506476e-09
cost at iteration 250 :  9.122838545614764
Newton iteration  0
err =  0.010918760292181086
Newton iteration  1
err =  3.0002301188172677e-05
Newton iteration  2
err =  3.1275580074778088e-09
cost at iteration 251 :  9.107660168236208
Newton iteration  0
err =  0.010901570478566271
Newton iteration  1
err =  2.9163587494163017e-05
Newton iteration  2
err =  2.847365996770924e-09
cost at iteration 252 :  9.092508000823408
Newton iteration  0
err =  0.01088472017117744
Newton iteration  1
err =  2.8382842196245974e-05
Newton iteration  2
err =  2.5988969336474037e-09
cost at iteration 253 :  9.077381974514356
Newton iteration  0
err =  0.010868189500704546
Newton iteration  1
err =  2.7655190947172904e-05
Newton iteration  2
err =  2.3778995973321764e-09
cost at iteration 254 :  9.062282022001302
Newton iteration  0
err =  0.010851960336147954
Newton iteration  1
err =  2.6976267739107517e-05
Newton iteration  2
err =  2.1807802599890927e-09
cost at iteration 255 :  9.047208077442363
Newton iteration  0
err =  0.010836016095607323
Newton iteration  1
err =  2.6342152106898755e-05
Newton iteration  2
err =  2.0044832744903863e-09
cost at iteration 256 :  9.03216007638003
Newton iteration  0
err =  0.010820341581602564
Newton iteration  1
err =  2.5749315261149977e-05
Newton iteration  2
err =  1.8464049019669568e-09
cost at iteration 257 :  9.017137955667316
Newton iteration  0
err =  0.010804922837127553
Newton iteration  1
err =  2.5194573725586567e-05
Newton iteration  2
err =  1.7043189951899303e-09
cost at iteration 258 :  9.002141653399084
Newton iteration  0
err =  0.010789747019568079
Newton iteration  1
err =  2.4675049307022844e-05
Newton iteration  2
err =  1.5763115263483126e-09
cost at iteration 259 :  8.987171108848893
Newton iteration  0
err =  0.010774802289776652
Newton iteration  1
err =  2.4188134439868054e-05
Newton iteration  2
err =  1.4607360117503682e-09
cost at iteration 260 :  8.972226262411315
Newton iteration  0
err =  0.010760077714450814
Newton iteration  1
err =  2.3731462078858236e-05
Newton iteration  2
err =  1.3561691815687223e-09
cost at iteration 261 :  8.957307055547618
Newton iteration  0
err =  0.010745563179758718
Newton iteration  1
err =  2.3302879477351736e-05
Newton iteration  2
err =  1.2613768675618173e-09
cost at iteration 262 :  8.94241343073645
Newton iteration  0
err =  0.010731249314932488
Newton iteration  1
err =  2.2900425291110337e-05
Newton iteration  2
err =  1.1752881634939582e-09
cost at iteration 263 :  8.927545331427662
Newton iteration  0
err =  0.010717127424264446
Newton iteration  1
err =  2.2522309516230512e-05
Newton iteration  2
err =  1.0969672547159895e-09
cost at iteration 264 :  8.912702701999114
Newton iteration  0
err =  0.010703189426732842
Newton iteration  1
err =  2.2166895884663442e-05
Newton iteration  2
err =  1.0255986352976929e-09
cost at iteration 265 :  8.897885487717275
Newton iteration  0
err =  0.010689427802012782
Newton iteration  1
err =  2.1832686362993954e-05
Newton iteration  2
err =  9.604673230285925e-10
cost at iteration 266 :  8.883093634699643
Newton iteration  0
err =  0.010675835542283288
Newton iteration  1
err =  2.151830747443586e-05
Newton iteration  2
err =  9.009460237232743e-10
cost at iteration 267 :  8.868327089880596
Newton iteration  0
err =  0.010662406108988382
Newton iteration  1
err =  2.1222498206759814e-05
Newton iteration  2
err =  8.464810795787929e-10
cost at iteration 268 :  8.853585800978708
Newton iteration  0
err =  0.01064913339404595
Newton iteration  1
err =  2.0944099282605742e-05
Newton iteration  2
err =  7.96586408505783e-10
cost at iteration 269 :  8.838869716466716
Newton iteration  0
err =  0.010636011684949834
Newton iteration  1
err =  2.0682043641563114e-05
Newton iteration  2
err =  7.508321746296746e-10
cost at iteration 270 :  8.824178785543264
Newton iteration  0
err =  0.010623035633293374
Newton iteration  1
err =  2.0435347951697033e-05
Newton iteration  2
err =  7.088359571098894e-10
cost at iteration 271 :  8.809512958106332
Newton iteration  0
err =  0.010610200226403915
Newton iteration  1
err =  2.0203105050455344e-05
Newton iteration  2
err =  6.702615984309886e-10
cost at iteration 272 :  8.794872184728579
Newton iteration  0
err =  0.0105975007615993
Newton iteration  1
err =  1.9984477172279958e-05
Newton iteration  2
err =  6.348078483870471e-10
cost at iteration 273 :  8.78025641663373
Newton iteration  0
err =  0.01058493282296307
Newton iteration  1
err =  1.9778689900939212e-05
Newton iteration  2
err =  6.022078684365181e-10
cost at iteration 274 :  8.765665605675103
Newton iteration  0
err =  0.010572492260166036
Newton iteration  1
err =  1.958502672355081e-05
Newton iteration  2
err =  5.722250135421281e-10
cost at iteration 275 :  8.751099704314528
Newton iteration  0
err =  0.01056017516934909
Newton iteration  1
err =  1.9402824156208822e-05
Newton iteration  2
err =  5.446455358956044e-10
cost at iteration 276 :  8.736558665603429
Newton iteration  0
err =  0.010547977875574177
Newton iteration  1
err =  1.923146733222954e-05
Newton iteration  2
err =  5.192814229874704e-10
cost at iteration 277 :  8.722042443164213
Newton iteration  0
err =  0.010535896916934864
Newton iteration  1
err =  1.907038604805638e-05
Newton iteration  2
err =  4.959608447190959e-10
cost at iteration 278 :  8.707550991173367
Newton iteration  0
err =  0.010523929030007702
Newton iteration  1
err =  1.891905116748715e-05
Newton iteration  2
err =  4.745318490291801e-10
cost at iteration 279 :  8.693084264345016
Newton iteration  0
err =  0.010512071136594888
Newton iteration  1
err =  1.8776971385383153e-05
Newton iteration  2
err =  4.5486010156903315e-10
cost at iteration 280 :  8.678642217915655
Newton iteration  0
err =  0.010500320331534438
Newton iteration  1
err =  1.8643690290093208e-05
Newton iteration  2
err =  4.368194105503116e-10
cost at iteration 281 :  8.66422480762989
Newton iteration  0
err =  0.010488673871558945
Newton iteration  1
err =  1.8518783697116382e-05
Newton iteration  2
err =  4.203002790227481e-10
cost at iteration 282 :  8.649831989726483
Newton iteration  0
err =  0.010477129165134757
Newton iteration  1
err =  1.8401857226762575e-05
Newton iteration  2
err =  4.0520430764477503e-10
cost at iteration 283 :  8.635463720925639
Newton iteration  0
err =  0.010465683763020274
Newton iteration  1
err =  1.8292544107491393e-05
Newton iteration  2
err =  3.9143902275094776e-10
cost at iteration 284 :  8.621119958416713
Newton iteration  0
err =  0.010454335349679222
Newton iteration  1
err =  1.819050316594789e-05
Newton iteration  2
err =  3.78923333266539e-10
cost at iteration 285 :  8.606800659846714
Newton iteration  0
err =  0.010443081735345167
Newton iteration  1
err =  1.8095417004308304e-05
Newton iteration  2
err =  3.6758224293467255e-10
cost at iteration 286 :  8.59250578330947
Newton iteration  0
err =  0.010431920848720777
Newton iteration  1
err =  1.8006990327490742e-05
Newton iteration  2
err =  3.573469164011389e-10
cost at iteration 287 :  8.578235287335051
Newton iteration  0
err =  0.010420850730261727
Newton iteration  1
err =  1.792494842911105e-05
Newton iteration  2
err =  3.481577536261013e-10
cost at iteration 288 :  8.563989130880133
Newton iteration  0
err =  0.010409869525946783
Newton iteration  1
err =  1.7849035787593058e-05
Newton iteration  2
err =  3.399536765081022e-10
cost at iteration 289 :  8.549767273318723
Newton iteration  0
err =  0.010398975481580032
Newton iteration  1
err =  1.7779014803610412e-05
Newton iteration  2
err =  3.326861907013549e-10
cost at iteration 290 :  8.535569674433505
Newton iteration  0
err =  0.010388166937496652
Newton iteration  1
err =  1.771466463184349e-05
Newton iteration  2
err =  3.263041792635424e-10
cost at iteration 291 :  8.521396294407083
Newton iteration  0
err =  0.010377442323663458
Newton iteration  1
err =  1.7655780113992115e-05
Newton iteration  2
err =  3.2076283763241824e-10
cost at iteration 292 :  8.507247093814817
Newton iteration  0
err =  0.01036680015521094
Newton iteration  1
err =  1.7602170797941442e-05
Newton iteration  2
err =  3.1602080798532066e-10
cost at iteration 293 :  8.49312203361659
Newton iteration  0
err =  0.010356239028215645
Newton iteration  1
err =  1.755366004547435e-05
Newton iteration  2
err =  3.120400275841814e-10
cost at iteration 294 :  8.479021075150364
Newton iteration  0
err =  0.010345757615867103
Newton iteration  1
err =  1.751008421189567e-05
Newton iteration  2
err =  3.0878307932244484e-10
cost at iteration 295 :  8.464944180125306
Newton iteration  0
err =  0.010335354664924816
Newton iteration  1
err =  1.747129189104011e-05
Newton iteration  2
err =  3.062159957712191e-10
cost at iteration 296 :  8.450891310615278
Newton iteration  0
err =  0.010325028992358562
Newton iteration  1
err =  1.74371432240774e-05
Newton iteration  2
err =  3.0430846535554714e-10
cost at iteration 297 :  8.436862429053274
Newton iteration  0
err =  0.010314779482390997
Newton iteration  1
err =  1.7407509270039755e-05
Newton iteration  2
err =  3.0303159886840744e-10
cost at iteration 298 :  8.422857498225401
Newton iteration  0
err =  0.010304605083590454
Newton iteration  1
err =  1.7382271427044022e-05
Newton iteration  2
err =  3.0235782817468916e-10
cost at iteration 299 :  8.408876481265745
Newton iteration  0
err =  0.01029450480635507
Newton iteration  1
err =  1.7361320901541812e-05
Newton iteration  2
err =  3.022615727882552e-10
cost at iteration 300 :  8.394919341651374
Newton iteration  0
err =  0.010284477720448798
Newton iteration  1
err =  1.734455821813535e-05
Newton iteration  2
err =  3.027201948375371e-10
cost at iteration 301 :  8.380986043197575
Newton iteration  0
err =  0.010274522952835
Newton iteration  1
err =  1.7331892789270654e-05
Newton iteration  2
err =  3.0371329287697595e-10
cost at iteration 302 :  8.367076550053328
Newton iteration  0
err =  0.010264639685620037
Newton iteration  1
err =  1.7323242502307707e-05
Newton iteration  2
err =  3.0522088788762664e-10
cost at iteration 303 :  8.353190826697203
Newton iteration  0
err =  0.010254827154192117
Newton iteration  1
err =  1.7318533354047955e-05
Newton iteration  2
err =  3.072268301869809e-10
cost at iteration 304 :  8.339328837933415
Newton iteration  0
err =  0.010245084645508936
Newton iteration  1
err =  1.731769912376214e-05
Newton iteration  2
err =  3.097146622120062e-10
cost at iteration 305 :  8.325490548888231
Newton iteration  0
err =  0.010235411496504617
Newton iteration  1
err =  1.73206810647966e-05
Newton iteration  2
err =  3.126729693447604e-10
cost at iteration 306 :  8.311675925006464
Newton iteration  0
err =  0.010225807092709474
Newton iteration  1
err =  1.732742763759152e-05
Newton iteration  2
err =  3.1608944962370677e-10
cost at iteration 307 :  8.29788493204824
Newton iteration  0
err =  0.010216270866902069
Newton iteration  1
err =  1.733789427018766e-05
Newton iteration  2
err =  3.1995462744244705e-10
cost at iteration 308 :  8.284117536086221
Newton iteration  0
err =  0.01020680229796515
Newton iteration  1
err =  1.7352043132509136e-05
Newton iteration  2
err =  3.242600057048941e-10
cost at iteration 309 :  8.270373703502608
Newton iteration  0
err =  0.010197400909811584
Newton iteration  1
err =  1.7369842956466188e-05
Newton iteration  2
err =  3.290010343147126e-10
cost at iteration 310 :  8.256653400986828
Newton iteration  0
err =  0.010188066270470756
Newton iteration  1
err =  1.7391268860904182e-05
Newton iteration  2
err =  3.3417325115878346e-10
cost at iteration 311 :  8.242956595533146
Newton iteration  0
err =  0.01017879799123705
Newton iteration  1
err =  1.7416302210589806e-05
Newton iteration  2
err =  3.397748357500292e-10
cost at iteration 312 :  8.22928325443864
Newton iteration  0
err =  0.010169595725921105
Newton iteration  1
err =  1.7444930497827835e-05
Newton iteration  2
err =  3.4580517083605715e-10
cost at iteration 313 :  8.215633345301129
Newton iteration  0
err =  0.010160459170286292
Newton iteration  1
err =  1.7477147234072907e-05
Newton iteration  2
err =  3.5226334855039417e-10
cost at iteration 314 :  8.202006836017606
Newton iteration  0
err =  0.010151388061481605
Newton iteration  1
err =  1.751295187589034e-05
Newton iteration  2
err =  3.591553442215056e-10
cost at iteration 315 :  8.188403694782872
Newton iteration  0
err =  0.010142382177609063
Newton iteration  1
err =  1.755234976180241e-05
Newton iteration  2
err =  3.664843495301553e-10
cost at iteration 316 :  8.174823890087946
Newton iteration  0
err =  0.010133441337422337
Newton iteration  1
err =  1.7595352071572173e-05
Newton iteration  2
err =  3.7425804862622527e-10
cost at iteration 317 :  8.161267390719201
Newton iteration  0
err =  0.010124565400060985
Newton iteration  1
err =  1.7641975803739507e-05
Newton iteration  2
err =  3.824829122975472e-10
cost at iteration 318 :  8.147734165757226
Newton iteration  0
err =  0.010115754264902079
Newton iteration  1
err =  1.769224377640475e-05
Newton iteration  2
err =  3.9117123453281395e-10
cost at iteration 319 :  8.13422418457646
Newton iteration  0
err =  0.010107007871518296
Newton iteration  1
err =  1.7746184635732924e-05
Newton iteration  2
err =  4.0033125469969907e-10
cost at iteration 320 :  8.12073741684424
Newton iteration  0
err =  0.010098326199735079
Newton iteration  1
err =  1.7803832894577003e-05
Newton iteration  2
err =  4.0998101829613647e-10
cost at iteration 321 :  8.107273832520715
Newton iteration  0
err =  0.010089709269664026
Newton iteration  1
err =  1.7865228983606987e-05
Newton iteration  2
err =  4.2013163852567323e-10
cost at iteration 322 :  8.09383340185873
Newton iteration  0
err =  0.010081157142043629
Newton iteration  1
err =  1.793041932269758e-05
Newton iteration  2
err =  4.308036900990181e-10
cost at iteration 323 :  8.080416095403713
Newton iteration  0
err =  0.010072669918452992
Newton iteration  1
err =  1.7999456406894193e-05
Newton iteration  2
err =  4.420146513294219e-10
cost at iteration 324 :  8.067021883994007
Newton iteration  0
err =  0.010064247741761732
Newton iteration  1
err =  1.8072398920311412e-05
Newton iteration  2
err =  4.5378584053413096e-10
cost at iteration 325 :  8.05365073876125
Newton iteration  0
err =  0.010055890796595975
Newton iteration  1
err =  1.8149311860283086e-05
Newton iteration  2
err =  4.661409257472278e-10
cost at iteration 326 :  8.040302631131096
Newton iteration  0
err =  0.010047599309983956
Newton iteration  1
err =  1.8230266682994816e-05
Newton iteration  2
err =  4.791053903184e-10
cost at iteration 327 :  8.026977532823976
Newton iteration  0
err =  0.010039373551996699
Newton iteration  1
err =  1.8315341474506054e-05
Newton iteration  2
err =  4.927075792914001e-10
cost at iteration 328 :  8.013675415856108
Newton iteration  0
err =  0.01003121383655408
Newton iteration  1
err =  1.8404621134716773e-05
Newton iteration  2
err =  5.069763919625628e-10
cost at iteration 329 :  8.000396252540833
Newton iteration  0
err =  0.010023120522354248
Newton iteration  1
err =  1.8498197590319304e-05
Newton iteration  2
err =  5.219457044743445e-10
cost at iteration 330 :  7.98714001548986
Newton iteration  0
err =  0.010015094013871174
Newton iteration  1
err =  1.8596170023404392e-05
Newton iteration  2
err =  5.376507411111589e-10
cost at iteration 331 :  7.97390667761529
Newton iteration  0
err =  0.010007134762483173
Newton iteration  1
err =  1.8698645123893623e-05
Newton iteration  2
err =  5.54130822123991e-10
cost at iteration 332 :  7.960696212131149
Newton iteration  0
err =  0.009999243267681233
Newton iteration  1
err =  1.8805737368903218e-05
Newton iteration  2
err =  5.714264511013578e-10
cost at iteration 333 :  7.947508592555636
Newton iteration  0
err =  0.00999142007849219
Newton iteration  1
err =  1.8917569326535872e-05
Newton iteration  2
err =  5.895813670581057e-10
cost at iteration 334 :  7.934343792713496
Newton iteration  0
err =  0.009983665794949054
Newton iteration  1
err =  1.9034271976626107e-05
Newton iteration  2
err =  6.086474896437384e-10
cost at iteration 335 :  7.921201786738643
Newton iteration  0
err =  0.009975981069690464
Newton iteration  1
err =  1.9155985070018674e-05
Newton iteration  2
err =  6.28673868835676e-10
cost at iteration 336 :  7.908082549076849
Newton iteration  0
err =  0.009968366609787881
Newton iteration  1
err =  1.9282857508541858e-05
Newton iteration  2
err =  6.497169604145605e-10
cost at iteration 337 :  7.89498605448894
Newton iteration  0
err =  0.009960823178635428
Newton iteration  1
err =  1.9415047753028568e-05
Newton iteration  2
err =  6.718383401351325e-10
cost at iteration 338 :  7.881912278054166
Newton iteration  0
err =  0.009953351598019265
Newton iteration  1
err =  1.955272426425585e-05
Newton iteration  2
err =  6.951029043824951e-10
cost at iteration 339 :  7.8688611951736505
Newton iteration  0
err =  0.009945952750389186
Newton iteration  1
err =  1.969606597618108e-05
Newton iteration  2
err =  7.195811061109158e-10
cost at iteration 340 :  7.855832781574667
Newton iteration  0
err =  0.009938627581249238
Newton iteration  1
err =  1.9845262796823227e-05
Newton iteration  2
err =  7.45348866507181e-10
cost at iteration 341 :  7.842827013314456
Newton iteration  0
err =  0.0099313771018026
Newton iteration  1
err =  2.0000516146452823e-05
Newton iteration  2
err =  7.724873027428906e-10
cost at iteration 342 :  7.829843866785014
Newton iteration  0
err =  0.009924202391663162
Newton iteration  1
err =  2.016203953044025e-05
Newton iteration  2
err =  8.010859720016777e-10
cost at iteration 343 :  7.81688331871779
Newton iteration  0
err =  0.00991710460196803
Newton iteration  1
err =  2.033005914524305e-05
Newton iteration  2
err =  8.312397433966234e-10
cost at iteration 344 :  7.803945346188948
Newton iteration  0
err =  0.009910084958523632
Newton iteration  1
err =  2.0504814526263714e-05
Newton iteration  2
err =  8.630503207097703e-10
cost at iteration 345 :  7.791029926624875
Newton iteration  0
err =  0.009903144765262523
Newton iteration  1
err =  2.0686559229127302e-05
Newton iteration  2
err =  8.96629958404208e-10
cost at iteration 346 :  7.778137037808137
Newton iteration  0
err =  0.009896285407998912
Newton iteration  1
err =  2.0875561553211083e-05
Newton iteration  2
err =  9.320982069230523e-10
cost at iteration 347 :  7.765266657883698
Newton iteration  0
err =  0.009889508358301418
Newton iteration  1
err =  2.107210529907077e-05
Newton iteration  2
err =  9.695859122376951e-10
cost at iteration 348 :  7.7524187653657135
Newton iteration  0
err =  0.009882815177749332
Newton iteration  1
err =  2.1276490575419525e-05
Newton iteration  2
err =  1.0092307867918985e-09
cost at iteration 349 :  7.739593339144463
Newton iteration  0
err =  0.00987620752248644
Newton iteration  1
err =  2.1489034630325355e-05
Newton iteration  2
err =  1.0511865409069397e-09
cost at iteration 350 :  7.726790358494299
Newton iteration  0
err =  0.009869687147961296
Newton iteration  1
err =  2.1710072730444545e-05
Newton iteration  2
err =  1.0956166410079054e-09
cost at iteration 351 :  7.7140098030812645
Newton iteration  0
err =  0.009863255914100777
Newton iteration  1
err =  2.1939959074371744e-05
Newton iteration  2
err =  1.1426997053618905e-09
cost at iteration 352 :  7.701251652972075
Newton iteration  0
err =  0.009856915790792505
Newton iteration  1
err =  2.217906773881653e-05
Newton iteration  2
err =  1.1926271033226805e-09
cost at iteration 353 :  7.688515888642841
Newton iteration  0
err =  0.009850668863638135
Newton iteration  1
err =  2.242779364812339e-05
Newton iteration  2
err =  1.2456082019461151e-09
cost at iteration 354 :  7.6758024909888745
Newton iteration  0
err =  0.00984451734018273
Newton iteration  1
err =  2.2686553576899238e-05
Newton iteration  2
err =  1.3018685072092917e-09
cost at iteration 355 :  7.663111441334938
Newton iteration  0
err =  0.009838463556459646
Newton iteration  1
err =  2.2955787158191435e-05
Newton iteration  2
err =  1.36165406582935e-09
cost at iteration 356 :  7.65044272144599
Newton iteration  0
err =  0.00983250998390936
Newton iteration  1
err =  2.3235957907955807e-05
Newton iteration  2
err =  1.4252301950942225e-09
cost at iteration 357 :  7.6377963135387805
Newton iteration  0
err =  0.009826659236764242
Newton iteration  1
err =  2.352755423391383e-05
Newton iteration  2
err =  1.49288632533243e-09
cost at iteration 358 :  7.625172200293843
Newton iteration  0
err =  0.009820914079738454
Newton iteration  1
err =  2.383109042637388e-05
Newton iteration  2
err =  1.5649353565276393e-09
cost at iteration 359 :  7.6125703648686915
Newton iteration  0
err =  0.009815277436272347
Newton iteration  1
err =  2.4147107604190455e-05
Newton iteration  2
err =  1.6417168268046813e-09
cost at iteration 360 :  7.599990790911121
Newton iteration  0
err =  0.009809752397008572
Newton iteration  1
err =  2.447617458732318e-05
Newton iteration  2
err =  1.7235993217923515e-09
cost at iteration 361 :  7.587433462574027
Newton iteration  0
err =  0.009804342228818135
Newton iteration  1
err =  2.48188886797379e-05
Newton iteration  2
err =  1.8109843092377124e-09
cost at iteration 362 :  7.57489836453053
Newton iteration  0
err =  0.009799050384092726
Newton iteration  1
err =  2.5175876288010577e-05
Newton iteration  2
err =  1.9043048029681677e-09
cost at iteration 363 :  7.562385481990354
Newton iteration  0
err =  0.009793880510463161
Newton iteration  1
err =  2.5547793369221164e-05
Newton iteration  2
err =  2.0040332416835827e-09
cost at iteration 364 :  7.549894800716963
Newton iteration  0
err =  0.009788836460704373
Newton iteration  1
err =  2.5935325613858137e-05
Newton iteration  2
err =  2.1106783873768954e-09
cost at iteration 365 :  7.537426307045587
Newton iteration  0
err =  0.009783922302982055
Newton iteration  1
err =  2.633918829296073e-05
Newton iteration  2
err =  2.2247944438722137e-09
cost at iteration 366 :  7.52497998790282
Newton iteration  0
err =  0.00977914233113634
Newton iteration  1
err =  2.67601256944044e-05
Newton iteration  2
err =  2.34697987725527e-09
cost at iteration 367 :  7.512555830826592
Newton iteration  0
err =  0.009774501075013096
Newton iteration  1
err =  2.7198910006745193e-05
Newton iteration  2
err =  2.477880909865271e-09
cost at iteration 368 :  7.500153823988005
Newton iteration  0
err =  0.009770003310599204
Newton iteration  1
err =  2.7656339532215763e-05
Newton iteration  2
err =  2.6181961669149398e-09
cost at iteration 369 :  7.487773956213732
Newton iteration  0
err =  0.009765654069783313
Newton iteration  1
err =  2.8133236050470487e-05
Newton iteration  2
err =  2.768677745677034e-09
cost at iteration 370 :  7.475416217010188
Newton iteration  0
err =  0.009761458649468441
Newton iteration  1
err =  2.863044111672924e-05
Newton iteration  2
err =  2.9301326690914217e-09
cost at iteration 371 :  7.463080596588835
Newton iteration  0
err =  0.009757422619622782
Newton iteration  1
err =  2.9148811054707186e-05
Newton iteration  2
err =  3.103427406174945e-09
cost at iteration 372 :  7.45076708589259
Newton iteration  0
err =  0.009753551829872264
Newton iteration  1
err =  2.9689210334377875e-05
Newton iteration  2
err =  3.289485296903913e-09
cost at iteration 373 :  7.43847567662398
Newton iteration  0
err =  0.009749852414065207
Newton iteration  1
err =  3.0252502951988608e-05
Newton iteration  2
err =  3.4892889093365995e-09
cost at iteration 374 :  7.426206361274561
Newton iteration  0
err =  0.009746330791994353
Newton iteration  1
err =  3.083954138060381e-05
Newton iteration  2
err =  3.7038746812260253e-09
cost at iteration 375 :  7.413959133155532
Newton iteration  0
err =  0.009742993667551856
Newton iteration  1
err =  3.145115253959402e-05
Newton iteration  2
err =  3.934329998704516e-09
cost at iteration 376 :  7.4017339864302665
Newton iteration  0
err =  0.009739848021939982
Newton iteration  1
err =  3.2088120135433247e-05
Newton iteration  2
err =  4.181787107736253e-09
cost at iteration 377 :  7.389530916147467
Newton iteration  0
err =  0.009736901100909051
Newton iteration  1
err =  3.275116260255677e-05
Newton iteration  2
err =  4.447408347274089e-09
cost at iteration 378 :  7.377349918276473
Newton iteration  0
err =  0.009734160393787737
Newton iteration  1
err =  3.3440905713660976e-05
Newton iteration  2
err =  4.732369029639533e-09
cost at iteration 379 :  7.36519098974291
Newton iteration  0
err =  0.009731633602747201
Newton iteration  1
err =  3.415784876741146e-05
Newton iteration  2
err =  5.037834855387816e-09
cost at iteration 380 :  7.35305412846606
Newton iteration  0
err =  0.009729328599216587
Newton iteration  1
err =  3.490232305345404e-05
Newton iteration  2
err =  5.364930720597321e-09
cost at iteration 381 :  7.340939333396536
Newton iteration  0
err =  0.009727253364437449
Newton iteration  1
err =  3.5674441092446135e-05
Newton iteration  2
err =  5.7146938102749805e-09
cost at iteration 382 :  7.328846604554824
Newton iteration  0
err =  0.009725415910288651
Newton iteration  1
err =  3.647403490225224e-05
Newton iteration  2
err =  6.088018614060089e-09
cost at iteration 383 :  7.316775943069449
Newton iteration  0
err =  0.009723824175419567
Newton iteration  1
err =  3.730058128833815e-05
Newton iteration  2
err =  6.485585172367059e-09
cost at iteration 384 :  7.304727351215091
Newton iteration  0
err =  0.00972248589117931
Newton iteration  1
err =  3.815311192487166e-05
Newton iteration  2
err =  6.9077601301855975e-09
cost at iteration 385 :  7.292700832449385
Newton iteration  0
err =  0.009721408410165031
Newton iteration  1
err =  3.903010576142496e-05
Newton iteration  2
err =  7.354486022898466e-09
cost at iteration 386 :  7.2806963914472345
Newton iteration  0
err =  0.009720598489301785
Newton iteration  1
err =  3.992936114590643e-05
Newton iteration  2
err =  7.825133778556886e-09
cost at iteration 387 :  7.268714034132369
Newton iteration  0
err =  0.009720062017840793
Newton iteration  1
err =  4.084784502995097e-05
Newton iteration  2
err =  8.318326199454976e-09
cost at iteration 388 :  7.256753767704086
Newton iteration  0
err =  0.009719803678429164
Newton iteration  1
err =  4.178151679406261e-05
Newton iteration  2
err =  8.831727192912407e-09
cost at iteration 389 :  7.244815600657051
Newton iteration  0
err =  0.009719826528935277
Newton iteration  1
err =  4.272512477712948e-05
Newton iteration  2
err =  9.361799844870752e-09
cost at iteration 390 :  7.232899542792949
Newton iteration  0
err =  0.009720131489545437
Newton iteration  1
err =  4.367197458877147e-05
Newton iteration  2
err =  9.90352798832348e-09
cost at iteration 391 :  7.221005605220383
Newton iteration  0
err =  0.009720716719626439
Newton iteration  1
err =  4.461367011296824e-05
Newton iteration  2
err =  1.0450111266433136e-08
Newton iteration  3
err =  1.398617220183724e-14
cost at iteration 392 :  7.209133800339945
Newton iteration  0
err =  0.009721576866879317
Newton iteration  1
err =  4.553983090504712e-05
Newton iteration  2
err =  1.0992655892544408e-08
Newton iteration  3
err =  1.3763130590648683e-14
cost at iteration 393 :  7.19728414181033
Newton iteration  0
err =  0.009722702171973602
Newton iteration  1
err =  4.6437794013292374e-05
Newton iteration  2
err =  1.1519873590516903e-08
Newton iteration  3
err =  1.457052684168393e-14
cost at iteration 394 :  7.1854566444915955
Newton iteration  0
err =  0.009724077413231105
Newton iteration  1
err =  4.729231447959694e-05
Newton iteration  2
err =  1.2017844522160257e-08
Newton iteration  3
err =  1.357974096275976e-14
cost at iteration 395 :  7.173651324358728
Newton iteration  0
err =  0.009725680679794829
Newton iteration  1
err =  4.8085287474875344e-05
Newton iteration  2
err =  1.2469883301809588e-08
Newton iteration  3
err =  1.3941319377051174e-14
cost at iteration 396 :  7.161868198379745
Newton iteration  0
err =  0.009727481969459016
Newton iteration  1
err =  4.879552672915221e-05
Newton iteration  2
err =  1.2856606152104496e-08
Newton iteration  3
err =  1.4278081006324533e-14
cost at iteration 397 :  7.150107284351836
Newton iteration  0
err =  0.009729441619304462
Newton iteration  1
err =  4.9398649099038134e-05
Newton iteration  2
err =  1.3156273224041452e-08
Newton iteration  3
err =  1.3837656796887847e-14
cost at iteration 398 :  7.138368600685554
Newton iteration  0
err =  0.009731508595989277
Newton iteration  1
err =  4.986713382851552e-05
Newton iteration  2
err =  1.3345534722105363e-08
Newton iteration  3
err =  1.43752149383691e-14
cost at iteration 399 :  7.126652166130711
Newton iteration  0
err =  0.009733618698777856
Newton iteration  1
err =  5.017064703507137e-05
Newton iteration  2
err =  1.3400667598194664e-08
Newton iteration  3
err =  1.4233773671138964e-14
cost at iteration 400 :  7.114957999433765
Newton iteration  0
err =  0.009735692763459997
Newton iteration  1
err =  5.0276745788699966e-05
Newton iteration  2
err =  1.3299386142870326e-08
Newton iteration  3
err =  1.4041397769829267e-14
cost at iteration 401 :  7.103286118918815
Newton iteration  0
err =  0.009737634999680512
Newton iteration  1
err =  5.015209941563411e-05
Newton iteration  2
err =  1.3023230273862849e-08
Newton iteration  3
err =  1.3582801165669328e-14
cost at iteration 402 :  7.091636541984913
Newton iteration  0
err =  0.009739331645368735
Newton iteration  1
err =  4.9764384222212626e-05
Newton iteration  2
err =  1.2560409613699563e-08
Newton iteration  3
err =  1.369468174446171e-14
cost at iteration 403 :  7.080009284514986
Newton iteration  0
err =  0.00974065017666878
Newton iteration  1
err =  4.908501598486364e-05
Newton iteration  2
err =  1.1908833874310026e-08
Newton iteration  3
err =  1.4185136414341598e-14
cost at iteration 404 :  7.068404360195021
Newton iteration  0
err =  0.00974143936038501
Newton iteration  1
err =  4.8092874893604476e-05
Newton iteration  2
err =  1.1078844602222444e-08
Newton iteration  3
err =  1.3497578061456434e-14
cost at iteration 405 :  7.0568217797465325
Newton iteration  0
err =  0.009741530466938741
Newton iteration  1
err =  4.6779140441579204e-05
Newton iteration  2
err =  1.0095040104838732e-08
Newton iteration  3
err =  1.4647977986138938e-14
cost at iteration 406 :  7.045261550085091
Newton iteration  0
err =  0.009740739956494174
Newton iteration  1
err =  4.515327368926383e-05
Newton iteration  2
err =  8.996510004716513e-09
cost at iteration 407 :  7.033723673420393
Newton iteration  0
err =  0.009738873893465777
Newton iteration  1
err =  4.3250028071159e-05
Newton iteration  2
err =  7.834996131122178e-09
cost at iteration 408 :  7.022208146329081
Newton iteration  0
err =  0.009735734215135955
Newton iteration  1
err =  4.1137054130954686e-05
Newton iteration  2
err =  6.670923693965339e-09
cost at iteration 409 :  7.010714958833902
Newton iteration  0
err =  0.009731126777683596
Newton iteration  1
err =  3.892199379656642e-05
Newton iteration  2
err =  5.567971502616132e-09
cost at iteration 410 :  6.999244093535542
Newton iteration  0
err =  0.009724870834959463
Newton iteration  1
err =  3.6756576150920496e-05
Newton iteration  2
err =  4.5876523144821494e-09
cost at iteration 411 :  6.987795524845162
Newton iteration  0
err =  0.00971680930961486
Newton iteration  1
err =  3.48328003798658e-05
Newton iteration  2
err =  3.785731734856403e-09
cost at iteration 412 :  6.976369218367619
Newton iteration  0
err =  0.00970681894911114
Newton iteration  1
err =  3.3363641807648354e-05
Newton iteration  2
err =  3.2109520410634428e-09
cost at iteration 413 :  6.964965130480298
Newton iteration  0
err =  0.009694819290242209
Newton iteration  1
err =  3.2542166755905926e-05
Newton iteration  2
err =  2.901885747580299e-09
cost at iteration 414 :  6.953583208143355
Newton iteration  0
err =  0.009680779348461156
Newton iteration  1
err =  3.248566661107678e-05
Newton iteration  2
err =  2.8752668463732654e-09
cost at iteration 415 :  6.942223388960291
Newton iteration  0
err =  0.009664721135660708
Newton iteration  1
err =  3.319231825939447e-05
Newton iteration  2
err =  3.1148330705041826e-09
cost at iteration 416 :  6.930885601492659
Newton iteration  0
err =  0.009646719478402193
Newton iteration  1
err =  3.4541097713293145e-05
Newton iteration  2
err =  3.580748325470575e-09
cost at iteration 417 :  6.919569765812022
Newton iteration  0
err =  0.00962689809315406
Newton iteration  1
err =  3.633522985543137e-05
Newton iteration  2
err =  4.228112859701795e-09
cost at iteration 418 :  6.908275794256148
Newton iteration  0
err =  0.009605422375915721
Newton iteration  1
err =  3.835857242420371e-05
Newton iteration  2
err =  5.012828726538204e-09
cost at iteration 419 :  6.8970035923446975
Newton iteration  0
err =  0.00958248977504741
Newton iteration  1
err =  4.041702839506683e-05
Newton iteration  2
err =  5.888626635918416e-09
cost at iteration 420 :  6.885753059801368
Newton iteration  0
err =  0.009558318856725322
Newton iteration  1
err =  4.235808147410718e-05
Newton iteration  2
err =  6.805323327299809e-09
cost at iteration 421 :  6.8745240916279915
Newton iteration  0
err =  0.00953313821159469
Newton iteration  1
err =  4.4074574504422617e-05
Newton iteration  2
err =  7.710939790023761e-09
cost at iteration 422 :  6.863316579181749
Newton iteration  0
err =  0.009507176206416591
Newton iteration  1
err =  4.550041240367954e-05
Newton iteration  2
err =  8.556158587950486e-09
cost at iteration 423 :  6.852130411211048
Newton iteration  0
err =  0.009480652315767274
Newton iteration  1
err =  4.660319849368785e-05
Newton iteration  2
err =  9.299015084230928e-09
cost at iteration 424 :  6.840965474818703
Newton iteration  0
err =  0.009453770448706099
Newton iteration  1
err =  4.737634589286554e-05
Newton iteration  2
err =  9.908360999795208e-09
cost at iteration 425 :  6.82982165632848
Newton iteration  0
err =  0.009426714382113193
Newton iteration  1
err =  4.783179862966938e-05
Newton iteration  2
err =  1.0365446718949274e-08
Newton iteration  3
err =  1.4050596565451031e-14
cost at iteration 426 :  6.818698842043592
Newton iteration  0
err =  0.009399645171044123
Newton iteration  1
err =  4.7993816202270816e-05
Newton iteration  2
err =  1.0663656162068951e-08
Newton iteration  3
err =  1.3178458760927205e-14
cost at iteration 427 :  6.807596918891882
Newton iteration  0
err =  0.009372700248258958
Newton iteration  1
err =  4.7893952232831266e-05
Newton iteration  2
err =  1.080684707040156e-08
Newton iteration  3
err =  1.4111951610100419e-14
cost at iteration 428 :  6.796515774959711
Newton iteration  0
err =  0.009345993848134834
Newton iteration  1
err =  4.756719491187188e-05
Newton iteration  2
err =  1.0806899474350026e-08
Newton iteration  3
err =  1.4007377762386024e-14
cost at iteration 429 :  6.785455299922468
Newton iteration  0
err =  0.00931961837814476
Newton iteration  1
err =  4.704915065642696e-05
Newton iteration  2
err =  1.0681039611110315e-08
Newton iteration  3
err =  1.3763700315135681e-14
cost at iteration 430 :  6.774415385379738
Newton iteration  0
err =  0.009293646395272564
Newton iteration  1
err =  4.637411048976975e-05
Newton iteration  2
err =  1.0449354696564205e-08
Newton iteration  3
err =  1.403056298483153e-14
cost at iteration 431 :  6.763395925107754
Newton iteration  0
err =  0.00926813290172403
Newton iteration  1
err =  4.557382607919982e-05
Newton iteration  2
err =  1.0132746137814827e-08
Newton iteration  3
err =  1.3485393971588039e-14
cost at iteration 432 :  6.75239681523993
Newton iteration  0
err =  0.009243117741110946
Newton iteration  1
err =  4.467682952317102e-05
Newton iteration  2
err =  9.751406700657248e-09
cost at iteration 433 :  6.741417954387316
Newton iteration  0
err =  0.009218627939091718
Newton iteration  1
err =  4.3708149991526546e-05
Newton iteration  2
err =  9.323814346922373e-09
cost at iteration 434 :  6.730459243708893
Newton iteration  0
err =  0.009194679887524499
Newton iteration  1
err =  4.268930491848253e-05
Newton iteration  2
err =  8.866163044919178e-09
cost at iteration 435 :  6.719520586942183
Newton iteration  0
err =  0.009171281314503777
Newton iteration  1
err =  4.1638469036240785e-05
Newton iteration  2
err =  8.392134460874286e-09
cost at iteration 436 :  6.708601890400773
Newton iteration  0
err =  0.009148433014795456
Newton iteration  1
err =  4.0570748282000334e-05
Newton iteration  2
err =  7.912911851157044e-09
cost at iteration 437 :  6.697703062947585
Newton iteration  0
err =  0.009126130337977203
Newton iteration  1
err =  3.949850589184493e-05
Newton iteration  2
err =  7.4373476352486464e-09
cost at iteration 438 :  6.68682401594873
Newton iteration  0
err =  0.009104364445456876
Newton iteration  1
err =  3.843170444874152e-05
Newton iteration  2
err =  6.972215193619555e-09
cost at iteration 439 :  6.675964663213248
Newton iteration  0
err =  0.009083123356465854
Newton iteration  1
err =  3.737824028682545e-05
Newton iteration  2
err =  6.5225037621976025e-09
cost at iteration 440 :  6.665124920922531
Newton iteration  0
err =  0.009062392806559407
Newton iteration  1
err =  3.634425600688165e-05
Newton iteration  2
err =  6.0917085020167406e-09
cost at iteration 441 :  6.654304707552929
Newton iteration  0
err =  0.009042156943816111
Newton iteration  1
err =  3.53344234933217e-05
Newton iteration  2
err =  5.682110899608971e-09
cost at iteration 442 :  6.643503943793962
Newton iteration  0
err =  0.00902239888649081
Newton iteration  1
err =  3.435219434994371e-05
Newton iteration  2
err =  5.2950251633383934e-09
cost at iteration 443 :  6.632722552464204
Newton iteration  0
err =  0.009003101164597987
Newton iteration  1
err =  3.340001759933598e-05
Newton iteration  2
err =  4.931014188476379e-09
cost at iteration 444 :  6.62196045842637
Newton iteration  0
err =  0.008984246065244207
Newton iteration  1
err =  3.2479526286996525e-05
Newton iteration  2
err =  4.590070210266184e-09
cost at iteration 445 :  6.611217588502981
Newton iteration  0
err =  0.008965815898659309
Newton iteration  1
err =  3.159169560556029e-05
Newton iteration  2
err =  4.27176425636873e-09
cost at iteration 446 :  6.600493871393054
Newton iteration  0
err =  0.008947793200256985
Newton iteration  1
err =  3.073697562830319e-05
Newton iteration  2
err =  3.975364624295658e-09
cost at iteration 447 :  6.5897892375914
Newton iteration  0
err =  0.008930160880689299
Newton iteration  1
err =  2.9915401815835417e-05
Newton iteration  2
err =  3.6999362129556096e-09
cost at iteration 448 :  6.579103619309782
Newton iteration  0
err =  0.008912902334462796
Newton iteration  1
err =  2.9126686362901595e-05
Newton iteration  2
err =  3.4444105194300977e-09
cost at iteration 449 :  6.568436950401283
Newton iteration  0
err =  0.008896001515751624
Newton iteration  1
err =  2.837029322317834e-05
Newton iteration  2
err =  3.207646777244422e-09
cost at iteration 450 :  6.5577891662876535
Newton iteration  0
err =  0.008879442988206554
Newton iteration  1
err =  2.764549934571673e-05
Newton iteration  2
err =  2.9884727518469274e-09
cost at iteration 451 :  6.547160203889558
Newton iteration  0
err =  0.008863211954658834
Newton iteration  1
err =  2.695144436831966e-05
Newton iteration  2
err =  2.785716011732653e-09
cost at iteration 452 :  6.536550001560378
Newton iteration  0
err =  0.008847294271154552
Newton iteration  1
err =  2.6287170691975027e-05
Newton iteration  2
err =  2.5982305699336286e-09
cost at iteration 453 :  6.525958499023008
Newton iteration  0
err =  0.008831676449089276
Newton iteration  1
err =  2.565165560392926e-05
Newton iteration  2
err =  2.4249103937957684e-09
cost at iteration 454 :  6.515385637309393
Newton iteration  0
err =  0.008816345648428977
Newton iteration  1
err =  2.5043836842105253e-05
Newton iteration  2
err =  2.264696795784614e-09
cost at iteration 455 :  6.504831358704022
Newton iteration  0
err =  0.00880128966424696
Newton iteration  1
err =  2.446263277779022e-05
Newton iteration  2
err =  2.116591407039426e-09
cost at iteration 456 :  6.494295606689477
Newton iteration  0
err =  0.008786496908508004
Newton iteration  1
err =  2.3906958207483914e-05
Newton iteration  2
err =  1.979656308709132e-09
cost at iteration 457 :  6.483778325895214
Newton iteration  0
err =  0.008771956388557135
Newton iteration  1
err =  2.3375736568348085e-05
Newton iteration  2
err =  1.8530134861426442e-09
cost at iteration 458 :  6.473279462048986
Newton iteration  0
err =  0.008757657683433683
Newton iteration  1
err =  2.286790925154169e-05
Newton iteration  2
err =  1.7358485554588326e-09
cost at iteration 459 :  6.462798961931089
Newton iteration  0
err =  0.008743590918807293
Newton iteration  1
err =  2.238244257825825e-05
Newton iteration  2
err =  1.6274048275969789e-09
cost at iteration 460 :  6.4523367733305275
Newton iteration  0
err =  0.008729746741415662
Newton iteration  1
err =  2.1918332888598976e-05
Newton iteration  2
err =  1.5269863735847087e-09
cost at iteration 461 :  6.4418928450042365
Newton iteration  0
err =  0.008716116293235576
Newton iteration  1
err =  2.147461012897602e-05
Newton iteration  2
err =  1.43394844118948e-09
cost at iteration 462 :  6.431467126637739
Newton iteration  0
err =  0.00870269118607923
Newton iteration  1
err =  2.1050340236654653e-05
Newton iteration  2
err =  1.3476994690272858e-09
cost at iteration 463 :  6.421059568808557
Newton iteration  0
err =  0.008689463476587837
Newton iteration  1
err =  2.0644626577187454e-05
Newton iteration  2
err =  1.2676956135162848e-09
cost at iteration 464 :  6.410670122951312
Newton iteration  0
err =  0.0086764256420704
Newton iteration  1
err =  2.025661063570987e-05
Newton iteration  2
err =  1.1934382258702844e-09
cost at iteration 465 :  6.4002987413243035
Newton iteration  0
err =  0.00866357055724329
Newton iteration  1
err =  1.988547213042698e-05
Newton iteration  2
err =  1.124469396643642e-09
cost at iteration 466 :  6.389945376978945
Newton iteration  0
err =  0.008650891471888219
Newton iteration  1
err =  1.953042867864067e-05
Newton iteration  2
err =  1.0603704807272663e-09
cost at iteration 467 :  6.379609983729496
Newton iteration  0
err =  0.00863838198954823
Newton iteration  1
err =  1.9190735124584782e-05
Newton iteration  2
err =  1.0007571372939796e-09
cost at iteration 468 :  6.369292516125287
Newton iteration  0
err =  0.008626036047407462
Newton iteration  1
err =  1.8865682619561323e-05
Newton iteration  2
err =  9.452783660656187e-10
cost at iteration 469 :  6.35899292942415
Newton iteration  0
err =  0.008613847896995182
Newton iteration  1
err =  1.855459751285212e-05
Newton iteration  2
err =  8.936125954387467e-10
cost at iteration 470 :  6.348711179567033
Newton iteration  0
err =  0.008601812086249364
Newton iteration  1
err =  1.8256840121189622e-05
Newton iteration  2
err =  8.454640896520835e-10
cost at iteration 471 :  6.338447223154226
Newton iteration  0
err =  0.008589923442326298
Newton iteration  1
err =  1.7971803416188874e-05
Newton iteration  2
err =  8.005660369466255e-10
cost at iteration 472 :  6.328201017422577
Newton iteration  0
err =  0.008578177055572738
Newton iteration  1
err =  1.7698911661347036e-05
Newton iteration  2
err =  7.586692259153827e-10
cost at iteration 473 :  6.317972520223959
Newton iteration  0
err =  0.008566568264410546
Newton iteration  1
err =  1.7437619027418947e-05
Newton iteration  2
err =  7.195488910173281e-10
cost at iteration 474 :  6.307761690004833
Newton iteration  0
err =  0.008555092641115317
Newton iteration  1
err =  1.7187408210648528e-05
Newton iteration  2
err =  6.829977444390831e-10
cost at iteration 475 :  6.297568485786663
Newton iteration  0
err =  0.008543745978498404
Newton iteration  1
err =  1.6947789064823022e-05
Newton iteration  2
err =  6.488267847339164e-10
cost at iteration 476 :  6.287392867147724
Newton iteration  0
err =  0.00853252427738159
Newton iteration  1
err =  1.6718297263730075e-05
Newton iteration  2
err =  6.168620957784158e-10
cost at iteration 477 :  6.277234794205239
Newton iteration  0
err =  0.008521423734867509
Newton iteration  1
err =  1.6498492997510444e-05
Newton iteration  2
err =  5.869444113357366e-10
cost at iteration 478 :  6.2670942275986805
Newton iteration  0
err =  0.008510440733386096
Newton iteration  1
err =  1.6287959716786377e-05
Newton iteration  2
err =  5.58928071419119e-10
cost at iteration 479 :  6.256971128474109
Newton iteration  0
err =  0.008499571830323493
Newton iteration  1
err =  1.6086302922318067e-05
Newton iteration  2
err =  5.326785626901691e-10
cost at iteration 480 :  6.246865458468624
Newton iteration  0
err =  0.008488813748430869
Newton iteration  1
err =  1.589314900556803e-05
Newton iteration  2
err =  5.080745088045093e-10
cost at iteration 481 :  6.236777179696367
Newton iteration  0
err =  0.008478163366732912
Newton iteration  1
err =  1.5708144142116658e-05
Newton iteration  2
err =  4.850015411856928e-10
cost at iteration 482 :  6.226706254734282
Newton iteration  0
err =  0.00846761771211579
Newton iteration  1
err =  1.5530953237275977e-05
Newton iteration  2
err =  4.63356813103022e-10
cost at iteration 483 :  6.2166526466096155
Newton iteration  0
err =  0.00845717395130059
Newton iteration  1
err =  1.5361258927872325e-05
Newton iteration  2
err =  4.4304467900862085e-10
cost at iteration 484 :  6.206616318786624
Newton iteration  0
err =  0.008446829383521783
Newton iteration  1
err =  1.51987606308227e-05
Newton iteration  2
err =  4.239788635456096e-10
cost at iteration 485 :  6.196597235155458
Newton iteration  0
err =  0.008436581433458499
Newton iteration  1
err =  1.5043173642220036e-05
Newton iteration  2
err =  4.0607817567795717e-10
cost at iteration 486 :  6.186595360020166
Newton iteration  0
err =  0.008426427644746807
Newton iteration  1
err =  1.4894228295865868e-05
Newton iteration  2
err =  3.8926952467193375e-10
cost at iteration 487 :  6.176610658088096
Newton iteration  0
err =  0.00841636567388891
Newton iteration  1
err =  1.4751669153511225e-05
Newton iteration  2
err =  3.734847279917522e-10
cost at iteration 488 :  6.166643094459568
Newton iteration  0
err =  0.008406393284434764
Newton iteration  1
err =  1.4615254252981814e-05
Newton iteration  2
err =  3.5866158274829005e-10
cost at iteration 489 :  6.156692634617778
Newton iteration  0
err =  0.008396508341696638
Newton iteration  1
err =  1.448475439755905e-05
Newton iteration  2
err =  3.447423737542744e-10
cost at iteration 490 :  6.146759244419429
Newton iteration  0
err =  0.008386708807608539
Newton iteration  1
err =  1.4359952477180007e-05
Newton iteration  2
err =  3.316749170462532e-10
cost at iteration 491 :  6.136842890085618
Newton iteration  0
err =  0.008376992736070019
Newton iteration  1
err =  1.4240642847812438e-05
Newton iteration  2
err =  3.194115629701764e-10
cost at iteration 492 :  6.126943538193289
Newton iteration  0
err =  0.008367358268459696
Newton iteration  1
err =  1.4126630729310577e-05
Newton iteration  2
err =  3.079061333481607e-10
cost at iteration 493 :  6.117061155666889
Newton iteration  0
err =  0.008357803629449125
Newton iteration  1
err =  1.4017731648684683e-05
Newton iteration  2
err =  2.971182664492504e-10
cost at iteration 494 :  6.107195709770312
Newton iteration  0
err =  0.008348327123164852
Newton iteration  1
err =  1.391377091951684e-05
Newton iteration  2
err =  2.870110165981683e-10
cost at iteration 495 :  6.097347168099662
Newton iteration  0
err =  0.008338927129416925
Newton iteration  1
err =  1.3814583144599115e-05
Newton iteration  2
err =  2.775494303438451e-10
cost at iteration 496 :  6.087515498575674
Newton iteration  0
err =  0.008329602100302473
Newton iteration  1
err =  1.3720011752742464e-05
Newton iteration  2
err =  2.687010170108033e-10
cost at iteration 497 :  6.077700669436927
Newton iteration  0
err =  0.008320350556976226
Newton iteration  1
err =  1.3629908567620305e-05
Newton iteration  2
err =  2.6043770035983574e-10
cost at iteration 498 :  6.06790264923323
Newton iteration  0
err =  0.0083111710865116
Newton iteration  1
err =  1.3544133397933361e-05
Newton iteration  2
err =  2.527321151711883e-10
cost at iteration 499 :  6.058121406819125
Newton iteration  0
err =  0.00830206233922868
Newton iteration  1
err =  1.3462553650346539e-05
Newton iteration  2
err =  2.455595922548239e-10
cost at iteration 500 :  6.048356911348104
Newton iteration  0
err =  0.008293023025738679
Newton iteration  1
err =  1.3385043980063087e-05
Newton iteration  2
err =  2.3889797742188397e-10
cost at iteration 501 :  6.038609132266453
Newton iteration  0
err =  0.008284051914680703
Newton iteration  1
err =  1.3311485942530303e-05
Newton iteration  2
err =  2.327258024430791e-10
cost at iteration 502 :  6.028878039307869
Newton iteration  0
err =  0.00827514783016297
Newton iteration  1
err =  1.324176768565226e-05
Newton iteration  2
err =  2.2702396529257512e-10
cost at iteration 503 :  6.019163602488128
Newton iteration  0
err =  0.008266309649628545
Newton iteration  1
err =  1.31757836491911e-05
Newton iteration  2
err =  2.2177551338124054e-10
cost at iteration 504 :  6.009465792099981
Newton iteration  0
err =  0.008257536301712302
Newton iteration  1
err =  1.311343428739471e-05
Newton iteration  2
err =  2.1696341770016532e-10
cost at iteration 505 :  5.999784578708088
Newton iteration  0
err =  0.008248826764298986
Newton iteration  1
err =  1.3054625813108487e-05
Newton iteration  2
err =  2.1257229809326458e-10
cost at iteration 506 :  5.990119933144624
Newton iteration  0
err =  0.00824018006266088
Newton iteration  1
err =  1.2999269947508895e-05
Newton iteration  2
err =  2.0858964754675e-10
cost at iteration 507 :  5.980471826504749
Newton iteration  0
err =  0.008231595267720217
Newton iteration  1
err =  1.2947283698630873e-05
Newton iteration  2
err =  2.0500130161193352e-10
cost at iteration 508 :  5.970840230142147
Newton iteration  0
err =  0.00822307149442724
Newton iteration  1
err =  1.2898589142865036e-05
Newton iteration  2
err =  2.0179570040141952e-10
cost at iteration 509 :  5.961225115665147
Newton iteration  0
err =  0.008214607900230295
Newton iteration  1
err =  1.2853113229292838e-05
Newton iteration  2
err =  1.9896095202508026e-10
cost at iteration 510 :  5.9516264549327325
Newton iteration  0
err =  0.008206203683628625
Newton iteration  1
err =  1.2810787592554716e-05
Newton iteration  2
err =  1.9648758284656608e-10
cost at iteration 511 :  5.942044220050898
Newton iteration  0
err =  0.008197858082843155
Newton iteration  1
err =  1.2771548375813476e-05
Newton iteration  2
err =  1.9436502540704425e-10
cost at iteration 512 :  5.932478383368895
Newton iteration  0
err =  0.00818957037456236
Newton iteration  1
err =  1.2735336074268272e-05
Newton iteration  2
err =  1.9258555058362096e-10
cost at iteration 513 :  5.922928917476177
Newton iteration  0
err =  0.008181339872747042
Newton iteration  1
err =  1.2702095380704025e-05
Newton iteration  2
err =  1.9113932040433875e-10
cost at iteration 514 :  5.913395795198665
Newton iteration  0
err =  0.008173165927550229
Newton iteration  1
err =  1.2671775044791708e-05
Newton iteration  2
err =  1.900179976653705e-10
cost at iteration 515 :  5.903878989596287
Newton iteration  0
err =  0.008165047924262675
Newton iteration  1
err =  1.2644327745162479e-05
Newton iteration  2
err =  1.892138001167686e-10
cost at iteration 516 :  5.894378473959487
Newton iteration  0
err =  0.008156985282420606
Newton iteration  1
err =  1.2619709965349516e-05
Newton iteration  2
err =  1.8872181034163823e-10
cost at iteration 517 :  5.884894221806873
Newton iteration  0
err =  0.008148977454812233
Newton iteration  1
err =  1.2597881882418634e-05
Newton iteration  2
err =  1.8853274024777358e-10
cost at iteration 518 :  5.875426206882051
Newton iteration  0
err =  0.008141023926771343
Newton iteration  1
err =  1.2578807265819067e-05
Newton iteration  2
err =  1.8864247169331393e-10
cost at iteration 519 :  5.86597440315173
Newton iteration  0
err =  0.008133124215322533
Newton iteration  1
err =  1.2562453380909902e-05
Newton iteration  2
err =  1.890441721183539e-10
cost at iteration 520 :  5.856538784802926
Newton iteration  0
err =  0.00812527786853264
Newton iteration  1
err =  1.254879089948866e-05
Newton iteration  2
err =  1.8973303071631942e-10
cost at iteration 521 :  5.847119326240894
Newton iteration  0
err =  0.008117484464797473
Newton iteration  1
err =  1.253779382312144e-05
Newton iteration  2
err =  1.9070385625075613e-10
cost at iteration 522 :  5.83771600208676
Newton iteration  0
err =  0.008109743612354246
Newton iteration  1
err =  1.2529439407198066e-05
Newton iteration  2
err =  1.9195337192986783e-10
cost at iteration 523 :  5.828328787175868
Newton iteration  0
err =  0.008102054948624145
Newton iteration  1
err =  1.2523708093553935e-05
Newton iteration  2
err =  1.9347627094272422e-10
cost at iteration 524 :  5.818957656555715
Newton iteration  0
err =  0.008094418139791814
Newton iteration  1
err =  1.2520583456278075e-05
Newton iteration  2
err =  1.9527082172051233e-10
cost at iteration 525 :  5.809602585484208
Newton iteration  0
err =  0.008086832880349774
Newton iteration  1
err =  1.252005214323651e-05
Newton iteration  2
err =  1.973338657118097e-10
cost at iteration 526 :  5.800263549428046
Newton iteration  0
err =  0.008079298892687463
Newton iteration  1
err =  1.2522103827768975e-05
Newton iteration  2
err =  1.9966296784063297e-10
cost at iteration 527 :  5.790940524061151
Newton iteration  0
err =  0.008071815926739475
Newton iteration  1
err =  1.2526731168988648e-05
Newton iteration  2
err =  2.0225740056897487e-10
cost at iteration 528 :  5.781633485263315
Newton iteration  0
err =  0.008064383759728776
Newton iteration  1
err =  1.2533929774933918e-05
Newton iteration  2
err =  2.051162981861608e-10
cost at iteration 529 :  5.7723424091189
Newton iteration  0
err =  0.008057002195855688
Newton iteration  1
err =  1.2543698169531945e-05
Newton iteration  2
err =  2.0823880407221202e-10
cost at iteration 530 :  5.763067271915612
Newton iteration  0
err =  0.00804967106607208
Newton iteration  1
err =  1.2556037765337035e-05
Newton iteration  2
err =  2.1162763254148652e-10
cost at iteration 531 :  5.753808050143435
Newton iteration  0
err =  0.008042390227981545
Newton iteration  1
err =  1.2570952845141713e-05
Newton iteration  2
err =  2.1528023440611236e-10
cost at iteration 532 :  5.744564720493715
Newton iteration  0
err =  0.008035159565601333
Newton iteration  1
err =  1.2588450540062183e-05
Newton iteration  2
err =  2.19201662556367e-10
cost at iteration 533 :  5.735337259858111
Newton iteration  0
err =  0.008027978989361778
Newton iteration  1
err =  1.2608540819728492e-05
Newton iteration  2
err =  2.2339403972087974e-10
cost at iteration 534 :  5.726125645328157
Newton iteration  0
err =  0.008020848435981729
Newton iteration  1
err =  1.2631236481362172e-05
Newton iteration  2
err =  2.2785938479751314e-10
cost at iteration 535 :  5.716929854194277
Newton iteration  0
err =  0.00801376786849668
Newton iteration  1
err =  1.265655314609669e-05
Newton iteration  2
err =  2.3260337412464485e-10
cost at iteration 536 :  5.707749863945526
Newton iteration  0
err =  0.00800673727625705
Newton iteration  1
err =  1.2684509253905609e-05
Newton iteration  2
err =  2.3763050802503653e-10
cost at iteration 537 :  5.698585652268866
Newton iteration  0
err =  0.007999756675014454
Newton iteration  1
err =  1.2715126070717521e-05
Newton iteration  2
err =  2.4294517762429433e-10
cost at iteration 538 :  5.689437197049102
Newton iteration  0
err =  0.007992826106975102
Newton iteration  1
err =  1.2748427687152409e-05
Newton iteration  2
err =  2.485554630297188e-10
cost at iteration 539 :  5.680304476368519
Newton iteration  0
err =  0.007985945641007066
Newton iteration  1
err =  1.2784441030841315e-05
Newton iteration  2
err =  2.5446898429135946e-10
cost at iteration 540 :  5.671187468506765
Newton iteration  0
err =  0.007979115372743435
Newton iteration  1
err =  1.2823195874343547e-05
Newton iteration  2
err =  2.606926238559937e-10
cost at iteration 541 :  5.662086151940796
Newton iteration  0
err =  0.00797233542483054
Newton iteration  1
err =  1.2864724847016013e-05
Newton iteration  2
err =  2.672353276259297e-10
cost at iteration 542 :  5.65300050534515
Newton iteration  0
err =  0.007965605947178054
Newton iteration  1
err =  1.2909063450312273e-05
Newton iteration  2
err =  2.74107773722818e-10
cost at iteration 543 :  5.643930507591898
Newton iteration  0
err =  0.00795892711719301
Newton iteration  1
err =  1.2956250069902644e-05
Newton iteration  2
err =  2.813205765179934e-10
cost at iteration 544 :  5.634876137751178
Newton iteration  0
err =  0.007952299140198101
Newton iteration  1
err =  1.3006325998573172e-05
Newton iteration  2
err =  2.888849659357409e-10
cost at iteration 545 :  5.625837375091374
Newton iteration  0
err =  0.007945722249692583
Newton iteration  1
err =  1.3059335445839894e-05
Newton iteration  2
err =  2.9681386544520436e-10
cost at iteration 546 :  5.616814199079996
Newton iteration  0
err =  0.007939196707749653
Newton iteration  1
err =  1.3115325558048631e-05
Newton iteration  2
err =  3.0512172066967435e-10
cost at iteration 547 :  5.607806589383843
Newton iteration  0
err =  0.007932722805532439
Newton iteration  1
err =  1.3174346434960832e-05
Newton iteration  2
err =  3.1382175503314044e-10
cost at iteration 548 :  5.598814525870272
Newton iteration  0
err =  0.007926300863622006
Newton iteration  1
err =  1.3236451140286027e-05
Newton iteration  2
err =  3.229308888063699e-10
cost at iteration 549 :  5.589837988607559
Newton iteration  0
err =  0.00791993123262582
Newton iteration  1
err =  1.3301695717980074e-05
Newton iteration  2
err =  3.324651266717776e-10
cost at iteration 550 :  5.580876957866417
Newton iteration  0
err =  0.007913614293626985
Newton iteration  1
err =  1.3370139197894228e-05
Newton iteration  2
err =  3.4244398906814325e-10
cost at iteration 551 :  5.571931414120602
Newton iteration  0
err =  0.007907350458765837
Newton iteration  1
err =  1.3441843600033422e-05
Newton iteration  2
err =  3.5288486788560735e-10
cost at iteration 552 :  5.563001338048367
Newton iteration  0
err =  0.007901140171855614
Newton iteration  1
err =  1.3516873936999256e-05
Newton iteration  2
err =  3.638091528859508e-10
cost at iteration 553 :  5.554086710534076
Newton iteration  0
err =  0.007894983908919864
Newton iteration  1
err =  1.3595298198831245e-05
Newton iteration  2
err =  3.7523779567046767e-10
cost at iteration 554 :  5.545187512669055
Newton iteration  0
err =  0.007888882178938204
Newton iteration  1
err =  1.367718734430061e-05
Newton iteration  2
err =  3.8719343013523494e-10
cost at iteration 555 :  5.536303725753673
Newton iteration  0
err =  0.007882835524427284
Newton iteration  1
err =  1.3762615266295192e-05
Newton iteration  2
err =  3.9970136935238735e-10
cost at iteration 556 :  5.527435331298962
Newton iteration  0
err =  0.007876844522166226
Newton iteration  1
err =  1.3851658759751167e-05
Newton iteration  2
err =  4.1278758662906064e-10
cost at iteration 557 :  5.51858231102845
Newton iteration  0
err =  0.007870909783886219
Newton iteration  1
err =  1.3944397461364093e-05
Newton iteration  2
err =  4.2647712264472146e-10
cost at iteration 558 :  5.509744646879991
Newton iteration  0
err =  0.007865031957033157
Newton iteration  1
err =  1.4040913784978192e-05
Newton iteration  2
err =  4.407999286588591e-10
cost at iteration 559 :  5.500922321008218
Newton iteration  0
err =  0.007859211725439227
Newton iteration  1
err =  1.414129282629913e-05
Newton iteration  2
err =  4.5578684482680946e-10
cost at iteration 560 :  5.492115315786662
Newton iteration  0
err =  0.007853449810081781
Newton iteration  1
err =  1.4245622247571616e-05
Newton iteration  2
err =  4.714686960569678e-10
cost at iteration 561 :  5.483323613810045
Newton iteration  0
err =  0.007847746969824118
Newton iteration  1
err =  1.4353992144917534e-05
Newton iteration  2
err =  4.878792651897743e-10
cost at iteration 562 :  5.474547197896926
Newton iteration  0
err =  0.007842104002147191
Newton iteration  1
err =  1.4466494864044308e-05
Newton iteration  2
err =  5.05053755354408e-10
cost at iteration 563 :  5.465786051092438
Newton iteration  0
err =  0.007836521743773838
Newton iteration  1
err =  1.4583224796828617e-05
Newton iteration  2
err =  5.230283390883619e-10
cost at iteration 564 :  5.457040156671103
Newton iteration  0
err =  0.007831001071461322
Newton iteration  1
err =  1.4704278134402176e-05
Newton iteration  2
err =  5.418417227525245e-10
cost at iteration 565 :  5.448309498139927
Newton iteration  0
err =  0.007825542902524175
Newton iteration  1
err =  1.4829752559794658e-05
Newton iteration  2
err =  5.615346262168241e-10
cost at iteration 566 :  5.439594059241327
Newton iteration  0
err =  0.007820148195435493
Newton iteration  1
err =  1.4959746901534555e-05
Newton iteration  2
err =  5.821485959323693e-10
cost at iteration 567 :  5.430893823956671
Newton iteration  0
err =  0.007814817950382491
Newton iteration  1
err =  1.5094360711366724e-05
Newton iteration  2
err =  6.037274939976459e-10
cost at iteration 568 :  5.422208776509639
Newton iteration  0
err =  0.007809553209530245
Newton iteration  1
err =  1.5233693778120226e-05
Newton iteration  2
err =  6.263164951838951e-10
cost at iteration 569 :  5.413538901369829
Newton iteration  0
err =  0.0078043550574769666
Newton iteration  1
err =  1.537784555891436e-05
Newton iteration  2
err =  6.499626483195095e-10
cost at iteration 570 :  5.404884183256695
Newton iteration  0
err =  0.007799224621228297
Newton iteration  1
err =  1.5526914508877222e-05
Newton iteration  2
err =  6.74715171860191e-10
cost at iteration 571 :  5.3962446071432195
Newton iteration  0
err =  0.00779416307028947
Newton iteration  1
err =  1.5680997316618e-05
Newton iteration  2
err =  7.006223676273083e-10
cost at iteration 572 :  5.387620158260255
Newton iteration  0
err =  0.007789171616364521
Newton iteration  1
err =  1.584018800747182e-05
Newton iteration  2
err =  7.277369282406848e-10
cost at iteration 573 :  5.379010822100767
Newton iteration  0
err =  0.007784251512838803
Newton iteration  1
err =  1.600457691846071e-05
Newton iteration  2
err =  7.561103899674247e-10
cost at iteration 574 :  5.370416584423978
Newton iteration  0
err =  0.007779404054052197
Newton iteration  1
err =  1.6174249511360867e-05
Newton iteration  2
err =  7.857960469032893e-10
cost at iteration 575 :  5.361837431260335
Newton iteration  0
err =  0.007774630574098501
Newton iteration  1
err =  1.6349285014340907e-05
Newton iteration  2
err =  8.168463552946523e-10
cost at iteration 576 :  5.353273348915815
Newton iteration  0
err =  0.007769932445300343
Newton iteration  1
err =  1.6529754863235735e-05
Newton iteration  2
err =  8.493155812188167e-10
cost at iteration 577 :  5.344724323977165
Newton iteration  0
err =  0.007765311076140512
Newton iteration  1
err =  1.671572092019712e-05
Newton iteration  2
err =  8.832562938235196e-10
cost at iteration 578 :  5.336190343316463
Newton iteration  0
err =  0.007760767908698464
Newton iteration  1
err =  1.690723343824332e-05
Newton iteration  2
err =  9.187191244633525e-10
cost at iteration 579 :  5.327671394096573
Newton iteration  0
err =  0.007756304415314915
Newton iteration  1
err =  1.710432874518552e-05
Newton iteration  2
err =  9.557545206308228e-10
cost at iteration 580 :  5.3191674637760435
Newton iteration  0
err =  0.007751922094658117
Newton iteration  1
err =  1.7307026609910238e-05
Newton iteration  2
err =  9.94408760495309e-10
cost at iteration 581 :  5.310678540114642
Newton iteration  0
err =  0.007747622466691406
Newton iteration  1
err =  1.751532725049221e-05
Newton iteration  2
err =  1.0347237832010814e-09
cost at iteration 582 :  5.302204611178359
Newton iteration  0
err =  0.007743407066873933
Newton iteration  1
err =  1.7729207949549903e-05
Newton iteration  2
err =  1.0767381359761941e-09
cost at iteration 583 :  5.293745665345269
Newton iteration  0
err =  0.007739277438887648
Newton iteration  1
err =  1.794861922410038e-05
Newton iteration  2
err =  1.1204816258772703e-09
cost at iteration 584 :  5.285301691310395
Newton iteration  0
err =  0.007735235126307767
Newton iteration  1
err =  1.817348051365307e-05
Newton iteration  2
err =  1.165976730542969e-09
cost at iteration 585 :  5.276872678091341
Newton iteration  0
err =  0.007731281662498428
Newton iteration  1
err =  1.840367531779366e-05
Newton iteration  2
err =  1.213234754305661e-09
cost at iteration 586 :  5.268458615033797
Newton iteration  0
err =  0.007727418558806001
Newton iteration  1
err =  1.8639045741658906e-05
Newton iteration  2
err =  1.2622550170965943e-09
cost at iteration 587 :  5.260059491816233
Newton iteration  0
err =  0.007723647290796462
Newton iteration  1
err =  1.88793863929605e-05
Newton iteration  2
err =  1.3130202873279214e-09
cost at iteration 588 :  5.25167529845525
Newton iteration  0
err =  0.007719969281987918
Newton iteration  1
err =  1.9124437562421617e-05
Newton iteration  2
err =  1.3654949555959318e-09
cost at iteration 589 :  5.243306025310007
Newton iteration  0
err =  0.007716385885252525
Newton iteration  1
err =  1.9373877643156313e-05
Newton iteration  2
err =  1.4196215305525536e-09
cost at iteration 590 :  5.234951663087291
Newton iteration  0
err =  0.007712898360883163
Newton iteration  1
err =  1.962731471352546e-05
Newton iteration  2
err =  1.475317226185749e-09
cost at iteration 591 :  5.226612202844605
Newton iteration  0
err =  0.00770950785168892
Newton iteration  1
err =  1.988427726035643e-05
Newton iteration  2
err =  1.5324675369995589e-09
cost at iteration 592 :  5.218287635994532
Newton iteration  0
err =  0.007706215353939754
Newton iteration  1
err =  2.014420396946258e-05
Newton iteration  2
err =  1.5909243586912085e-09
cost at iteration 593 :  5.209977954307041
Newton iteration  0
err =  0.007703021684454247
Newton iteration  1
err =  2.0406432573351692e-05
Newton iteration  2
err =  1.6504995100852887e-09
cost at iteration 594 :  5.201683149911779
Newton iteration  0
err =  0.007699927442642521
Newton iteration  1
err =  2.0670187738890633e-05
Newton iteration  2
err =  1.710957520083692e-09
cost at iteration 595 :  5.193403215299303
Newton iteration  0
err =  0.0076969329675687145
Newton iteration  1
err =  2.09345680062105e-05
Newton iteration  2
err =  1.772012365723289e-09
cost at iteration 596 :  5.1851381433209305
Newton iteration  0
err =  0.007694038289053742
Newton iteration  1
err =  2.1198531828448528e-05
Newton iteration  2
err =  1.8333207316496864e-09
cost at iteration 597 :  5.176887927187565
Newton iteration  0
err =  0.007691243072542982
Newton iteration  1
err =  2.1460882818521602e-05
Newton iteration  2
err =  1.8944752194511004e-09
cost at iteration 598 :  5.168652560466924
Newton iteration  0
err =  0.007688546557093541
Newton iteration  1
err =  2.1720254346073952e-05
Newton iteration  2
err =  1.9549982557593685e-09
cost at iteration 599 :  5.160432037078894
Newton iteration  0
err =  0.007685947485827505
Newton iteration  1
err =  2.1975093732839245e-05
Newton iteration  2
err =  2.014337234123511e-09
=================Mesh Quality Check=================
max quality: 1.37841365416341  min quality: 0.5553651552984573  mean quality: 0.858318409251656
===================================================